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FOREWORD 


This document presents results of the study "Methods for Conparative 
Evaluation of Propulsion System Designs for Supersonic Aircraft." The NASA 
technical representative was Dr. Edward A. Willis. In addition to the authors 
noted, significant contributions to this study and report were made by 
Ellwood Bonner, aerodynamics ; Henry K. Chin and Louis C. Young, propulsion . 
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METHODS FOR COMPARATIVE EVALUATION 
OF PROPULSION SYSTEM DESIGNS FOR 
SUPERSONIC AIRCRAFT 


By Ray M. Tyson, Ronald Y. Mairs, Floyd D. Halferty, Jr., 
Bruce E. Moore, David Chaloff, and Arnold W. Knuds en 

Los Angeles Aircraft Division, Rockwell International 


INTRODUCTION 


The National Aeronautics and Space Administration is conducting a 
continuing program of advanced supersonic technology studies with the objective 
of developing an adequate technology base to seaport development of future 
supersonic cruising aircraft. It is recognized in this program that one of 
the more sensitive problems in the synthesis of a successful supersonic cruis- 
ing aircraft is that of airframe/engine integration. This process must inves- 
tigate and properly manage the interactions between the technical disciplines 
of external aerodynamics, internal aerodynamics, engine cycle design, acous- 
tics, mass properties, and structural design; and, it must be responsive to 
the practical considerations of fabrication, maintenance, and operation. 

The results of a recently completed study, reference 1, of the effects of 
nacelle size and nacelle shape on the drag, weight, and wing camber plane 
warping of a supersonic transport illustrated the sensitivity of these param- 
eters to relatively small changes in nacelle shape. The resultant shape of a 
nacelle is dependent on the geometry of the engine (inlet area, mounting pro- 
visions, accessory location, nozzle area, etc.) since this establishes certain 
control points in the design of the nacelle. It is important, therefore, that 
the engine designer be aware of this sensitivity to engine geometry, and be 
provided with some guidelines for favorable geometry relationships. It is 
probable that some engine geometry control can be achieved by the designer 
with no penalty in engine performance, although on a total system basis some 
engine performance degradation could be accepted in trade for reduced drag. 

Although considerable effort has been expended on the problem of airframe/ 
engine integration, it has been mostly in the nature of point designs. The 
study of reference 1 produced results for two specific nacelle shapes which 
resulted from installation of a dry turbojet engine and a duct heating turbo- 
fan engine. A comparison of these results shows the superiority of one nacelle 
shape over the other, but gives no information directly applicable to other 
installations having differing nacelle shapes. This report therefore treats 
nacelle shape and size in a parametric fashion so that a range of propulsion 
systems can be readily compared on a consistent basis. To meet the user's 
needs, it was clear that methodology faster and more convenient than the 
traditional aircraft-preliminary-design process would be required- -even at 
some cost in terms of accuracy. Therefore, the approach was taken of organiz- 
ing a relevant, existing set of nacelle drag data (reference 2) , together with 
supplementary data points as required to cover the parametric range, into a 
conputer table -lookup program. The program then yields supersonic wave and 
friction drag increments as function of size and shape parameters for a 



representative supersonic cruise airplane configuration (reference 3). The 
drag code, combined with linear sensitivity factors (derived from perturbation 
studies of the reference 2 airplane), provides the desired rapid approximate 
methodology for corrparing alternative propulsion system designs. 

The methods of analysis and major results of this study are described 
herein in the "STUDY PROCEDURE" section. User’s information for the code, 
program listings and mathematical details are presented in the Appendix. 

SUNNARY OF RESULTS 

The present work is an extension of a previous study performed for NASA 
Langley Research Center (contract NAS1-13906) and documented in reference 2. 

In that program, a baseline airplane was defined. Under the current contract, 
the baseline was revised slightly as described on page 71 for consistency in 
validating the approximate method; the revised baseline is used when perturbat- 
ing and comparing airplanes with other engines. The baseline and revised 
baseline airplanes are described in Table 1. The baseline airplane was based 
on the NASA modified SCAT 15F vehicle described in reference 3. Parametric 
data were generated showing the effects of variations of nacelle shape on 
cruise drag for a range of shapes that reasonably cover engine designs applicable 
to supersonic cruising aircraft. Generally, it was found that nacelles shaped 
such that the maximum cross-sectional area occurred at or near the nozzle exit 
and having little or no boattail resulted in the lowest wave drag. In fact, 
nacelle shapes were found that produce favorable interference effects (drag 
reduction) of such magnitude as to nearly offset the friction drag of the 
nacelle. These results are valid only for vehicles of this general configura- 
tion and nacelle location. Different vehicle configurations or nacelle 
locations could results in different "best" shapes. In considering possible 
trades of reduced drag through design changes in the engine for some penalty 
in engine weight and specific fuel consumption (SFC) , it is necessary to have 
visibility of the net impact of all three effects on the total airplane in 
order to make a comparative evaluation. Therefore, sensitivity data were 
developed for the effects of changes in drag, propulsion system weight, takeoff 
thrust, and SFC on the takeoff gross weight as a figure of merit. Results of 
the weight sensitivity trades showed that the airplane gross weight is highly 


TABLE 1.- BASELINE AIRPLANES 



NASA Langley Study 

REVISED BASELINE 

Design Mission F&nge, km (n mi) 

7 408 

(4 000) 

7 408 

( 4 - 000) 

Design Cruise Mach Nunber 

2.4 

2.4 

2.4 

2.4 

Payload (292 passengers), kg (lb) 

27 682 

(61 028) 

27 682 

(61 028) 

Balanced Field Length, m (ft) 

3 190 

(10 500) 

3 190 

(10 500) 

Engines (41 

VSCE 502B 

VSCE 502 B 

VSCE 502B 

VSCE 502B 

Takeoff Gross Weight, kg (lb) 

322 046 

(712 188) 

316 783 

(698 375) 
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sensitive to both drag and engine SFC at supersonic cruise. A cne-drag-count 
change (approximately 1 percent of airplane drag) results in a 1-percent 
takeoff gross weight change; a 1 -percent change is SFC also results in a 1- 
percent change in takeoff gross weight. Changes in drag or SFC at other flight 
conditions and changes in propulsion system weight had relatively small effects 
on takeoff gross weight. 

The follow-on program, described in this report, was intended to render 
the above-mentioned parametric data into a convenient, useable form. The 
objective was to develop a reasonably accurate method for the rapid, prelimin- 
ary evaluation of the effects of variations in propulsion system design 
parameters on the total system performance of an integrated engine/airframe 
system. The figure of merit used was the airplane takeoff gross weight to 
perform a design reference mission. The effort was organized around the 
following five tasks: 

(1) Estimation of supersonic cruise drag increments reflecting nacelle 
shape and size (in the form of a computer table look-up program) 

(2) Estimation of propulsion system installation weight 

(3) Estimation of airplane takeoff gross weight 

(4) Validation of the approximate method 

(5) Reporting 


Estimation of Supersonic Cruise Drag 

A computer table look-up program was developed (see appendix) which yields 
the incremental wave and friction drags of nacelles as functions of five 
nacelle geometry variables and airplane mach nunber. The drag increments 
are for the total vehicle relative to the vehicle with nacelles removed. The 
five nacelle shape parameters used as inputs to the program are: 

A c Inlet capture area 

Aj^ Nacelle maximun cross-sectional area 

A^ Nozzle exit area (supersonic cruise position) 

Distance from inlet cowl leading edge to maximum cross-sectional 
area 

L Nacelle total length 

^REF ^ e ^ erence w i n g area 

It has been found that the table look-up results correlate best with more 
detailed analyses when the maximum cross-sectional area and its position are 
based on the area that occurs at the intersection of straight lines originating 
from the inlet and nozzle and whose slopes nearly match the slopes of the 
actual nacelle. A sample output from the computer program is shown in table 2. 
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TABLE 2 


SAMPLE OUTPUT 






NACELLE * tCMETKY 


CAP TUP f: AREA 

2.79 

SO 

M 

MAXIMUM AKl A 

4. IF 

SC 

M 

Nvzztr ah' 

3,*£ 

so 

M 

LOC. OF MAX. AF c A 

l .77 

M 


TOTAL LtM-TH 

K »3t> 

M 



3 D.oe - 50 -rn~ 

45.00 SO FT I 
37.5-' SO FT) 
20. 4C FT ) 

33.94 FT ) 


WING REFERENCE AREA 920.'3 SC H 


1 0000.00 SC TTT 


INCREMENTAL NACELLE t F AG COEFFICIENTS 
MACH 1.? CDF - C.y0«.*'3 

MAth *.?; 

M*f.H ].2‘. 


XMAX/L* .0;i. 


AN/4C= 1.2 S' 


COWa 0.0004C 
CUW* C.CCCC7 


CDC* 0.00043 
CIO* C.CC054 
CDO* 0.00093 


AMAX/AC* 1.5GC L/OC* 5 


INCREMENTAL NaC.ILU I R AG COEFFICIENTS - LTV REFERENCE VEHICLE 

mac m i.-> u.f* c.:;o7; tun * -i.c:-i>4i coo* 0.00031 

MAC 1 .'.32 C.t»F «• C.CCOl.C CDW* -C.00C1F CDO* 0.~TCC*? 

MAE H COO* C. 00032 
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This task also included preparation of design guidelines for the 
installation of the engine in the nacelle so that consistency in definition 
of design control points for external nacelle shaping is achieved. Guide- 
lines for establishing inlet and inlet cowl shape, engine envelope definition 
(including provisions for wiring, pitching, power takeoff, engine accessories , 
aircraft accessories, fluid reservoirs, air bleed ducts, and engine moults) , 
structural allowances, and engine cowl and nozzle fairing shapes were defined. 
A saiple of the nacelle shape buildup is shown in figure 1. 


Estimation of Propulsion System Installation Weight 

A simplified procedure was developed for the prediction of nacelle struc- 
ture weight. Weight estimation of aircraft structure is a complex process and 
requires more design detail than will ordinai .lv be performed in the type of 
preliminary studies being cons i de red here; therefore, the procedure was keyed 
to gross elements of the propulsion system installation ami yields only 
approximate weights. The inportant aspect of having a well-defined procedure, 
even though considerable tolerance in the results must be accepted, is that 
consistency is achieved in making conparative analyses. 


Estimation of Airplane Takeoff Gross Weight 

A method has been defined for the determination of the inpact of the 
propulsion system installation (cruise drag, SFC, weight) on the total system 
performance utilizing takeoff gross weight as the figure of merit. This 
method was based' on results of sensitivity studies performed for Langley 
Research Center. Utilizing these sensitivity values, the drag and weight 
increments from tasks 1 and 2, and SFC’s from engine performance estimates, 
this procedure yields the airplane takeoff weight required to accoopfish 
the design mission. The baseline vehicle is the vehicle defined in 
reference 2 with the VSCH 502B engine. The total change in vehicle 
takeoff gross weight due to propulsion changes may be determined from 
the equation: 


TCX> new TOGW baseline X ^FC X R C D X R WT X R FNE TO 


where f, R" factors are the relative takeoff gross weight factors for each of 
the propulsion changes obtained from a linear sensitivity analysis of the 
baseline system. The drag factor is based on a supersonic cruise increment. 
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Legend 

(T) Engine accessories. Encapsulated for 
cooling. Includes all tenperature - 
limited components. 

• 

( 2 ) 'jigine lube reservoii 

Engine peripheral hardware. 

Includes: engine fluid lines (anti-ice air, fuel, 
lube oil. hydraulic, drains, etc), 
variable geometry mechanisms, 
electrical harnesses, 
instnmientation. 

Local protrusions will occur beyond this 
envelope. 

( 4 ) Conpressor bleed manifold. 

(5) Engine power takeoff for aircraft 
accessories drive. Angle gearbox and 


power transmission shaft to wing - 
mounted accessories drive gearbox. 



(6) Engine and nacelle support structure 
on wing, inside pylon. 

© Engine cowl. Nonstructural, hinged 
from pylon. Includes small doors 
for local access. 

(8) Main mount (front) 

(5) Stabilizer mount (rear) 

Dimension in meters (inches) 
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Validation of the Approximate Method 


Standard preliminary design procedures were applied in the installation 
of four representative engines selected from the NASA Supersonic Cruise 
Aircraft Research (SCAR) program engine studies, in the baseline supersonic 
transport airplane. Drag and weight estimates were made utilizing conventional 
procedures. The airplanes were then sized to the design mission utilizing an 
automated reiterative process. The results of this task provide a more exact 
evaluation of the selected engines than is obtainable with the approximate 
method and thus serve as a reference for its evaluation. 

The four engines considered were the refined (January 1S76) Pratt and 
Whitney Aircraft (PWA) VSCE 5Q2B and VCE 112C, and the General Electric Company 
(GE) GE21/J10 B1 and GE21/J11 B3. These engines were chosen as representative 
exanples which would exercise the approximate method over a sizeable range, 
in order to determine its limits of validity. This can be done in a consistent 
fashion by comparing the approximate and detailed results shown in Table 3 
for each engine . (On the other hand, meaningful engine -to -engine comparisons 
cannot be nude on the basis of table 3 because the example engines do not 
necessarily reflect a consistent set of basic technology assumptions, noise 
characteristic-, or state of evolution within the SCAR program. Hence, these 
and similar results discussed later in this report should not be interpreted 
as being indicative of the final outcome of the ongoing SCAR engine studies.) 


TABLE 3. - COMPARISON OF APPROXIMATE AND DETAILED VEHICLE TAKEOFF GROSS WEIGHTS 


Engine 

Weight Based on 
Detailed Analysis 
kg (lb) 

Weight Based on 
Sensitivities 
kg (lb) 

VSCE 502B 

316 783 (698 37S) 

Revised Baseline 

VSCE 502B (refined) 

320 146 (70S 790) 

320 046 (705 568) 

VCE 112C 

402 625 (887 622) 

401 092 (884 258) 

GE21/J1Q B1 

514 450 (1 134 149) 

463 708 (1 022 283) 

GE21/J11 B3 

629 306 (1 387 359) 

510 136 (1 124 638) 


By comparing detailed and sensitivity results for each engine in table 3, 
it is immediately clear that a good level of agreement has been reached. The 
relative error (normalized by the gross weight based on detailed results) 
is shown in Figure 2 as a function of the total incremental change in TOGW 
(normalized by the baseline value). As might be expected from theoretical 
considerations, the error is negligible for small perturbation; in fact, 
it does not exceed 2 percent of the takeoff gross weight until the increment 
itself is in excess of 30 percent. It is inportant to note also that the 
error is consistent, i.e. always of the same sign (the approximate method 
underpredicts). Thus, even among highly-dissimilar engines, the correct 
ranking is preserved. With these facts in mind, it is concluded that the 
approximate method is in fact a reliable and reasonably accurate tool for 
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such purposes as engine evaluation and comparison, over a range of about 
♦ 30 percent from the baseline TOGK. Considerably larger increments also 
could be accepted (temporarily) as intermediate steps in an optimization 
study, provided that the final case of interest is within the +_ 30 percent 
band. 


The user should nevertheless observe several cautions in applying these 
results. As a general practice, it is desirable to check the "final result" 
of a study by detailed methods. This is strongly recommended for cases 
approaching or passing beyond the accuracy band. The sensitivity values 
("R" factors) are to some extent dependent upon the engines sizing criteria, 
the assumed mission profile and flight rules. The user should therefore 
review these items carefully before beginning a study and generate a more 
appropriate set of "R" factors if significant differences are noted. More 
fundamentally, it should be recognized that the wave drag data is strictly 
applicable only to the reference 3 airplane configuration and geometrically 
similar scaled versions thereof. Trend results with nacelle shape for 
different airplanes of the same general arrangement are believed to be 
represent at ire, although detailed agreement would not be expected. The use 
of the present data for airplanes having significantly different shape, 
proportions or nacelle treatment is not indicated. Doubtful cases should 
be checked at several points to validate the data and/or establish correct ions. 



1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 


Relative takeoff gross weight calculated by approximate method 
Figure 2.- Error characteristics of approximate method 
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Conclusions and Recommendations 

Engine shape, airflow lapse rate with mach number , thrust lapse rate 
with mach number, SFC and noise characteristics have large impacts on vehicle 
takeoff gross weight. As an example of engine shape effects, a comparison 
of cross sectional area variation of nacelles with the VSCE 502B and VCE 
112C engines is shorn in figure 3. The only significant difference in shape 
is that the VCE 112C has a smaller nozzle exit area. This results in drag 
and takeoff gross weight changes as shown in table 4. Thus, a nozzle that 
is 0.043 b (1.7 in) smaller results in takeoff gross weight increment of 
3 200 kg (7 100 lb) just due to the nacelle drag change. 

Engine airflow lapse rate with mach nuaber directly affects inlet 
capture area. For example, the GE21/J10 B1 has approximately 16-percent 
lower supersonic cruise airflow relative to takeoff airflow than does the 
VSCE 502B. The smaller capture area results in approximately 4-percent 
lower inlet recovery at static conditions and therefore reduced takeoff thrust. 
In addition, the smaller capture area results in a more rapid increase of 
nacelle cross-sectional area with nacelle length, and therefore higher 
drag. 


Engine takeoff thrust and thrust lapse rate with mach number have 
significant effects on engine size required to meet takeoff distance require- 
ments. For exanple, the VCE 112C has 6-percent lower takeoff thrust at 
static conditions and 20-percent lower thrust at mach 0.3 (at reduced power 
to meet noise requirements) than the VSCE 502B for a given static takeoff 
airflow. This resulted in am increase in engine size of approximately 15 
percent to meet balanced field length requirements. 

A change of 1 percent in SFC at supersonic cruise results in a 1-percent 
change in vehicle takeoff gross weight or about 3 200 kg (7 100 lb) . 

Engine exhaust noise characteristics have a significant impact on 
vehicle takeoff gross weight. All four engines were assumed to employ thrust 
cutback at the takeoff noise measurement point. However, all the engines 
did not take full advantage of the extra ground attenuation while the air- 
craft was still cm the ground. For exanple, the GE21/J11 B3 has 26-percent 
lower thrust-per-unit airflow than the VSCE 502B. Thus, the GE21/J11 B3 
yields approximately 6 db lower sideline noise at takeoff on the ground, 
but it must be sized larger to meet the takeoff distance requirement. 

The sensitivity method has been shown to be a valid method for prelimin- 
ary assessment of propulsion system modifications, and it is therefore 
recommended to be used for this purpose. Continued airframe/propulsion 
integration studies and coordination effort between engine and airframe 
manufacturers in the aforementioned high-sensitivity areas are also recomnended. 
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A 

c 


Nacelle characteristics 

PWA VSCE 502B 

PWA VCE 112C 

A^ 2.93 sq m (31.5 sq ft) 

2.93 sq m (31.5 sq ft) 

1 10.3 m (33.25 ft) 

11.2 m (36.86 ft) 

Ag 0.17 sq m (0.57 sq ft) 

0.15 sq m (0.50 sq ft) 



Figure 3. - Normalized nacelle cross-sectional area variation 
for VSCE 502B and VCE 112C. 






TABLE 4. - DRAG AND TAKEOFF GROSS WEIGHT INCRE4EMTS DUE TO 
CHANGES IN ENGINE SHAPE 



ACL, Supersonic Cruise 
Nacelle Drag 

Takeoff Gross Weight 

Engine 

Increment Relative 

Increment Relative 

to VSCE 502B 


to Nacelles Off 

kg (lb) 

VSCE 502B 

0.00046 

Base 

VCE 11 2C 

0.00055 

3200 (7100) 


SYMBOLS 


A Area, sq m (sq ft or sq in) 

BLB Boundary layer bleed 

BLC Boundary layer control 

BP Basepoint 

C Coefficient or Chord, m (ft or in) 

d Diameter, m (ft or in) 

D Drag, daN (lb) 

db Decibel 

F Thrust, kg (lb) 

K Drag-due-to-lift factor 

I Length, m (ft or in) 

L Lift, N (lb) 

M Mach number 

R Relative TOGW factor 

S Area, sq m (sq ft or sq in) 

SFC Specific fuel consumption, kg/hr/N (lb/hr/lb) 

T Thrust, N (lb) 

TOGf Takeoff gross weight, kg (lb) 

V Velocity, m/sec (ft/sec) 

W Weight, kg (lb) 

X Nacelle station, m (ft or in) 

A Increment 

Subscripts 

AMAX Maximum cross-sectional area 

B Base 
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c Capture 

CD Drag coefficient 

D Drag 

. F Friction 

f Fuel 

i Inlet throat 

K Indicates lift coefficient at minimum drag 

L Lift 

LO Loftoff 

MAX Maximum 

n Nozzle exit 

ne Net effort 

P Profile 

REF Reference 

R toot 

SFC Specific fuel consumption 

SUB Subsonic 

SUPER Supersonic 

TO Takeoff 

W Wave 

WT Weight 

0 Freestream 

1 Critical engine failure 

STUDY PROCEDURE 
Approach 

The general approach of this study included using the baseline airplane, 
parametric nacelle drag results, and takeoff gross weight sensitivities devel- 
oped in the NASA Langley Research Center contract of reference 2. A nacelle 
drag table look-up conputer program and guidelines for determining nacelle 
shape were developed to allow estimation of supersonic cruise drag. A method 
to assess the propulsion system installation weight was defined. A method of 
determining vehicle takeoff gross weight using vehicle sensitivities to propul- 
sion changes was developed. The method was verified by analyzing in detail 
four selected propulsion systems. 

Because of the dependence of this study on the baseline airplane and 
ground rules of the study of reference 2, the definition of the baseline air- 
plane is included. In this report, descriptions of the airplane configurations 
used are as follows: 

(1) The reference airplane is the NASA-modified SCAT 15F arrow wing 
supersonic transport (defined in reference 3) , 
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(2) The basepoint vehicle is the reference modified only as required to 
install the Pratt and Whitney Aircraft (PWA) VSCE 50 2B engine, 

(3) The baseline airplane is the basepoint resized to the design require 
ments on a standard-plus- 8° C day. 

The structure design and operational eupty weight of the reference 
airplane were assigned to meet all design criteria. Weight and aerodynamic 
characteristics of the study airplanes were derived by increments from the 
reference configuration. 


Baseline Airplane Definition 


Basepoint airplane .- The "basepoint" airplane for this study is based on the 
NASA modified SCAT 15F arrow wing reference configuration as described in 
reference 3. The propulsion system of this airplane has been replaced with 
PWA variable stream control engines (VSCE 502B) having 408 kg/sec (900 lb/sec) 
airflow each and witn axisymnetric variable geometry inlets designed for 
mach 2.4 cruise conditions. The resuiting basepoint vehicle is shown in fig- 
ure 4. This airplane has a gross weight of 336 973 kg (742 890 lb), a range 
of 7471 km, (4034 n.mi.), and a balanced field length of 3017 m (9898 ft). 

All performance and sizing calculations were made using the Rockwell 
Vehicle Sizing and Performance Evaluation Program (VSPEP) . This conputer 
program is a design tool capable of scaling a known basepoint vehicle accord- 
ing to specified values of several different design parameters. These include 
vehicle gross weight (or fuel weight) , thrust-to-weight ratio (or engine size) , 
wing- loading (or wing area), and payload or fixed equipment weight and volume. 
Performance may be determined at specified gross weight, or alternatively, a 
search routine permits automatic sizing of the vehicle gross weight such that 
a specified radius or range of the design mission is satisfied. Vehicle per- 
formance is calculated internally from a set of subroutines progranmed accord- 
ing to a detailed performance analysis model. The subroutines are general in 
nature and permit calculation of a wide variety of mission profiles. Several 
mission profiles may be calculated simultaneously. Takeoff and landing dis- 
tances and maneuvering capability may also be determined. Figure 5 illustrates 
the evaluation process. 

Typical mission legs which may be calculated include warmup, taxi, takeoff, 
climb, descent, cruise, and loiter operations. Climb and descent performance 
are determined by nunerical integration of the equations of motion along a 
specified flight schedule. Internally generated schedules are also available, 
including minimum time and minimum fuel flight paths as defined by the energy 
method. Constraints on the allowable flight regime are included. Cruises and 
loiters may be determined at fixed or optimum speeds and altitudes. Numerical 
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Figure 4. - Basepoint airplane 


FOLDOUT FRAMfeJ 1 ? 





PRECEDING PAGE BLANK NOT FILMED 



\o 



searches axe used to determine optimal speeds and altitudes at the beginning 
and end of each of these legs. 

Data input to the VSPEP for the AST basepoint vehicle include: 

* heights broken down by major component , along with scaling 
information on the wing, tails, fuselage, and engines . 

* Drags broken down by major coqxment and by type (e.g. , friction 
drag, wave drag, drag due to lift, base drag). 

* Installed propulsion dat , including thrust and fuel flow as 
functions of speed, altitude, and power setting. 

* Dimensional data such as lengths, areas, and volumes for major 
coaponents and the total vehicle. 

Performance items calculated by the VSPEP on the basepoint and baseline 
vehicles for this study consist of the following: 

(1) Design mission range 

(2) Alternate mission range 

,3) Takeoff distance with FAR 56 (Federal Aviation Regulation, 
part 36) noise requirements 

(4) Balanced field takeoff distance 

(5) Thrust -to -drag ratio at mach 2.32, 18 300 m (60 000 ft) 

(6) Thrust- to -drag ratio at mach 1.2 during the clinb leg 

A description of each of these performance items is given in the follow- 
ing paragraphs. Because engine data were provided for a standard-plus -8°C 
(14.4°F) day, all airplane performance characteristics were computed for that 
atmospheric condition. 

Design mission .- A profile of the design mission is shown in figure 6. This 
mission consists mainly of a mach 2.32 cruise. Fuel reserves as reconmended 
in reference 4 are calculated for an a’temate airport located 460 km (250 n.mi.) 
from the destination airport. 

The design mission consists of: 

(1) Warmup and takeoff - 10 minutes at die pt „r plus 1 minute at 
maximum power. 
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Figure 6. -Reference mission 



(2) Clint) - Maximum power clint> and accelerate to cruise altitude and 
mach nimber. 

(3) Cruise - Cruise at mach 2.32 at altitude for best cruise range. 

(4) Descent - Descend and decelerate to aach 0.5 aid 457 a (1500 ft) 
using idle power. 

(5) Approach and land - Descend to mach 0.3 at sea level using 
idle power. 

(6) Taxi - 5 minutes at idle power. 

(7) Reserve allowance - 5-percent of total fuel used in all previous legs. 

(8) Reserve clint) - Climb to subsonic cruise conditions . 

(9) Reserve cruise - Subsonic cruise at mach n uaber and altitude for 
best range. 

(10) Reserve descent - Descend and decelerate to holding altitude and 
mach nunber using idle power. 

(11) Reserve hold - Loiter for 30 minutes at 3048 m (10 000 ft) at the 
mach nunber for best endurance . 

(12) Reserve approach and land - Descend to sea level using idle power. 

Alternate mission .- A profile of the alternate mission is shown in figure 7. 

The first half of the alternate is identical to the first half of the design 
mission. At the point corresponding to the nCipoint of the design mission, a 
failure is assumed to occur ir. the most critical engine. At this point, the 
airplane descends and continues to cruise subsonically with one engine wind- 
milling. The fuel reserve remaining at the aid of this mission is equal to the 
reserve fuel as calculated for the design mission. 

The alternate mission consists of: 

(1) Warmup and takeoff - Same as design mission. 

(2) Climb - Same as design mission. 

(3) Cruise - Same as design mission. 
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(4) Descent - Descend and decelerate to subsonic cruise conditions using 
idle power, following failure of nst critical engine. 

(5) Cruise - Subsonic cruise at mach number and altitude for best range 
with one engine inoperative. 

(6) Descend and land - Descend to sea. level using idle power. 

(7) Reserve - Allow total reset - . ? fuel equal to that calculated for design 
Mission legs 7 through 12. 

Balanced field takeoff. - Takeoff distance is calculated over a 10.7 m (35 ft) 
obstacle. It is assmed that a maximum usable lift coefficient of 0.555 is 
available for climbcut. Balanced field length involves three requirements : 

(1) Distance for a normal takeoff is calculated with all engines 
(throttled if necessary so that FAR 36 noise requirements are not exceeded) and 
this distance is multiplied by 1.15. 

(2) Distance is calculated for a takeoff when an engine fails at the 
critical speed and the airplane continues the takeoff. In this instance, the 
throttles may be advanced after the engine failure if they are not already at 
maxi bib power (without regard to noise requirements) . 

(3) Distance is calculated when an engine fails at the critical speed and 
the takeoff is aborted. The airplane is stopped by applying brakes and cutting 
the remaining engines to idle. The critical speed for engine failure is deter- 
mined by varying the speed at which engine failure occurs (i.e. , Vj) until the 
accelerate -continue distance is equal to the accelerate- stop distance (i.e., 
segments B + C = D + Eas shown in figure 7). The balanced field length is 
then defined as the greatest of items (1), (2), and (3). 

Thrust- to-drag ratio .- The thrust -to- drag (T/D) ratio is calculated using max- 
inun available thrust at 2.32 mach, 18 300 m (60 000 ft). Drag is that for 
level flight at the same conditions. Airplane weight is that at the start of 
the supersonic cruise as calculated for the design mission. The thrust-to-drag 
ratio is also calculated for the point in the cl irtt -accelerate leg at which 
mach 1.2 is reached. In this case the altitude and vehicle weight are the 
actual values during the cliob at which the vehicle teaches mach 1.2. 

Baseline airplane .- Hie "baseline" airplane for this study is a resized version 
of the aforementioned "basepoint." Resizing was accomplished by exercising 
the VSPEP for a matrix of thrust -to-weight and wing loading values, 
and allowing the program to search for the gross weight, in each case, that 
satisfies the design mission range requirement of 7408 km (4000 n.mi.). Plots 
of the results are shown in figures 9 through 11. The parameters shown include 
vehicle gross weight as well as those performance items for which requirements 
must be met. 
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The balanced field length requirement is plotted on the airplane gross 
weight plot in figure 9. This allows a "baseline" airplane to be chosen which 
is defined as the miniatm gross weight vehicle that meets or exceeds the 
following pe rfo r m ance requirements : 

Design mission range 7408 km (4000 n. mi.) with 292 passengers 

Balanced field length 3200 m (10 500 ft) 

Minimus T/D during cliab 1.2 
or cruise 



A - Distance tp to critical engine failure 

B - 3- engine acceleration distance from Vj to 

C - 3-engine lift-off to barrier distance 

D - Distance gained after engine failure before full 

brake application 

E - Stopping distance 

Vj - Critical engine failure speed 

- Lift-off velocity 

Figure 8. - Balanced field length definition 
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Figure 9. - Gross wight versus thrust -to-weight and wing loading. 
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Figure 10. -Balanced field length versus thrust- to-weight and wing loading 
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Since the thrust- to-drag requirements are exceeded for all cases, 
figure 10, only the balanced field length requirements are included in figure 
9. 


The resulting "baseline" airplane has a gross weight of 323 046 kg 
(712 188 lb); a thrust-to-weight ratio of 3.16 n/kg (0.323 lb/lb) based on 
installed, static takeoff thrust; and a wing loading of 354 kg/sq m (72.5 lb/ 
sq ft) based on gross wing area. (As discussed on page 71, the baseline was 
revised slightly for consistency in validating the approximate methods, but 
the disc jss ion presented here is for the NASA Langley stud)' baseline.) 

Further airplane design and performance characteristics for both the "base- 
point" and "baseline" airplanes are shown in table 5 . Design and alternate 
mission sumnaries are shown in tables 6 through 9 for the baseline airplane. 

In these tables, the first leg of the alternate mission includes the first four 
legs of the design mission while leg five of the alternate includes reserves 
for legs nine tni'. igh 15 of the design mission. The path followed during the 
climb -accelerate leg is a mini nun fuel path calculated internally by the VSPEP 
program. This path as calculated for the baseline airplane is shown in 
figure 12. 

Propulsion .- Because many of the current and recently completed supersonic 
cruising aircraft studies have used axisymmetric inlets, a mixed-conpression, 
axisyumetric inlet was defined for use in the basepoint aircraft for this 
study. The inlet diameter is 1.93 meters (76 inches), and capture area is 
2.926 square meters (4 536 square inches). For takeoff, the basic centerbody 
is held in the transonic position, but the fore and aft conical segments are 
translated aft to create a centerbody auxiliary inlet. The auxiliary inlet 
opening in the centerbody is 10 percent of capture area. Inlet pressure recov- 
eries and spillage, bypass, and BLC drags were estimated, and their effects 
were included in installed propulsion performance. Engine accessories were 
assumed to be located in the wing. The nacelle drawing is shown in figure 13. 

The engine performance data available for the VSCE 502B engine included 
the effects of an inlet recovery schedule, nozzle external drags (base plus 
boattail) , 0.45 kilogram-per-second (1.0 pound-per- second) high-pressure com- 
pressor air bleed, and 149 kilowatts (200 horsepower) power extraction. 

Installed performance data were computed by modifying the engine data to 
include the effects of changes in inlet pressure recovery and inlet drags 
(spillage, bypass, and boundary layer control). Because the amount of engine 
data available was not sufficient to compute aircraft mission performance, 
additional installed performance data were generated by calculating corrected 
thrust and fuel flow parameters and by extrapolating based on trends of engines 
with similar characteristics. Fortunately, these techniques were required only 
at flight conditions where the airplane flies for a short duration. Thus, any 
possible errors due to data extrapolation should have minimal effect on air- 
plane performance. All data were for standard-plus-8°C day. 
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TABLE 5. -AIRPLANE CHARACTERISTICS 



Basepoint 

Baseline 

Engine Airflow, kg. /sec. (lb. /sec.) 

408 

(900) 

395 

(869) 

Thrust -to -weight, n./kg. (lb. /lb.) 

3.14 

(0.3208) 

3.16 

(0.323) 

Reference Wing Area, sq.m. (sq. ft.) 

926 

(9 969) 

827 

(8 902) 

Gross Wing Area, sq.m. (sq. ft.) 

1 022 

(10 996) 

912 

(9 819) 

Wing Loading (gross area), kg. /sq.m, (lb./sq. ft.) 

330 

(67.56) 

354 

(72.53) 

TOGW, kg. (lb.) 

336 973 

(742 890) 

323 042 

(712 188) 

Fuel Weiglt, kg. (lb.) 

158 475 

(349 834) 

151 643 

(334 754) 

Max. Wing Fuel, kg. (lb.) 

207 950 

(459 050) 

175 484 

(387 382) 

Design Range, km. (n.mi.) 

7 475 

(4 034) 

7 412 

(4 000) 

Eng. CXit Range, km. (n.mi.) 

6 285 

(3 392) 

6 259 

(3 378) 

FAR 36 T.O. Dist., m. (ft.) 

2 498 

(8 194) 

2 648 

(8 661) 

Bal. Field T.O. Dist. , m. (ft.) 

3 017 

(9 896) 

3 190 

(10 466) 

Thrust -to- Drag e 2.32 M/18 300 m. (60 000 ft.) 


1.747 


1.716 

Thrust-to-Drag 8 1.2M/Climb 


2.005 


2.019 

Initial Cruise L/D 


9.649 


9.541 

Initial Cruise SFC (installed) kg./hr./daN. 

(lb. /hr. /lb.) 

1.393 

(1.366) 

1.385. 

(1.358) 





UNITS 


TIME, 


min. 

■CTE9KM 

10.0 

10.0 

0 

0 

1.3 

11.3 

10 

10 

12.8 

24.1 

288 

299 

82.9 

107.0 

3 405 

3 704 

82.6 

189.7 

3 395 

7 100 

17.3 

207.0 

296 

7 396 

1.5 

208.6 

12 

7 409 

mm 

213.6 

0 

7 409 

0.0 

213.6 

0 

7 409 

0.7 

214.3 

5 

7 414 

10.2 

224.5 

158 

7 573 

10.8 

235.3 

181 

7 754 

9.5 

244.9 

123 

7 877 

30.0 

274.9 

0 

7 877 

3.7 

278.6 

33 

7 911 

3 












TABLE 7. - BASELINE DESIGN MISSION SUMMARY - ENGLISH UNITS 


LEG. NO. OPERATION 

WEIGHT, 

lbs. 

ALTITUDE, 

ft. 

MACH NO. 

FUEL USED, 
lbs. 

TIME 

min. 

TOTAL TIME 
min. 

RANGE 

n.m. 

TOTAL 
RANGE n.m. 

1 

INITIAL WEIGHT 
WO 8 TO 

712 188 
703 653 

0 

0.305 

8 

534 

10.0 

10.0 

0 

0 

2 

CL TO 1500 

699 022 

1 500 

0.500 

4 

631 

1.3 

11.3 

5 

5 

3 

CLB-ACC 

644 931 

54 943 

2.320 

54 

090 

12.8 

?A.\ 

155 

161 

4 

CRUISE 

528 936 

58 761 

2.320 

115 

994 

82.9 

107.0 

1 838 

2 000 

5 

CRUISE 

432 684 

62 889 

2.320 

96 

252 

82.6 

189.7 

1 833 

3 833 

6 

DESCEND 

429 318 

1 500 

0.500 

3 

365 

17.3 

207.0 

159 

3 993 

7 

DES-LAND 

428 800 

0 

0.300 


518 

1.5 

208.6 

7 

4 000 

8 

TAXI -All 

427 374 

0 

0.0 

1 

425 

5.0 

213.6 

0 

4 000 

9 

SPCT ALL 

413 134 

0 

0.0 

14 

240 

0.0 

213.6 

0 

4 000 

IP 

CL TO 1500 

411 882 

1 500 

0.500 

1 

251 

0.7 

214.3 

2 

4 003 

11 

CLB-ACC 

400 052 

40 476 

0.950 

11 

830 

10.2 

224.5 

85 

4 089 

12 

CRUISE 

394 775 

40 773 

0.950 

5 

276 

10.8 

235.3 

98 

4 187 

13 

DESCEND 

392 853 

10 000 

0.470 

1 

921 

9.5 

244.9 

66 

4 253 

14 

LOITER 

378 632 

10 000 

0.454 

14 

221 

30.0 

274.9 

0 

4 253 

15 

DES-LAND 

377 448 

0 



0.300 

1 

184 

3.7 

278.6 

18 

4 271 

L 


TOTAL FUEL USED 

- 334 739 










TABLE 8. - BASELINE ALTERNATE MISSION SUNMARY - INTERNATIONAL UNITS 


:.G. NO. OPERATION WEIGHT, 

kg. 


INITIAL WEIGHT 323 046 

1 DES LEGS 1-4 239 921 

2 DES-DEC 239 539 

3 CRUISE 194 718 

4 DESCEND 193 850 

5 RESERVE 171 204 


ALTITUDE, MACH NO. 
m 


0 

8 212 
9 461 
0 
0 


.0 

.900 

.900 

0.300 

0.0 


FUEL USED, TIME TOTAL TIME RANGE TOTAL RANGE 
kg. min. min. km. km. 


83 121 107.0 107.0 3 704 3 704 

381 7.9 115.0 189 3 894 

44 821 136.7 251.7 2 248 '6 142 

868 10.2 262.0 113 6 256 

22 646 65.0 327.0 502 6 758 


Total Fuel Used = 151 838 


TABLE 9. - BASELINE ALTERNATE MISSION SUNMARY - ENGLISH UNITS 


LEG NO. OPERATION 

WEIGHT, 

lbs. 

INITIAL WEIGHT 

712 188 

1 DES LEGS 1-4 

528 936 

2 DES-DEC 

528 095 

3 CRUISE 

429 281 

4 DESCEND 

427 368 

5 RESERVE 

377 441 


ALTITUDE, MACH NO. fUEL USED TIME TOTAL TIME RANGE TOTAL RANGE 
ft. lbs. min. min. n.m. n.m. 


0 

26 944 
31 040 
0 
0 


183 251 
841 
98 813 
1 913 
49 926 


107.0 

7.9 

136.7 

10.2 

65.0 


107. 
115. 
251. 
262. 
327.0 


2 000 
102 
1 213 
61 
271 


2 000 
2 102 
3 316 
3 378 
3 649 


Total Fuel Used - 334 746 
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Figure 13. • Basepoint nacelle. 



At takeoff, _.*e engine power setting was scheduled so that the aircraft 
meets FAR 36-traded noise levels. The SAE exhaust jet noise prediction method 
was used with the modification to overall sound pressure level recommended by 
Bushell (reference 5). This modification has no effect on perceived noise 
level at static conditions, but it results in approximately 4 decibels (db) 
higher noise than the standard SAE method at mach 0.3. A 1.5 db reduction in 
sideline noise was assumed due to sideline shielding while the airplane is on 
the grcuid. An 8-decibel reduction in noise level due to the uuannular nozzle 
effect was assumed for all flight conditions and power settings. Information 
from PKA indicates that coannular configurations reduce noise by 7 to 9 decibels 
when the difference between core velocity and bypass velocity is 152 meters per 
second (500 feet per second) or more, with the core stream having the lower 
velocity. 

Mass properties .- The basepoint vehicle weight summary is given in table 10. 

The NASA reference vehicle weight sumaary (reference 3) from which the base- 
point was derived is also shown. The differences between the weights of the 
two vehicles are in the engines and nacelles. The basepoint vehicle has 
VSCE 502B 408 kg/ sec (900 lb/sec) airflow engines in lieu of the 363 kg/sec 
(800 Ib/sec) engines in the NASA reference vehicle. 

The VSCE 502B bare engine weight including nozzle and thrust reverser was 
supplied by PWA. Weight increments of 22.7 kg (50 lb) for residual fluids and 
22.7 kg (50 lb) for miscellaneous engine/airfrarae interfacing provisions were 
added to the bare weight to obtain an installed weight. Table 11 shows the 
installed engine weight suranary. 

The basepoint nacelle weight estimate is based on the nacelle drawing, 
figure 13. For the weight evaluation, the nacelle was divided into three 
sections: forward of the engine front face (inlet cowl), aft of the front face 

(.engine cowl), and inlet spike. The engine cowl weight was estimated at 
34.2 kg/sq m (7 lb/sq ft) of wetted area. This weight includes all the nacelle 
structure that supports and surrounds the engine and was derived from prior 
Rockwell International studies of a similar type. The inlet cowl and spike 
weights were calculated using statistical weight estimating equations obtained 
from the technical report SEG-TR-67-1, Preliminary design Methodology for Air- 
Induction Systems (reference 6). Engine mount weights were calculated sta- 
tistically at 1.5-percent of the engine weight. The mount weights are included 
with the nacelle weight. 

The weight summary of the basepoint nacelle is presented in table 12. 
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TABLE 10. - VEHICLE WEIGHT SUMARY 


ITB! 

! NASA REFERENCE VBHCLE 

f" BASEPOINT I 

Kg 

LB 

Kg 

LB 

Wing 

37 805 

83 347 

37 805 

83 347 

Horizontal Tail 

2 391 

5 271 

2 391 

5 271 

Vertical Tail 

2 148 

4 735 

2 148 

4 735 

Riselage 

24 636 

54 314 

24 636 

54 314 

Landing Gear 

13 158 

28 965 

13 138 

28 965 

Nacelle 

8 625 

19 015 

7 410 

16 336 

Structure Total 

(88 743 

(195 647 

(87 528 

(192 968 

Engines 

27 139 

59 832 

^ 24 494 

^ 54 000 

Thrust Reverse rs 

4 809 

10 601 



Miscellaneous Systems 

807 

1 780 

807 

| 1 780 

Fuel System-Tanks and P limbing 

2 622 

5 781 

2 622 

S 781 

Propulsion Total 

(3E 377 

( 77 994 

(27 923 

(61 561 

Surface Controls 

4 527 

9 981 

4 527 

9 981 

Instilments 

] 542 

3 400 

1 542 

3 400 

Hydraulics 

2 540 

5 600 

2 540 

5 600 

Electrical 

2 291 

5 050 

2 291 

5 050 

Avionics 

1 ZlM 

2 690 

1 220 

2 690 

Furnishings and Equipment 

11 590 

25 ill 

11 390 

25 111 

Air Conditioning 

3 720 

S 200 

3 720 

8 20 G 

Anti-icing 

95 

210 

95 

210 

Systems and Equipment Total 

(27 325 

( 60 242 

(27 325 

(60 242 

Weight Empty 

151 445 

333 883 

142 776 

314 771 

Crew and Baggage-Flight , 

306 

67S 

306 

675 

-Cabin, 

744 

1 640 

744 

1640 

Unusable Riel 

1 059 

2 335 

1 059 

2 335 

Engine Oil 

361 

795 

361 

795 

Passenger Service 

4 015 

8 852 

4 015 

8 852 

Cargo Containers 

1 343 

2 960 

1 343 

2 960 

Operating Weight 

1S9 273 

351 140 

150 604 

332 028 

Passengers, (292) 

21 854 

48 180 

21 854 

48 180 

Passenger Baggage 

5 828 

12 848 

5 828 

12 848 

Zero Fuel Weight 

186 955 

412 168 

178 286 

393 056 

Mission Fuel 

1 S 8 680 

349 832 

158 681 

349 834 

Design Gross Weight 

345 635 

762 000 

336 973 

742 890 
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TABLE 11. - BASEPOINT ENGINE WEIGHT 


Item 

Weij^it/Vehicle 

kg 

lb 

Engines (including nozzle 8 thrust reverser) (4) 

24 312 

53 600 

Residual Fluids 

91 

200 

Miscellaneous Provisions 

91 

200 

Engines as Installed 

24 494 



Aerodynamics .- Friction drag estimates were made for a fully turbulent, hydrai- 
lically smooth condition using the inconpressible Von - Ka rman - Schoenherr method 
(reference 7) in conjunction with the adiabatic coirressibility correction of 
Sommer and Short (reference 8). Conponent characteristic lengths (e.g. , the 
distance from the inlet lip to the exhaust nozzle exit, the exposed mean aero- 
dynamic chord of planar surfaces, etc.) and the altitude along the mission 
climb profile were used to evaluate length Reynolds numbers. Flat plate values 
were increased by 3-percent to account for form losses. 


TABLE 12. - BASEPOINT NACELLE WEIGH 


Item 

Weight/Vehicle 


kg 

lb 

Nacelles 

Engine Cowl 

2782 

6132 

Inlet Cowl 

1299 

2864 

Spike 

2961 

6528 

Engine Mounts 

368 

812 

Total Nacelle 

7410 

16 336 
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The wave drag due to thickness was estimated as a function of mach maker 
using supersonic area rule theory (references 9 and 10) in conjunction with a 
transparent wing simulation, an inlet mass flow ratio of one, and *he nozzle 
exit area held fixed at its supersonic cruise position. The effect of inlet 
spillage and nozzle position is included in the installed thrust. All results 
reported here are based on the use of a 51-aach-plane (AX « 0.02 L (#) ) , 

13- roll-angle (A# = 15°) analysis. Basepoint configuration results for 
increased solution mesh density did not indicate any appreciable change. 

Supersonic cruise t rimed drag- due- to- lift characteristics are assmed to 
be equal to the reference configuration of reference 3 and consequently inde- 
pendent of wing and engine size and nacelle shape. A different desijyi wing 
twist and canker is required for each case to realize this performance. The 
lifting efficiency may be conservative for some of the more favorably shaped 
nacelles of the parametric drag study in that any increased benefit that may 
be realized from favorable nacelle thickness/wing lift interference over and 
above that of the reference configuration is neglected. Conversely , for the 
less favorably-shaped nacelles, the analysis may be somewhat optimistic. At 
off design conditions , the above assumption is necessary because the required 
analysis is beyond the scope of the contract effort. 

A conparison of the VSCE 502B (408 kg/sec, 900 lb/sec airflow) nacelle of 
figure 13 to that of the reference configuration nonafterbuming single spool 
turbojet with variable geometry' turbine (363 kg/sec, 800 lb/sec airflow) of 
reference 3 is presented cm figure 14. The basepoint nacelle is 1.95 meters 
(6.4 feet) shorter and has a 0.14 meter (0.46 feet) smaller maxinun diameter. 
The relative cross-sectional shape of the two nacelles is presented in fig- 
ure 15. The basepoint total configuration normal cross-sectional area distrib- 
ution is shown in figure 16. 

Estimated total and nacelle incremental skin friction and wave drag 
characteristics (relative to nacelles off) for the basepoint configuration are 
presented in table 13. The wave drag results are for the case in which the 
nozzle e;:it planes are the sane as the reference configuration. A slightly 
higher drag results (AC^ = 0.00006 at mach 2.7) if the inlet planes are 
matched. 

The friction, wave, and total drag increments of the basepoint nacelle are 
conpared to those of the reference nacelle in figure 17. The basepoint con- 
figuration has a slightly smaller installation drag in spite of 12.5-percent 
greater airflow because of the more favorable nacelle shape (no boattail) as 
shown in figure 14. It was subsequently determined that a further reduction 
of 0.5 count could be realized by meridial contour optimization. 
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Figure 17. -Basenoi nt incremental nacelle drag versus Mach number 





TABLE 13. - RASEFOINT CONFIGURATION ESTIMATED PROFILE AND WAVE 

DRAG CHARACTERISTICS S nt?r . = 929 sq m (10 000 sq ft) 

Ktr 



Altitude 

Aircraft 

Nacelle 

M 

0 

m 

ft 

% 

% 

% 


0.4 

457 

1 500 

0.0061 

— 

0.00065 

— 

0.8 

6 400 

21 000 

0.00572 

— 

0.00062 

— 

1.2 

10 455 

34 300 

0.00545 

0.00365 

0.00060 

-0.00017 

i.4 

11 521 

37 800 

0.00522 

0.00316 

0.00058 

-0.00018 

1.8 

13 594 

44 600 

0.00490 

0.00254 

0.00055 

-0.00019 

2.32 

16 764 

55 000 

0.00450 

0.00222 

0.00050 

-0.00018 

2.7 

18 288 

60 000 

0.00418 

0.00217 

0.00046 

-0.00014 


The aerodynamic charact istics used in resizing the basepoint. wing and 
engine size to produce the baseline configuration used for all parametric 
nacelle drag studies were established as follows. 

Fully turbulent friction levels were adjusted for difference in surface 
area and length Reynolds number of the wing and nacelle. The wave drag varia- 
tion of the basepoint configuration as a function of wing and engine size were 
parametrically evaluated for input to the sizing program. The results are pre- 
sented in figure 18. The effect of engine size was essentially nil at this 
scale for the nacelle shape under consideration. 

The trimmed drag due to lift characteristics were assumed to be independent 
of wing size and equal to the reference configuration. The specific levels used 
are presented in figures 19 through 21 and were taken directly from reference 3. 

Sizing of the basepoint configuration produced the study baseline (table 5) 
which had a 12-percer.t smaller wing size and a 3.5-percent smaller engine size. 
The associated normal cross- sect ic.ial area distribution is presented in fig- 
ure 22. A summary of the component surface areas and reference lengths is 
presented in table 14, and table 15 p 'esents baseline drags. 
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TABLE 14. - BASELINE CONFIGURATION SURFACE AREA 
AND LENGTH SUMMARY 



Swet 

sq.m, (sq.ft.) 

Length 
m. (ft.) 

Fuselage 

786 (8 450) 

96 (315) 

Wing 

1 505 (16 987) 

7.65-39.4 (25.1-129.) 

Nacelles (4) 

276. (3 088) 

12.7 (35.1) 

Center Line Vertical 

20.1 (219) 

4.9 (16.2) 

Wing Verticals 

91. (992) 

7.9 (25.9) 

| Horizontal 

89.5 (921) 

5.8 (18.9) 


TABLE 15. - BASELINE CONFIGURATION ESTIMATED SKIN 
FRICTION AND WAVE DRAG CHARCTERISTICS 


= 929 sq.m. (10 000 sq. ft.) 


M o 

AL 

m. 

TITUDE 

(ft.) 

AIRt 

% 

31AFT 

I % 

NACE 

ac d 

F 

LLE 

AC tw 

0.4 

457 

(1 500) 

0.00568 


0.00065 



0.8 

6 400 

(21 000) 

0.00537 


0.00061 

— 

1.2 

10 455 

(34 300) 

0.00508 

0.00339 

0.00058 

-0.00009 

1.4 

11 521 

(37 800) 

0.00489 

0.00305 

0.00056 

I 

1.8 

13 594 

(44 600) 

0.00455 

0.00237 

0.00052 

-0.00011 

2.32 

16 764 


0.00420 

0.00209 

0.00049 

-0.00012 

2.7 

18 288 


0.00392 

0.00200 

0.00045 

-0.00011 


. rjj 
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r- • r 


•■[,an t e: not filmed 
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Estimation of Supersonic Drag 


Parametric drag analysis .- The parametric nacelle wave drag analysis utilized 
the baseline configuration described in the previous section. The installation 
of the propulsion system followed several general ground rules in order to pre- 
serve the basic arrangement concepts and provide consistent conparisons concern 
ing the effect of nacelle size variations. They are: 

(1) Nacelle overhang of the wing trailing edge and vertical nacelle-wing 
separation was limited to the reference configuration values for structural 
reasons . 

(2) The longitudinal and lateral separation distance between the inboard 
and outboard nacelles was preserved in order to maintain inlet flow quality. 

(3) The reference configuration philosophy of locating the nacelle volume 
in a region of decreasing wing thickness was maintained. 

(4) The maximum boat tail angle considered was 10 degrees. 

The outboard nacelle is moved inboard and forward as required along the 
midchord (approximate maxinun thickness) line of the wing until its trailing 
edge overhang does not exceed 3 meters (Id feet). The inboard nacelle is 
shifted laterally by the same amount holding the longitudinal distance between 
the inboard and outboard nacelle inlet plr es the same as the reference 
configuration. 

The nacelle parametric variables considered in the present analysis were 
the ratio of nozzle area to capture area An/-\;» the ratio of maximum cross- 
sectional area to capture area An^/Ac, the relative axial position of maximum 
area X^YX'/f » the ratio of nacelle length to capture diameter, l/dc, and the 
nacelle absolute capture area A summary of the number of variations and 
variable range analyzed is presented in table 16. For purposes of computation, 
the nacelles were assumed to be axisymmetric and the inlet, maximum area, and 
nozzle planes to be connected by straight lines. 
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TABLE 16. - NACELLE PARAMETER VALUES 


PARAMETER 

VALUES 

Mach nuaber 

1.2, 2.32 

V*c 

1.0, 1.2S, l.S, 2.0 

A MAX /A c 

1.0, 1.2S, 1.5, 2.0 

X AMAX /1 

0.4, 0.6, 0.8, <1.0 * 

A c 

1.86, 2.79, 3.72 sq.m. 
(20, 30, 40 sq.ft.) 

l/ d c 

S.S and 7.0 

* Maxinun value considered corresponds to a boat tail angle 

1 of ten degrees. 
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The parametric nacelle friction drag analysis is based cm the use of fully 
turbulent flat plate levels in conjunction with the expression for surface 
areas (for four nacelles) : 



( 



» 2 


S , d 
ref c 



The largest deviation between the exact and approximate express’ on occurs for 
*AMAX/f approaching 0.4 and An/A^. approaching 2.0 with the former resulting in 
10-percent greater area. It will be subsequently found that these differences 
are negligible in terms of the total installation drag for such cases. The 
parametric nacelle friction results are presented in reference 2. 

Nacelle normalized cross- sectional area parametric extremes of the present 
study are presented in figure 23. Jtaximum- to -capture area ratio of 1 to 2 at 
40, 60, and 80 percent of the nacelle length are shown for nozzle-to-capture 
area rations of 1 and 2. Wave drag results are discussed in detail in 
reference 2. 

Briefly, the incremental nacelle wave drag is a strong function of the 
ratio of maximum- to-capture cross-sectional area, Amav/Ac, hoattail area, and 
to a somewhat lesser extent relative axial position of maximum cross-sectional 
area, Nacelle shapes with negative wave drag exist because of favor- 

able total system thickness interferences associated with the location of grow- 
ing nacelle cross-sectional area in a region of decreasing wing thickness. The 
nacelle geometric variable behavior and sensitivity are unchanged by mach num- 
ber, nacelle capture area, or nacelle fineness ratio. The incremental wave 
drag results are, in general, weak functions of the latter two variables for 
efficient installations. 

Detailed nacelle wave drag variations with freestream mach number were 
defined for a range of levels covering high-positive, zero, and negative 
installation increments. These characteristics correspond to nacelles with 
large maximum cross- sect ional area relative to ‘he capture and nozzle area 
cylindrical, and near- truncated conical shapes, respectively. Figure 24 illus- 
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trates the mach number difference for these extremes for the fineness ratio 5.5, 
medium-size nacelle. Examination of the results indicate that weak to moderate 
mach number variations are associated with small nacelle installation drags. 
Conversely, strong compressibility variations are exhibited for inefficient 
installations. The large benefit at transonic speeds is somewhat illusory as 
the thrust must be progressively penalized for nozzle contraction with 
decreasing mach numbers. 

Drag table look-up computer program .- A table look-up computer program was 
developed ( appendix ) which yields the incremental wave and friction drags of 
nacelles as functions of nacelle geometry variables and airplane mach number. 

The drag increments are for the total vehicle relative to the vehicle with 
nacelles removed. The nacelle shape parameters used as inputs to the program 
are: 

Inlet capture area 

Nacelle maximum cross-sectional area 
Nozzle exit area (supersonic cruise position) 

Distance from inlet cowl leading edge to maxinun cross-sectional 

Nacelle total length 
Reference wing area 
Ihe output of this program includes for the nacelle of interest: 

(1) The aforementioned input data 

(2) Drag coefficients at mach 1.2, mach 2.32, and the input mach number 
for friction (CTF) , wave (CDW) , and total (CDO) drags 

(3) The nondimensional parameters of position of maximum cross-sectional 
area ( x AMAx/^)> nozzle- to-capture area ration (An/ Ac) , maximum- to-capture area 
ration (A|^/ A c )> an <l fineness ration ( £ /d c ) 

In addition, incremental drag coefficients of the reference airplane nacelle 
(reference 3) are printed. A sample output is shown in table 2. It has been 
found that the table look-up results correlate best with more detailed analyses 
when the maximum cross-sectional area and its position are based on the area 
that occurs at the intersection of straight lines originating from the inlet 


(1) 

A c 

(2) 

i W 

(3) 

A 

n 

(4) 

*MAX 

(5) 

l 

(6) 

S REF 
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which is tangent to the maximum slope of the forebody and from the nozzle which 
goes through the actual maximum area and whose slopes nearly match the slopes 
of the actual nacelle, as illustrated in figure 25. This method most closely 
approximates the parametric nacelle shape drag analysis because, in that study, 
nacelle shapes were defined by two straight-line segments. Using this 
method will result in nacelle drag increments at supersonic cruise within 0.5 
drag count of the drag resulting from a detailed analysis. 

Nacelle shape estimation. - Nacelle external shapes are determined by such 
installation items as engine accessories, compartment cooling, shrouds and 
insulation, aircraft accessories, engine clearance, engine mount geometry, 
nacelle structure, boundary layer gutters, etc. The engine configurations sup- 
plied by the engine manufacturers usually include only the engine case outline 
and nozzle dimensions. A method has been established to determine the engine 
external envelope and aircraft structure and equipment space allowances. 
Guidelines are presented for determining the engine buildup envelope, engine 
cowl, nozzle fairing, and inlet and inlet cowl shapes. Those installation 
iters which have the largest effect on nacelle shape are then discussed. An 
exam} ? of the nacelle shape buildup is presented in figure 1. 

Engine buildup envelope: The procedure to establish the engine buildup 

envelope is: 

(1) Establish fan and gas generator case outline. 

(2) Add 5 cm (2 in.) constant to all surfaces of the preceding outline to 
provide for wiring, plumbing, etc. 

(5) Add 2.5 cm (1 in.) constant additional to outline for variable com- 
pressor geometry mechanisms where applicable. 

(4) Establish mechanical power extraction drive station and radial loca- 
tion for engine accessories drive and for aircraft accessories power takeoff. 
Depending on engine configuration and accessory design, accessories may be on 
the engine or in the pylon or wing. For engine accessories where encapsulation 
is required for cooling, provide 0.595 cu m (21 cu ft) of volume proximate to 
engine accessories drive of item (4). Dimensions of the capsule may be varied 
for best packaging, but the capsule thickness at the gearbox should be 0.305 m 
(12 in.) minimum. For nonencapsulated engine accessories, provide 0.51 cu m 
(18 cu ft) of volume proximate to engine accessories drive of item (4). Arrange- 
ment of the accessories package may vary, but the minimum thickness at the 
gearbox must be 25 cm (10 in.). The dimensions of items (2) and (3) and (4) 
are additive and will usually establish the maximum radial dimensions of 
the gas generator section of the engine. All other engine and aircraft equip- 
ment in this portion of the nacelle should be contained within the volume of 
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revolution established by this line of radii. For minimum aircraft drag, 
locate aircraft accessories within pylon or adjacent wing. The engine oower 
takeoff pad and angle drive gearbox will be within the pylon and this will not 
impact the nacelle mold line. 

(5) For engine fluid reservoirs, add 0.0566 cu m (2 cu ft'' on left or 
right side of engine. The radial dimension is 15.2 cm (6 in.) additive to 
dimension of item (2) or (3) . 

(6) Equally space four compressor high-pressure bleed ports around com- 
pressor rear frame. The 12.7 cm (5 in.) diameter collector manifold (with 
flexible sections) will interconnect the ports and connect to the airframe duct 
in the pylon. The collector manifold diameter is additive to the allowat 

item (2). Low-pressure bleed ports may be provided in place of or in adci ’ 
to the preceding. The same space allowance must be made for these. Or. or 
more engine and inlet anti-icing air ducts will be routed from the blee<. uu 
fold forward to the engine front frame. These ducts will be 10.2 cm (4 m.; 
diameter and will be additive to the allowance of items (2) or (3). 

(7) Determine location of engine mounting stations by engine internal 
structure. The nacelle/wing relationship should be considered in the placement 
of the engine load- carrying frames and the mounting provisions on the frames 
for best weight effectiveness of the total system. 

(8) Main mounts - transfer thrust, side, and vertical loads: Provide 

15. 2 x 15.2 x 12.7 cm (6 x 6 x 5 in.) radial space additive to items (2) or (3) 
at two positions on engine main mount frame, circumferentially spaced greater 
than or equal to 90 degrees. 

(9) Stabilizer mount - transfeis vertical loads: 1 "ovide 10.2 x 10.2 

x 20.3 cm (4 \ 4 x 8 in.) radial space, additive to items (2) and (3). 

(10) Locate local protrusions of miscellaneous engine equipment beyond 
envelope of items 12) and (3) to occur at random locations. These will be 
relatively small and will not exceed the maximum envelope noted in item (4). 

Engine Cowl: Hie engine cowl shape may be determined by the following: 

(1) The cowl inner skin mold line must maintain a minimum 2.5 cm (1 in.) 
clearance from all points on the engine buildup envelope developed in the 
preceding discussion. 
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(2) Cowl structural requirements will vary depending on the arrangement 
and location of engine mount points. Where mount points are in the proximity 
of the pylon structure, the cowl can be made nonstructural ; i.e., sufficient 
to withstand internal and external airloads and flight dynamic forces. This 
will require 5 cm (2 in.) of structure (constant) from the engine front face 
station to the cowl- to- nozzle fairing interface. Where mount points are widely 
separated from the pylon structure a structural cowl must be provided to trans- 
fer the engine loads. Cowl thickness in the load paths will be 7 to 10 cm (3 to 
in.). In these areas the structure can intrude into unoccupied space in the 
volume of revolution developed in the preceding items (4) and (5), but 2.5 cm 
(1 in.) clearance must be maintained from adjacent engine buildup equipment. 

Areas of the cowl outside the load paths can be 5 cm (2 in.) thickness. 

Nozzle fairing: The nozzle fairing . ust fair smoothly into the engine 

cowl mold line developed in the preceding and fair smoothly to the base diameter 
dictated by the nozzle. The nozzle fairing leading edge step height from the 
engine case will vary depending upon the engine services (hydraulics, pneumatics, 
fuel, secondary' airflow, etc) required to pass through it. This step height may 
require adjustments to the cowl outer mole line as it approaches the nozzle 
fairing interface. 

Inlet and inlet cowl: Inlet and cowl shape can be determined by the fol- 

lowing method: 

(1) Establish inlet length and capture area based on appropriate nacelle 
design methodology and external constraints. 

(2) Establish inlet flow path area geometry. 

(3) Define inlet cowl external lines. Fair from inlet lip to engine 
cowl. The faired mold line should provide minimun rate of cross-sectional area 
increase. 

(4) Establish requirement for the following airflow t afhs appropriate to 
inlet geometry' and engine cycle used: 

(a) Auxiliary' air inlet 

(b) Bypass air 

(c) Engine secondary air 
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(d) Boundary layer bleed 

(i) Centerbody or ranps 
(ii) Cowl inner wall 

(5) Determine cowl wall thickness by requirements of structural integrity 
plus space required for flow paths and door mechanisms associated with require- 
ments established in item (4). A minirun thickness of 15.2 cm (6 in.) is 
suggested for the cowl wall from the inlet throat aft to the engine front face. 
This .nay be varied locally, but internal lines should be maintained. Thickness 
of the cowl structure will vary from approximately 0.16 cm (1/16 in.) at the 
inlet lip to the throat thickness established in the preceding. 

Major nacelle shape elements: Three major elements establish the engine 

external envelope: 

(1) The allowance over the total surface of the gas generator of space 
for engine variable geometry- mechanism, plumbing, wiring, etc. 

(2; The space required for the engine accessory- gearbox and associated 
accessories. This package establishes the location and magnitude of the nacelle 
maximun cross-sectional area. 

« 5) The location of the main engine mounts on the engine as defined by 
the engine manufacturer. Khcre the main mounts are placed at the compressor 
front or midframe, sufficient structure is available in the adjacent nacelle, 
pylon, and wing to carry the multidirectional loads, and a simple, nonstructural 
cowl may be used. Where the main mounts are placed at the turbine frame, it is 
necessary to consider the cowl as a structural cylinder with penalties to the 
nacelle sire and weight. 

All other elements of the engine installation fall within the envelope 
defined by tne preceding. 


list imat ion of Naccl lc/lnlet Weight 

Weight estimation of nacelle and inlet systems is a complex process and 
requires design detai’ not normally performed in the type of preliminary- 
studies being considered here. The estimating procedure described in this sec- 
tion uses a simplified approach producing a first-order-type weight estimate 
keyed to gross definitions of the nacelle/ inlet package. The procedure defined 
will provide the capability of maintaining consistency between nacelle weight 
estimates while making comparative analyses. 
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Two dimensional inlets.- To estimate the weights of engine nacelles with 
two-dimensional (2-D) inlets, the nacelle package is divided into the following 
canponents: 

(1) Engine cowl 

(2) Inlet cowl 

(3) Ramps 

(4) Air induction special features 

(a) Bypass system 

(b) Auxiliary inlet 

(c) Secondaiy air provisions 

(d) Inlet controls 

(5) Engine mounts 

The methods used to estimate the weights of the nacelle components are 
primarily based on a prior Rockwell inlet study for the Boeing SST. This 
study was conducted for Boeing and consisted of designing a 2-D inlet as a 
contender to be compared to Boeing’s ax i symmetric inlet design in the inlet 
selection for the SST. Unit weights used to estinate wights of the nacelle 
components were derived from data developed for this study. 

Engine cowl: The engine cowl is defined as the total nacelle structure 

aft of the engine front face, including all structure that supports and sur- 
rounds the ermine. The engine cowl weight is estimated at 34.2 kg/sq m 
(7.0 lb/sq ft) of nacelle external wetted area. 

Inlet cowl: The inlet cowl is defined as the total nacelle/inlet structure 
forward of the engine front face, exclusive of the variable-geometry ramps and 
special air induction features. The inlet cowl weight is estimated at 24.4 kg/ 
sq m (5.0 lb/sq ft) of wetted area. 

Ramps: The ranps are defined to be the movable panels, including an actua- 

tion system, used to vary the inlet geometry in a 2-D variable-geometry inlet. 
Weight of the variable-geometry ranps is estimated at 48.8 kg/sq m (10.0 lb/ 
sq ft) of movable ramp planform area. 
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Air induction special features: The bypass system consists of inlet air 

bypass doors, including actuation provisions. The system weight is estimated 
at 39.1 kg/sq m (8.0 lb/sq ft) of door area. 

The auxiliary' inlet is defined as the auxiliary air inlet doors and inlet 
actuation system. The weight of this system is estimated at 29.3 kg/sq m 
(6.0 lb/sq ft) of door area. 

The secondary’ air provisions provide inlet air to the engine compartment 
for engine compartment cooling. The weights of these provisions are estimated 
with the following equation: 

/» \ 1/2 /* \ 

” = 22 ’(li?) or hT ' 50 \lk) lb 

where K is engine design airflow, kg/sec (Ib/sec). 

The inlet controls are defined as the system provided to monitor the inlet 
conditions and transmit position signals to the movable inlet systems. The 
weight of this system is estimated at 22.7 kg (50.0 lb) per inlet. 

Hngine mounts: The engine mounts are the fittings used to support the 

engine in the nacelle. The weights of these fittings are estimated at 1.5 per- 
cent of the engine weight. 

.Vxi symmetric inlets .- To estimate the weights of engine nacelles with 
ax i symmetric inlets, the nacelle package is divided into the following compo- 
nents; this breakdown is similar to the one described for a 2-D inlet: 

I ) Ungine cowl 

(2) Inlet cowl 

(3) Spike 

(4) Hngine mounts 
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The procedures for weight estimation of the engine cowl and mounts are 
the same as those described for the 2-D inlet/nacelle. The methodologies to 
estimate the weights of the inlet cowl and spike were obtained from the Air 
Force Technical Report SEG-TR-67-1, Preliminary Design Methodology for Air 
Induction Systems. 

Inlet Cowl: The inlet cowl is defined as the total nacelle/ inlet structure 
forward of the engine front face, exclusive of the inlet spike and its systems. 
The inlet cowl weight is determined by the following statistical equation. 

wr = 0.159 (N) [(A c ) 0 ' 5 L (P 2 )] 0 ' 731 kg or 
7.435 (N) [(A )°' 5 L (P^] 0 ’ 731 lb 

where: 

N = number of inlets 

A c = capture a.ea per inlet - sq m (sq ft) 

L = subsonic duct length per inlet - m (ft) 

P 7 = maximum steady- state static pressure at engine face at 
supersonic cruise mach - kg/sq m (psia) 

Spike: The spike is defined to be the center body structure, including 

its systems and actuation. Weight of the spike is estimated with the following 
statistical equation; 


hT = K (N) A 

c 

where: 

K = 252.9 kg/sq m (51.8 lb/sq ft) 

N = number of inlets 

A c = capture area per inlet = sq m (sq ft) 
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Estimation of Airplane Takeoff Gross Weight 

Weight sensitivity analysis. - In considering possible trades of reduced 
drag through design changes in the engine envelope for some penalty in engine 
weight and performance, it is necessary to have visibility of the net impact of 
all these effects on the total airplane system. To evaluate these effects, the 
sensitivities of the airplane takeoff gross weight to variations of propulsion 
system parameters were determined. These sensitivity data were obtained by 
conducting design trades on the baseline airplane for variations of the follow- 
ing items: 

(1) Incremental nacelle drag 

(2) Propulsion system weight 

(3) Engine specific fuel consunption 

(4) Engine sizing condition thrust 

In each case, the parameter of interest was varied independently and the 
airplane resized to the design mission range of 7408 km (4000 n mi) while 
maintaining thrust -to- weight and wingloading values equal to those for the base 
line vehicle. 

Incremental nacelle drag: Several variations of nacelle drag were invest! 

gated. These were chosen as representative of the combined wave and friction 
drag variations as found in the nacelle shape analysis to allow use of the 
trade data for any nacelle geometry analyzed in this program. 

The results of this trade are shown in figure 26, which shows relative 
takeoff gross weight (TOGW) versus nacelle drag at roach 2.32 for several varia- 
tions of the drag increment at mach 1.2. 

Propulsion system weight trades: Airplane TOGW was calculated for several 
propulsion system weight increments. Incremental propulsion weight, in this 
case, is defined as a percent of the sun of the engine, nacelle, and miscellane 
ous propulsion systems (198 kg, (44S lb) per nacelle) weights. The results 
of this trade are shown in figure 27, which plots relative TOGW versus 
propulsion weight increment. 

Engine specific fuel consunption trades: Four separate trades were 

performed with SFC increments applied independently to the following mission 
segments : 

(1) Maximum power climb legs only 

(2) Supersonic cruise legs only 

(3) Subsonic cruise and loiter legs only 

(4) The entire mission 

The results of this trade are presented in figure 27 as relative TOGW versus 
percent change in SFC. 
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Relative TOGW Relative TOGW 
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Engine sizing thrust: Figure 28 presents relative TOGW versus percent 

change in thrust at the engine sizing condition. In this trade, it is assuned 
that the thrust available at the engine sizing condition varies without 
change in propulsion characteristics at other flight conditions. In ihe 
current study, all engines were sized at takeoff. 

Application of sensitivities . - To demonstrate the method of applying the 
vehicle sensitivities, an example using the VCE 112C engine is in the follow- 
ing paragraphs. The VCE 112C is discussed in detail later wider "Validation 
of the Approximate Method." Additional exanples are also given therein. 

The baseline airplane characteristics and sensitivities were originally 
computed using a friction drag that was approximately one count too high 
and an ambient tenperature increment that was incorrect. Thus the baseline 
airplane should have been somewhat lighter. The baseline airplane was 
recouputed and resulted in a takeoff gross weight of 316 783 kg (698 375 lb), 
a propulsion system weight (four nacelles) of 30 882 kg (68 081 lb) , and drag 
coefficients of 0.0040 and 0.0030 at mach 1.2 and 2.32, respectively. Figure 
26 lias been revised relative to that shown in reference 2 for this reason. 

A nacelle drawing (figure 29) was made with engine accessories located cm 
the engine without encapsulation for cooling. Weight for the nacelle was 
estimated to be 8750 kg (19 298 lb), which includes 198 kg (445 lb) for 
miscellaneous propulsion systems. The revised baseline nacelle weighed 
7880 kg (17 020 lb). Thus, the VCE 112C nacelle is 12 percent heavier than 
the baseline. From figure 27, the relative TOGW ratio, R,^., for this change 
i'. 1.059. 

In order to maintain the takeoff distance, a new engine must have the 
same effective thrust -to- weight ratio between 0.0 and 0.3 mach (approximate 
liftoff speed) as the baseline. The effective thrust occurs at approximately 
mach 0.25. Because the VCE 112C has a significantly different thrust lapse 
rate with mach number than the baseline (as shown in figure 30) the effective 
thrust is 17 percent lower than the baseline. Extrapolating figure 28 to a 
thrust increment of 17 percent yields a relative TOGW ratio due to takeoff 
thrust, of 1.10. 

Figures 31 and 32 show installed performance of the 100-percent size VCE 
112C and VSCE 502B propulsion systems. Because the Breguet range factor, 

M x L/ D/ SFC, maximizes near minimum SFC, the airplane will tend to fly 
at or near minimum SFC; cruise altitude will be adjusted to achieve this. 
Therefore, the SFC increment may be taken at the minimum of each engine. 

This assumption is slightly optimistic because the lift-drag ratio will also 
change and will affect the operating point. Thus, the VCE has about 6.3- 
percent higher SFC at supersonic cruise and 4.6 percent lower SFC at subsonic 
cruise than the baseline. From figure 27, this results in a relative TOGW 
due to a change in supersonic cruise SFC, R^, of 1.071 and due to a change 
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Sizing point thrust increment ~ percent 


Figure 28. - Sizing point thrust sensitivity trade 
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Figure ?9.- VCE 112C nacelle. 
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Figure 30.- VSCE 502B and VCE 112C takeoff thrust. 
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Figure 31 - Installed performance comparison of baseline 
and MCE 112C at Mach 2.32, 19 800m (65 000 ft). 
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Figure 32. * Baseline and VCE 112C performance, Mach 0.9, 13,700m (45 000 ft). 




in subsonic cruise SFC of 0.994. 

The nacelle shape of the VCE 112C is similar to that of the baseline 
but has a smaller nozzle area. The drags used in the mission analysis 
computer program are indicated in table 17. From figure 26, the relative 
TOGW ratio due to change in drag, R^, is 1.021. 


TABLE 17. - COMPARISON OF BASELINE AND VCE 112C NACELLE DRAG INCREMENTS 


Mission 

analysis conputei program values 


d< T)’ 

fC D’ 


mach 1.2 

mach 2.32 

VSCE 502B (baseline) 

0.00040 

0.00030 

VSCE 112C 

0.00085 

0.00048 


The ratios R^, Rp^» ^FC’ are t * ien “^iLiplied together to 


obtain the total relative TOG!*.' ratio, 

SlTAL = \l X R FNE X X ^^SUPER 


, ^ 


= 1.266 


The new TOG'. is 


obtained by multiplying Rj-gr^ by the baseline 


316 783 kg (698 375 lb). Thus, it is estimated that a vehicle 
performance requirements using the VCE 112C engine would weigh 
(884 258 lb). 


TOGV of 

meeting the 
401 092 kg 


Validation of the Appioximate Method 

Standard preliminary resign procedures were applied in the installation 
of four candidate engines in the baseline supersonic transport airplane. 

Drag and weight estimates were made utilizing conventional procedures. The 
airplanes were then sized to the design mission utilizing an autt mated 
reiterative process. The results of this task provide a more exact 
evaluation of the candidate engines thar. is obtainable with the more approxi- 
mate methods resulting from the earlier task, and thus serve as a reference 
for evaluation of the applicability of the approximate methods. 
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Propulsion. - Four engines were selected so that the sensitivity Method could 
be validated for a wide range of engine types: 

(1) PWA VSGE 502B duct-burning turbofan 

(2) PWA VCE 112C variable-cycle engine 

(3) GE GE21/J10 B1 low bypass turbojet 

(4) GE GE21/J11 B3 double-bypass variable-cycle engine 

Characteristics of these engines are summarized in table 18. The engine 
and installations are discussed in the following paragraphs. While there 
may be different weight and performance margins and noise and technology 
as simp t ions for the four engines, the data were used as supplied by the 
engine manufacturers without modification for these differences. Hence, 
results discussed later should not be interpreted as representative of 
the final SCAR engine studies. (Engine studies are currently still 
under way.) 


TABLE 18.- ENGINE SUMMARY 


ENGINE 



vsce so:b 

VCE 112C 

GE21/J10 B1 

GE21/J11 B5 

I 

Design airflow, 

4ns 

498 

518/541 

518/582 

kg/ sec (lb/ sec) 

(900) 

(900) 

(700/750) 

(700/840) 

Bypass ratio 

,3 

2.5 

0.1 

0.5 

Sea- level static 

26 500 

25 000 

25 500 

18 500 

takcof f thrust , 

(59 500) 

(5b 200) 

(52 300) 

(41 100) 

daX (lb)* 





i Takeo f f spec i f i c 

74 

^1 

75 

55 

i thrust, daN/kg/ 

(~b) 

(72) 

(77) 

(56) 

sec (lb/lb/scc)* 





jPry weight. 

6077 

0191 

7280 

5964 

; kg (lb) 

(13 400) 

(13 b50) 

(16 050) 

(15 150! 

i Overall length. 

0.8 

7.9 

7.0 

6.9 

! m (in.) 

(2bb) 

(510) 

(275) 

(275) 

Maximun diameter. 

2.24 

2.19 

1.96 

2.01 

m (in.) 

(88) 

(86) 

(77) 

(79) 

*As installed for 

this study, standard plus 8 

3 C day 
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VSCE S02B: Refined performance data (dated January 1976) for the VSCE 502B 

engine were used to calculate revised installed propulsion performance data. 

All performance installation effects used were the same as for the baseline 
data. The noise calculation procedure was identical to that for the baseline. 
Installed performance changes were very small relative to the baseline. Per- 
formance data at important flight conditions are shown in figures 30, 33, and 34 

The nacelle was revised slightly compared to the baseline, to include 
engine accessories (unencapsulated) and more realistic structure allowances and 
boundary layer diverter, as shown in figure 35- This engine configuration lent 
itself well to establishing an efficient nacelle shape. A well-defined waist at 
the compressor m.j frame provided space for the required accessory geartox vol- 
une without forcing the nacelle maximim diameter much beyond the nozzle diameter 
Location of the engine main mount at the engine . r ont frame enabled transfer of 
the mount loads directly into the pylon/wing structure. 

VCE 112C: Installed performance data for the VCE 112C (408 kg/sec (900 lb/ 

sec)) were calculated in the same manner as the baseline. Because the super- 
sonic cruise airflow is the same as that of the VSCE 592B, the same capture area 
was used. The only significant difference in procedures was that only a 4 dB 
reduction in noise due to the coannular effect was used (instead of 8 dB) 
because the exhaust characteristics and nozzle configuration of the VCE 1’2C are 
such that an 8 dB reduction could not be achieved. This results in some thrust 
reduction at nvich 0.5 takeoff power (figure 30) while the airplane is on the 
ground in order to stay within FAR 56 noise requirements. Thrust is reduced 
even further at the takeoff noise measurement point. Figures 33 and 34 compare 
installed performance of the VSCE 502P and VCE 1'2C. 

The installation of this engine into a nacelle (figure 29) is quite similar 
to tint of the VSCE 502B. The only significant changes arc a longer engine and 
a slightly reduced nozzle diameter. 

CE21/J10 Rl: Installed performance of the (3:21/110 R1 was calculated in a 

manner similar to the baseline except that a 2-0 mixed conprcssion inlet with a 
capture area of 2.0" sq m (5 208 sq in.) was used . The CE2I/J10 BI has appioxi- 
matcly 16-perccnt lower supersonic cruise airflow relative to takeoff airflow 
than does the VSCE 502B. If an axisymmctric inlet had been used and sized for 
;■ personic cruise, the static takeoff inlet recover)' would have been 4 percent 
lower than the VSCE. 502R. Ti.e 2-0 inlet has more throat area variation capabil- 
ity and larger auxiliary doors than the ax i symmetric inlet. Thus, the takeoff 
recovery is actually slightly higher than that for the VSCE a!) 2 B/a x i synnic trie 
inlet. Pei formancc data at important flight conditions are presented in 
figures 36 through 38. The noise calculation procedure was identical to that 
for the baseline. An exhaust noise reduction of 8 dB at all flight conditions 
was used due to the coannular noise reduction effect. 


80 














Fnc net propulsive effort - 1000 daN (1000 lb) 



STANDARD + 8° C DAY 



Figure 56.- GEZ^JIO Bl and CE21/J11 83 takeoff thrust. 
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The near-cylindrical configuration of this engine penalized the nacelle 
shape when the required volune for the encapsulated engine accessories package 
was added (figure 39). The package was shaped to minimize the added cross 
section. Location of the main mount at the turbine frame forced the engine cowl 
to be designed as a structural element and added more to the required cross 
section. 


GE21/.J11 B5: Relative to the GE21/.J10 Bl, the GE21/.J11 B3 has even higher 

takeoff airflow and the same supersonic cruise airflow, thus creating a greater 
takeoff airf low/inlet matching problem. Therefore, the 2-D inlet was used. 
Because the engine operates without augmentation at takeoff and the exhaust velo- 
cities are low, the sideline noise is quite low while the airplane is on the 
ground (approximately b dB below FAR 56 requirements ) . However, this signifi- 
cantly reduces takeoff thrust (figure 36) and increases engine size to meet 
takeoff distance requirements. Some thrust cutback is still required at the 
takeoff noise measurement point. Installed performance is compared to the 
GE21/J10 Bl in figures 37 and 38- 

This engine is similar in configuration to the GE21/J10 Bl except that the 
main engine mount was located at the compressor midframe. This enabled the cowl 
to be made nonstructural and reduced somewhat the nacelle maximum cross section. 

A 1 J other details of the engine and inlet lines development for the two engines 
arc identical, as shown in figure 40. 

Mass properties. - height estimates were made for the four candidate engine 
installations. The engine weights, including nozzles and thrust reversers, were 
furnished by t ho engine manufacturers. Weight increments for residual fluids and 
miscellaneous engine/airframe interfacing provisions were added to the manufac- 
turer's quoted weights to obtain the installed engine weight. The weight sum- 
maries of the engines arc presented in tables 19 and 20. The nacelle/inlet 
weights for these engine installations were calculated from their respective 
nacelle lines development layout drawings figures 29, 35, 39, and 40. Summar- 
ies of the nacelle/ inlet weights arc shown in tables 21 and 22. The differ- 
ence in nacelle weights between the PWA YSCE 502B, and VCE 112C engine instal- 
lations is primarily due to the longer engine cowl length of the VCE 112C instal- 
lation. For a given inlet capture area, an axisymmctric inlet would weigh less 
than a 2-D inlet. However, the 2-D inlet/nacelles for the GE GE21/.J10 Bl and 
CE21/.J11 B5 are of similar weight as the axisymmctric FWA engine nacelles. 

This results from the smaller inlet capture areas and engine dimensions cf the 
GE eng ires. 


Aerodynamics . - Normalized cross-sectional area shapes of the four candidate 
nacelles of this study arc presented in figures 41 and 42. The estimated 
nacelle incremental skin friction, wave drag, and total drag characteristics 
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TABLE 19.- ENGINE WEIGHTS, INTERNATIONAL UNITS 



k i log rams/ Eng i ne 


Pratt § Whitney 

General Electric 

Engines 

VSCE 502B 

VCE 112C 

GE21/J10 B1 

GE21/J11 B3 

Bare engine (including 
nozzle § thrust 
reverser) 

6077 

6191 

7280 

5964 

Residual fluids 

23 

23 

23 

23 

Miscellaneous provisions 

23 

23 

23 

23 

Engine as installed 

6123 

6237 

7326 

6010 


TABLE 20. ENGINE WEIGHTS, ENGLISH UNITS 



Pounds/engine 


Pratt § Whitney 

General Electric 

Engines 

VSCE 50 2B 

VCE 1 J2C 

GE 21/J10 B 1 

GE21/J11 B3 

Bare engine (including 
nozzle § thrust 
reverser) 

13 400 

13 650 

16 050 

13 150 

Residual fluids 

50 

50 

50 

50 

Miscellaneous provisions 

50 

50 

50 

50 

Engine as installed 

13 500 

13 750 

16 150 

13 250 


PRECEDING PAGE BLANK NOT FILM** 
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TABLE 21.- NACELLE/ INLET WEIGHTS, INTERNATIONAL UNITS 


Nacelle/ inlet 

Ki lograms/nacelle 

Axisymmetric 

Two-dimensional 

VSCE 50 2B 

VCE 112C 

GE21/J10 B1 

GE21/J11 B3 

Engine cowl 

816 

1124 

706 

760 

Inlet cowl 

3"0 

565 

566 

598 

Spike 

^40 

735 



Ramps 



51' 

517 

*Air induction features 





Bypass 



49 

49 

Auxiliary inlet 



50 

30 

Secondary air provisions 





Inlet controls 



45 

45 

Engine mounts 

92 

95 

110 

90 

Total nacelle/inlet 

2018 

2515 

2047 

2113 

* Included with spike weight in axisynmetric inlets. 


TABLE 22.- NACELLE/ INLET WEIGHS, ENGLISH UNITS 



Pounds/nacelle 


Axisynmetric 

Two-dimensional 

Naccl le/ inlet 

VSCE. 50 2B 

VCE 112C 

GEZ1/J10 B1 

GE21/J11 B3 

Engine cowl 

1798 

2477 

■ 

1673 

Inlet cowl 

816 

804 


1319 

Spike 

1632 

1616 



Ramps 



1140 

1140 

*Air induction features 





Bypass 



107 

107 

Auxiliary inlet 



67 

67 

Secondary air provisions 



55 

55 

Inlet controls 



100 

100 

Engine mounts 

205 

206 

242 

199 

Total nacelle/ inlet 

1 

4449 

5105 

4515 

4658 

*Included with spike weight in axisymmetric inlet 

s 
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for all of the nacelles are presented in table 23. Inducted in the table 
are the incremental drags obtained from the computer table look-up program 
of the parametric drag analysis results. The differences between the estimated 
drags and the parametric results are attributed to the following considerations: 

(1) The parametric study, which is the basis of data for the table look-up 
program, was made with a resized wing with the nacelles relocated further inboard 
and forward in accordance with the ground rules of the study. 

(2) The current nacelles deviate from a linear radius connection of the 
inlet face, maximum area, and nozzle stations. 

(3) In some cases, the ratio of the distance to the maximum area to the 
total length (X^^y/I) falls outside of the parametric study envelope, thereby 
requiring extrapolation. 

(4) Drag coefficients for intermediate mach numbers are obtained by using 
a cubic curve fit based on the total drag (OK)) increments at mach 1.2 and 
2.32. However, since wave drag does not follow such a sinplified solution, 
interpolated or extrapolated drags will deviate from an estimated value by 
varying amounts. In the case of the reference nacelle, the deviation was 
-0.00017 at mach 1.4. 

Performance and sizing . - Performance calculated for the aircraft having the 
four selected propulsion systems installed includes all items as described 
earlier under "Baseline Airplane Definition." In addition, performance has been 
calculated for an "economic" mission as described in the following paragraphs. 

As on the design mission, the economic mission is calculated for a standard- 
plus-8 0 C (14.4° F) day. 

A profile of the economic mission is shown in figure 43 . This mission 
consists mainly of a mach 2.32 cruise as in the design mission; however, it is 
preceded by a subsonic climb and cruise totaling 741 km (400 n mi). Fuel 
reserves are calculated just as in the design mission for an alternate airport 
located 463 km (250 n mi) beyond the destination airport. The economic mission 
is an off-design mission in that the airplane, as sized to 7408 km (4000 n mi) 
range on the design mission, carries a reduced payload equal to 55 percent of 
the design payload and fuel is then off-loaded to yield a total economic mission 
range of 4630 km (2500 n mi) plus fuel reserves. 

The economic mission consists of the following legs: 

(1) Warmup and takeoff - 10 minutes at idle power plus 1 minute at 
maximum power 
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TABLE 23. - DRAG COMPARISON 


NASA LARC CR132374 Reference Configuration 


Estimated 

Parametric 



M 




AC r 

r, P 

A %. 

A \ 

1.2 

0.00072 

-0.00041 

0.00031 

0 . 00)881 

-0.00034 

0.00031 

1.4 

0.00070 

-0.00017 

0.(88)52 



o: 0003S 

1.8 

li.OOOCiS 

-0.(88)22 

0.00045 



0.0(8)40 

2.52 

O.IHHKjO 

-0.00018 

0.00042 

0.00057 

-0.00014 

0.(88)43 

1 “ 

(l.(8B)5b 

-0.00015 

0.00041 



0.00044 

Base-point VSO: S02B 


estimated 

Parametric 

M 

A \ 

A \ 


A S 

A \ 

AC °„ 

1.2 

o.ooooo 

-0.(888)7 

0.(8)055 

0.00055 

0.(888)2 

0.00057 

1.4 

0.00058 

0.(88)11 

0.00069 




1.8 

0.188154 

-0.(888)4 

0.(88)50 



0.0(8)46 

2.52 

n. 0(8)50 

-0.(888)8 

0.00042 

0.00049 

-0.0(88)7 

0.00042 

2 . 7 

0.(88)47 

-0.(888)9 

0.(88)58 




Baseline VSCli 502B 


list i mated 

Parametric 



M 

A \ 

W 

AC »„ 

A S- 

AC °* 


1.2 

0.(8(058 

-0.(888)9 

0.(88)49 

0.00055 

0.00002 

0.0005S 

1.4 

0.00050 

0.(88(05 

0.00059 




'..8 

0.(88)52 

-0.(88)11 

0.00041 



0.0004S 

2.52 

0.(88)40 

0.(88)12 

0.(8)037 

0.(88)47 

-0.(888)5 

0.00042 

2.7 

0.00045 

-0.00011 

0.00034 




j l’5K VSO! 502B (Revised) 



estimated 



Parametric 


M 


‘S, 

AC '-„ 

A S» 

P 

A \ 

A \ 

1.2 

0.00055 

0.(88)14 

0. 00009 

0.0(8)52 

0.00021 

0.00073 

1.4 

0.00055 

0. 0(8)52 

0.00085 




1.8 

0.00050 

0.00011 

0.0000 1 



0.00053 

2.52 

0. 0(8)40 

0.00000 

0.188)40 

0.00047 

-0.00001 

0,00046 

2.7 

0.00045 

-0.00004 

0.00039 
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TABLE 23. - Concluded 


GE 21/J10 B1 


Estimated 

— 

Parametric 

M 

AC »r 

A \ 

A S> 

% 



l.Z 

0.0004S 

0.00161 

0.00206 

0.00019 

0.00153 

0.00202 

1.4 

0.00044 

0.00164 

0.002% 




1.8 

0.00041 

0.00097 

0.00138 



0.00123 

2.32 

0.00038 

0.00060 

0.00098 

0.00044 

0.00051 

0.00095 

2.7 

0.00035 

0.00043 

0.00078 




a 21/Jll B3 

Estimated 

Pamctric 

M 


‘S, 


A S 

■t 


iC »o 

1.2 

0.00048 

0.00235 

0.00283 

0.00049 

0.00227 

0.00276 

1.4 

0.00047 

0.00223 

0.00270 




1.8 

0.00044 

0.00140 

0.00184 



0.00157 

2.32 

0.00040 

0.00083 

0.00123 

0.00044 

0.00071 

0.00115 

2.7 

0.00038 

0.00061 

0.00099 







P5K VLE 11 

:C 




Estimated 

Pamctric 

H 

A \ 

A \ 

A \ 


Ac -i, 

A \ 

1.2 

0.0006C 

0.00028 

0.00088 

0.00057 

0.00032 

0.00089 

1.4 

0.00058 

0.00043 

0.00101 




1.8 

0.000S4 

0.00017 

0.00071 



0.00063 

2.32 

0. 00050 

0.00004 

0. IKK) 54 

0.00050 

0.00005 

0.00055 

2.7 

0.0004b 

0.00000 

0.00046 
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(2) Climb - Maximum power climb to cruise altitude and mach number 

(3) Cruise - Subsonic cruise at mach number and altitude for best range 

(4) Climb - Maximum power climb and accelerate to cruise altitude and 
mach number 

(5) Cruise - Supersonic cruise at mach 2.32 at altitude for best cruise 

range 

(o) Descend - Descend and decelerate to 457 m (1500 feet) altitude 
using idle power 

(7) Approach and land - Descend to sea level using idle power 

(8) Taxi - 5 minutes at idle power 

(9) Reserve - Fuel reserves for an alternate airport located 463 km 
(250 n mi) from destination airport 

For each of the four selected propulsion systems, a "basepoint" airplane 
was developed. This basepoint is similar to the basepoint described earlier, 
in that it is the NASA reference airplane but with the selected propulsion sys- 
tem installed. This airplane is then scaled using the Vehicle Sizing and Per- 
formance Evaluation Program (VSPEP) to yield a baseline airplane. Again, this 
baseline is similar to the baseline described earlier, in that it is the mini- 
mum gross weight airplane that meets all specified performance requirements . 

To gain some insight as to what could be expected for each of the four re- 
sized airplanes, an engine scale trade was performed on each. The results of 
these trades are presented in figures 44 and 45 which show design mission range 
and balanced field length versus engine airflow at sea level static conditions. 
The gross wing area was maintained constant at 1022 sq m (10 996 sq ft) and 
airplane gross weight was maintained at that for the basepoint with the first 
selected propulsion system installed or 377 632 kg (744 350 lb). Results of 
these trades indicate that the airplanes having the GE21 engines could be ex- 
pected to result in a higher gross weight to meet the design requirements . In 
addition the VCE 112C and the GE21/J11 engined airplanes can be expected to pay 
an engine size penalty to meet the balanced field length requirement thereby 
increasing the gross weight of those vehicles. 

In the case of the first selected propulsion system, which uses a revised 
VSCE 502B engine installation, sizing was performed maintaining takeoff thrust - 
to-weight (based on sea-level static installed thrust) at 3.13 n/kg (0.32 lb/lb). 
Wingloading was then varied and the VSPEP program allowed to search for the 
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TOGW = 337 640 kg (744 350 lb) 
Sj^gp * 969 sq m (9969 sq ft) 



Sea level static airflow - kg/sec (lb/sec) 


Figure 45. - Balanced field length versus engine size. 
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gross weight yielding 7408 km (4000 n mi) on the design mission. These results 
were then cross-plotted to obtain the wingloading for a minimum gross weight 
baseline airplane. Plots of gross weight and balanced field length versus 
wingloading (based on gross wing area) 'or a family of airplanes having a design 
range of 7408 km (4000 n mi) are shown in figure 46. All other performance 
requirements were easily met. 

For the last three selected propulsion systems , it was originally intended 
to maintain thrust- to- weight at 3.13 n/kg (0.32 Ib/lb) and wing loading at that 
value obtained for the VSCE 502B baseline, which is 345 kg/sq m (70.7 lb/sq f t), 
and si^tly scale gross weight to yield the required design range. However, dut 
to considerable differences in thrust lapse rate at takeoff for each engine, 
this method would not suffice to maintain balanced f ; eld length near the required 
distance. For this reason, the thrust -to-weight was varied in each of the last 
three cases while maintaining wingloading constant at 345 kg/sq m. Plots of 
gross weight and balanced field length versus thrust- to- weight ratio are 
presented in figures 47 through 49 for families of airplanes having a design 
mission range of 7408 km (4000 n mi.). 

Airplane characteristics for both the basepoint and the resized baseline 
vehicles are presented in tables 24 through 27 for each of the fair selected 
propulsion systems. Design mission summaries are shown for each baseline air- 
plane in tables 28 through 35. 

Performance for each baseline airplane was also calculated for the economic 
mission. Payload for this mission was 55 percent of that carried on the design 
mission or 15 256 kg (33 565 lb). Fuel was then off-loaded from the baseline 
to yield 4630 km (2500 n rai.) on the economic mission. Characteristics of each 
of the four baseline airplanes on the economic mission are shown in tables 36 
through 59. 

The results of the economic mission were used to conpute direct operating 
cost (DOC) and return on investment (ROI) for the VSCE 5026 and the VCfc 112C. 

The DOC and ROI calculations were supplied by Pratt and Whitney Aircraft and 
are presented in table 40. Input to the DOC and ROI calculations included the 
following airplane and mission data supplied by Rockwell: airframe weight, 

engine weight, mission fuel, TOQf, block time, and engine design thrust. These 
data are included for each of the four airplanes in tables 36 through 39. 

F.conomic ...ission summaries for each baseline airplane are presented in tables 
41 through 48. 

Sensitivity method verification.- To verify that the sensitivity method 
of determining aircraft takeoff gross weights is valid for preliminary studies, 
takeoff gross weights were estimated using the sensitivities for the four 
engines. These results were then compared with the results of the detailed 
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Balanced field length, m x l(f (ft x 10 ) Takeoff gross weight, kg x 10 (lb x 10 



Wing loading kg/sq m (Ib/sq ft) 


Figure 46.- Sizing with VSCE 502B engines. 
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Balanced field length, m x 10 (ft x 10 *) Takeoff gross weight, kg x 10’ s (lb x 10 



Thrust -to -weight n/kg (Ib/lb) 
Figure 47. - Sizing with VCE 112C engines. 
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Balanced field length, m x 10‘ 3 (ft x 10* 3 ) Takeoff grow weight, kg x 10' 5 (lb x 10 
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TABLE 24. - AIRPLANE CHARACTERISTICS WITH REFINED VSCE 502B ENGINES 


Engine airflow, kg/sec (lb/ sec) 

Thrust -to -weight, n/kg (lb/lb) 

Reference wing area, sq m (sq ft) 

Gross wing area sq m (sq ft) 

Wing loading (gross area), kg/sq m (lb/.?q ft) 

TOGW, kg (lb) 

Fuel weight, kg (lb) 

Max wing fuel, kg (lb) 

Design range, km (n mi) 

Eng out range, km (n mi) 

Normal T.O. dist m (ft) 

Bal field length, m (ft) 

Thrust -to -drag 2.32 m/18 300 m (60 000 ft) 
Thrust -to -drag 1.2 m/climb 

Initial cruise L/D 

Initial cruise SFC (installed) kg/hr/daN 
(lb/hr/lb) 


Basepoint 

Baseline 

408 

( 

900) 

387 

( 

854) 

3.13 

( 

0.3199) 

3.13 

( 

0.320) 

926 

( 

9969) 

841 

( 

9049) 

1022 

( 

10 996) 

927 

( 

9981) 

331 

( 

67.69) 

345 

( 

70.71) 

337 632 

(744 350) 

320 141 

(705 790) 

1S8 683 

(349 835) 

148 772 

(327 986) 

208 222 

(459 050) 

180 061 

(396 966) 

7610 

( 

4109) 

7410 

( 

4001) 

6553 

( 

3538) 

6384 

( 

3447) 

2635 

( 

8645) 

2741 

( 

8991) 

3076 

( 

10 092) 

3197 

( 10 489) 



1.921 



1.899 



1.937 



1.930 



9.626 



9.566 

1.391 

( 

1.364) 

1.389 

( 

1.362) 






TAB Li; 25. - AIRPLANI (IftRACTERISTICS WITH VCE 1 12C ENGINES 




Basepoint 

Baseline 

Engine airflow, kg/ sec (lb/sec) 


408 

( 900) 

600 

( 1323) 

Thrust -to-weight, N/kg (lb/lb) 


2.94 

( .3004) 

3.64 

( .372) 

Reference wing area, sq m (sq ft) 


926 

( 9969) 

1057 

( 11 380) 

Gross wing area, sq m (sq ft) 


1 022 

( 10 996) 

1166 

( 12 552) 

Wing loading (gross area) , kg/sq m 

( lb/ sq ft) 

332 

( 68.02) 

345 

( 70.7) 

TOGW, kg (lb) 


339 272 

(747 966) 

402 618 

(887 622) 

Fuel weight, kg (lb) 


158 683 

(349 835) 

190 934 

(420 938) 

Max wing fuel, kg (lb) 


208 222 

(459 050) 

253 954 

(559 863) 

Design range, km (n mi) 


7419 

( 4006) 

7410 

( 4001) ' 

Eng out range, km (n mi) 


6504 

( 3512) 

6434 

( 3472) 

Normal T.O. dist m (ft) 


3341 

( 10 960) 

2715 

( 8907) 

Bal field length, m (ft) 


3957 

( 12 982) 

3206 

( 10 519) 

Thrust -to-drag at 2.32 M/18 300 m 

(60 000 ft) 


1.220 


1.483 

Thrust -to -drag at 1.2 M/climb 



1.730 


2.050 

Initial cruise L/D 



9.605 


9.601 

Initial cruise SFC (installed) kg/hr/daN 

1.470 

1.442 

1.491 

( 1.462) 

(lb/hr/ lb) 











TABLE 26.- AIRPLANE CHARACTERISTICS WITH GE21/J10 B1 ENGINES 


Engine airflow, kg/sec (lb/sec) 
Thrust-to-weight, n/kg (lb/lb) 

Reference wing area, sq m (sq ft) 

Gross wing area, sq m (sq ft) 

Wing loading (gross area), kg/sq m (lb/sq ft) 

TOGW, kg (lb) 

Fuel weight, kg (lb) 

Max wing fuel, kg (lb 

Design range, km (n mi) 

Eng out range, km (n mi) 

Normal, T.O. dist m (ft) 

Bal field length, m (ft) 

Thrust -to -drag at 2.32 M/18 300 m (60 000 ft) 
Thrust -to- drag at 1.2 M/climb 

Initial cruise L/D 

Initial cruise SFC (installed) kg/hr/daN 
(lb/hr/lb) 


Basepoint 

Baseline 

340 

( 750) 

624 

( 1375) 

2.71 

( .2769) 

3.31 

( .338) 

926 

( 9969) 

1351 

( 14 540) 

1022 

( 10 996) 

1490 

( 16 038) 

335 

( 68.68 

345 

( 70.71) 

342 556 

(755 206) 

514 441 

(1 134 149) 

158 683 

(349 835) 

255 446 

(563 161) 

208 222 

(459 050) 

366 785 

(808 622) 

5045 

( 2724) 

7410 

( 4001) 

4939 

( 2667) 

6295 

( 3399) 

3348 

( 10' 983) 

2736 

( 8977) 

3909 

( 12 823) 

3197 

( 10 490) 


1.084 


1.286 


1.235 


1.694 


8.923 


9.514 

1.568 

( 1.538) 

1.457 

( 1.429) 





TABLE 27.- AIRPLANE CHARACTERISTICS WITH GE21/J11 B3 ENGINES 



Basepoint 

Baseline 

Engine airflow, kg/scc (lb/ sec) 

381 

( 840) 

1093 

( 2409) 

Thrust-to-weight, n/kg (lb/lb) 

2.16 

( .2210) 

3.33 

( .340) 

Reference wing area, sq m (sq ft) 

926 

( 9969) 

1653 

( 17 787) 

Gross wing area, sq m (sq ft) 

1022 

( 10 996) 

1823 

( 19 619) 

Wing loading (gross area) , kg/sq m (lb/sq ft) 

330 

( 67.68) 

345 

( 70.71) 

TOGW, kg (lb) 

337 557 

(744 186) 

6/ -..6 (1 387 359) 

Fuel weight, kg (lb) 

158 683 

(349 835) 

312 330 

(688 569) 

Max wing fuel, kg (lb) 

204 139 

(450 050) 

496 237 (1 094 016) 

Design range, km (n mi) 



7412 

4002) 

Eng out range, km (n mi) 



6245 

( 3372) 

Normal T.O. dist m (ft)* 



2706 

( 8878) 

Bal field length, m (ft) 



3205 

( 10 515) 

Thrust- to-drag at 2.32 M/18 300 m (60 000 ft) 




1.511 

Thrust -to-drag at 1.2 M/climb 




2.078 

Initial cruise L/D 




9.412 

Initial cruise SFC (installed) kg/hr/daN 



1.447 

( 1.419) 

(lb/hr/lb) 





*T.O. distance using maximum power. This results in noise 

levels below FAR 36. 





TABLE 28. - REFINED VSCE 5Q2B BASELINE DESIOt MISSION SUNMARY, INTERNATIONAL UNITS 



TOT 

OPERATION 

WEIGHT 

ALTITUDE MACH NO. 

PUEL USED 

TIME TOTAL TIME 

RANGE TOTAL RANGE 


NO. 


KG. 

METERS 


KG. 

MIN. 

MIN. 

w. 

KM. 

DO 


iNlVlAL MClGHT 

320140.9 








g| 

1 

MU ( TO 

314337.4 

0.0 

0.299 

3803.4 

10.000 

10.000 

0.0 

0.0 

2 

CL T015C 

314196.4 

497.2 

0.900 

2141.1 

1.337 

11.337 

10.T2 

10.72 

Sg 


EE1=S3C 

290090.4 

17221.9 

2.320 

24146.0 

12.525 

23.862 

207.30 

290.02 

h » 

K S3 

4 

CRUISE 

238697.2 

18402.3 

2.320 

91393.2 

01 .462 

105.324 

3407.20 

3709.22 

5 

CRUISE . 

149807.7 

19796.2 

2.320 

42769.4 

80.963 

186.280 

3306.39 

7091.57 


<T 

OESCEND 

194324.4 

497.2 

0.9 00 

1563.4 

17.600 

203.967 

305.49 

7397.02 


7 _ 

OESHLANO 

194084.7 

0.0 

0.300 

239.6 

1 .607 

209.574 

13.42 

7410.45 


8 

TAXI-ALL 

193449.6 

0.0 

0.0 

639.2 

5.Q00 

210.576 

0.0 

7410.45 


___ 

5PCT ALL 

1071 14.9 

0.0 

0.0 

6334.6 

0.0 

210.574 

0.0 

7410.49 


_|©_ 

CLTO 15C 

186929.2 

497.2 

0.900 

505.0 

0.739 

211.313 

5.66 

7414.11 


11 

CL 8— ACC 

181900.4 

11990.1 

0.900 

4628.8 

8.236 

219.549 

124.00 

7540.20 


IF" 

CRUISE 

174047.6 

12092.4 

0.900 

2852.4 

13.921 

233.071 

219.39 

7759.99 


13 

DESCEND 

170109.9 

3048.0 

0.466 

898.9 

9.375 

242.446 

119.41 

7070.00 

' 

14 

LOITER 

171402.2 

3040.0 

0.494 

6287.3 

30.000 

272.446 

0.0 

7879.00 


"is 

OESH.ANO 

171399.3 

0.0 

0.300 

542.8 

3.077 

276. 323 

35.00 

7414.09 


TOTAL FUEL USED* 140781.2 



TABLE 29.- REFINED VSCE 5Q2B BASELINE DESIGN MISSION SUMMARY, ENGLISH UNITS 


NASA CR-132374 AST KF. CONFIG. WITH VSCt 3028-02 ENGIN8S 
LEWIS BASELINE T/W-.32 W/S*Tt. ROES ■4000 NM 


LEG 

NO. 

OPERATION 

WEIGHT 

POUNDS 

ALTITUDE MACH NO. 
FEET 

FUEL USED 
ROUNDS 

TIME TOTAL TIME 
MIN. MIN. 

range fomwam 

N. M. N. N. • 


INITIAL WEIGHT 

705789.4 








1 

MU t TO 

697404.8 

o 

• 

© 

0.299 

8385.1 

10.000 

10.000 

t 

0.0 

0.0 

2 

CL T015C 

602684.6 

1500.0 

0.500 

4720.2 

1.337 

11.337 

5.79 

3.79 

3 

CL B -ACC 

639AM .8 

56500.9 

2.320 

53232.7 

12.525 

23.1*2 

- Tir. iT" 

1*0.92 

A 

CRUISE 

526149.2 

60374.9 

2.320 

113302.6 

81 .442 

105.324 

Mlf.TL 

. 2000.42 

5 

CRUISE 

A31858.7 

64816.8 

2.320 

94290.5 

10.963 

186.200 

1028.45 

3029.00 

6 

DESCEND 

428411.9 

1500.0 

0.5 00 

3446.7 

17.480 

20S.94T 

1*4.93 

3994.00 

7 

OES— LAND 

627883.7 

0.0 

0.300 

528.2 

1 .407 

205.574 

T.n _ 

400l.lt 

b 

TAXI-ALL 

626663.3 

0.0 

0.0 

1400.4 

5.000 

210.574 

0.0 

4001.23 

V 

5PCT ALL 

*12517.6 

0.0 

0.0 

13965.4 

0.0 

210.574 

0.0 

*001.29 

1C 

CL TO 15C 

4U226.6 

1500.0 

0.500 

1291.4 

0.739 

211.313 

3.0* 

4004.UL 

11 

CL 0 -ACC 

601021.7 

39337.6 

0.900 

10204.8 

8.236 

219.549 

*7.00 

*071.21 

12 

CRUISE 

396733.2 

39673.3 

0.900 

6288.4 

13.521 

233.071 

118.44 

*109.77 

13 

DESCEND 

392640.6 

10000.0 

0.466 

1892.6 

9.375 

242.44* 

64.41 

3t9±m2iL 

lA 

LOITER 

378979.5 

10000.0 

0.4 54 

13841.1 

30.000 

272.44* 

0.0 

*294.2* 

15 

OES-LAND 

377782.7 

0.0 

0.300 

1194.7 

3.877 

29*. 323 

18.93 

4279.17 




TOTAL FUEL 

USED* 

328007.1 







TABLE 30. - VCE 112C BASELINE DESIGN MISSION SUMMARY, INTERNATIONAL UNITS 



LEG 

NO. 

OPERATION 

WEIGHT 

KG. 

ALTITUDE MACH NO. 
METERS 

FUEL USED 
KG. 

TIME TOTAL TIME 

MIN. MIN. 

RANGE TOTAL RANGE 

KM. KM. 

Og 

S| 

1 

INITIAL WEIGHT 
MU ( TO 

40261*. 3 
390212.2 

0.0 

0.299 

4406.2 

10.000 

10.000 

0.0 

0.0 

2 

Cl TO ISC 

396092.0 

457.2 

0.500 

2120.1 

1.334 

11.334 

10.60 

10.60 

|S 

3 

CIB-ACC 

370960.4 

18063.6 

2.320 

25131.6 

t2.99l 

24.324 

320.56 

331.16 

fe O 

4 

CAUISE 

301835.9 

19678.5 

2.320 

69124.4 

80.672 

104.996 

3374.17 

3705.33 


5 

CRUISE 

244642.1 

21247.6 

2.320 

57193.9 

81.124 

186.120 

3397.44 

7102.82 


6 

OE SCENO 

242332.6 

457.2 

0.500 

2309.5 

16.797 

202.917 

295.10 

7347.91 


7 

OES-LANO 

241980.9 

0.0 

0.300 

351.7 

1.523 

204.440 

12.82 

7410.73 


• 

TAXI-ALL 

240996.6 

0.0 

0.0 

984.3 

5.000 

2C9.440 

0.0 

7410.73 


9 

5PCT ALL 

232915.6 

0.0 

0.0 

8081.1 

0.0 

209.440 

0.0 

7410.73 


10 

CL TO 15C 

232263.8 

457.2 

0.500 

651.7 

0.619 

210.059 

4.78 

7415.51 


n 

CL 6-ACC 

226368.2 

12684.3 

0.900 

5895.6 

9.457 

219.516 

130.89 

7546.40 


12 

CRUISE 

222708.3 

12798.6 

0.900 

3660.0 

13.063 

232.579 

211.96 

7758.35 


13 

DESCEND 

221396.2 

3048.0 

0.442 

1312.0 

9.346 

241.925 

120.16 

7878.52 


14 

L31TFR 

212458.7 

3048.0 

0.432 

8937.6 

30.000 

271.925 

0.9 

7878.52 


15 

OES-LANO 

211690.4 

0.0 

0.300 

768.3 

3.533 

275.458 

31.74 

7910.25 





total FUEL 

USEO« 

190927.5 







TABLE 31. - VCE 112C BASELINE DESIGN MISSION SUWARY, ENGLISH UNITS 


LEG 

OPERATION 

WEIGHT 

ALTITUOE MACH NO. 

EUEL USED 

TIME TOTAL TIME 

RANGE 

TOTAL RANGE 

MO . 


POUNOS 

ffrr 


POUNOS 

MIN. 

MIN. 

N. M . 

N. M . 


INITIAL WEIGHT 

887621.6 








1 

wu C TO 

877907.6 

0.0 

0.299 

9713.9 

10.000 

10.000 

0.0 

0.0 

2 

CL T015C 

873233.5 

1500.0 

0.500 

4674.1 

1.334 

11.334 

5.73 

5.73 

3 

CL 8-ACC 

817827.9 

592 6 3.7 

2.320 

55405.6 

12.991 

24.324 

173.09 

178.81 

4 

CRUISE 

665434.4 

64562.1 

2.320 

152393.4 

80.672 

104.996 

1821.87 

2000.66 

5 

CRUISE 

539343.6 

69710.1 

2.320 

126090.9 

81.124 

186.120 

1834.47 

3635.15 

6 

DESCEND 

534252.1 

1500.0 

0.500 

5091.5 

16.797 

202.917 

159.34 

3994.49 

7 

DES-LAND 

533476.7 

0.0 

0.300 

775.3 

1.523 

204.440 

6.92 

4C01.4I 

8 

TAXI-ALL 

531306.7 

0.0 

0.0 

2170.1 

5.000 

209.440 

0.0 

4001.41 

9 

5PCT ALL 

513490.9 

0.0 

0.0 

17815.7 

0.0 

20«.440 

0.0 

4001.41 

10 

CL TO 1 5C 

512054.2 

1500.0 

0.500 

1436.7 

0.619 

.10.059 

2.58 

4003.99 

11 

CL 8-ACC 

499056.6 

41515.1 

0.900 

12997.6 

9.457 

219.516 

70.67 

4074.66 

12 

CRUISE 

490987.8 

41990.3 

0.900 

8068.8 

13.063 

232.579 

114.45 

4189. 11 

13 

OESCENO 

488095.2 

10000.0 

0.442 

2892.6 

9.346 

241.925 

64.88 

4253.96 

14 

LOITER 

468391.3 

10000.0 

0.432 

19703.9 

30.000 

271.925 

0.0 

4253.98 

15 

OES-LANO 

466697.6 

0.0 

0.300 

1693.7 

3.533 

275.458 

17.14 

4271.12 




TOTAL EUEL 

USED- 

420924.0 







TABLE 32. - GE21/J10 


,r ELINE DESIGN MISSION SUMMARY, INTERNATIONAL UNITS 


III. :>tKATlL)M 

M . 




INITIAL WEIGHT 



i 

WU t Til 


81 


LL TC16L 


9 

Ut -Af.C 




CFblft 

y o 
C E» 

•n 

CRUISE 



*> 

1 fcSClNL 


7 

bfcS-LAND 


! 

taxi-all 


«# 

5PC T ALL 


1 

CLTu 15C 


1 

CLl-Aft 


1< 

CRUlSf 


1 4 

bfcSCEMh 


I- 

lUlTeh 


15 

I* S-LANl> 


WLJUHT ALTITUDE HaCM Nb. Pl*H USED 


KO. 

METERS 


KG. 

11 **4*1. 3 




V~9639.1 

C • U 

0.299 

4802.2 


457.2 

0.5 v- 

3016. & 

*»blo66 *2 

17538.4 

2.320 

44954 .'t 

38012*1.0 

18713.0 

2 .3 21 

81542.7 

3-8072.2 

19901.2 

2.3 20 

71453.3 

3:6-96.3 

457.2 

0 .5 0 ) 

2575.9 

3.570: .4 

o.:- 

0.307 

395 . t 

3f 4o77.< 

c.c 

O.J 

1323.2 

29ol89 .r 

*• • ' 

0.0 

10488.2 

293079 .o 

457.2 

0.5 00 

1109.5 

279368 .5 

17288 .S 

0.950 

13711.0 

277377.4 

1 726o.o 

0.9 M 

1991.1 

?75oo6»o 

3048.0 

0.3*8 

1710.8 

259781.0 

3048.0 

0.3 50 

15885.0 

?5B9Pe.7 

-j.a 

0.3 00 

792.9 


TIME TOTAL TIME RMIGE TOTAL RANGE 


MIN. 

MIN. 

KM. 

KM. 

10.000 

10.000 

0.0 

0.0 

1.280 

1 1. 280 

10.16 

10.16 

16.868 

30. 148 

477. r 15 

487.31 

79.936 

107.084 

3217.80 

3705.19 

80.827 

187.911 

3380.66 

7085.85 

16.055 

205.995 

310.60 

7396.45 

1.649 

207.614 

13.94 

7410.38 

5.000 

212.614 

0.0 

7410.38 

0.0 

212.614 

6.6 

7410.38 

0.492 

213. 106 

3.84 

7414.22 

10.537 

223.612 

199.12 

7604.34 

6.362 

229.974 

108.96 

7713.30 

12.115 

242.089 

163.93 

7877.23 

30.000 

272.089 

0.0 

7877.23 

3.480 

275.569 

30.86 

7908.09 


TOTAL FULL U&EO- 


259452.1 



TABLE 33. - GE21/J10 B1 BASELINE DESIGN MISSION SUMMARY, ENGLISH UNITS 


LFG 

URFRATIUM 

NE IGHT 

ALT I TUI F 1 

RaO no. 

FULL USE! 

TIME TOTAL TIME 

RANGE TOTAL RANGE 

NO. 


ROUNDS 

RttT 


ROUNDS 

MIN. 

MIN. 

N * M » 

N. M. 


TRTTTal we T(tRT 

113414V. 0 








1 

WU L TO 

1123562.? 

C. j 

0.299 

10987.0 

16.000 

10. 000 

0.3 

b.O 

2 

CL TU15C 

1116911.3 

I5C 

a .sc: 

6651 .0 

1.28C 

11.283 

5.48 

9.4b 

“5“ 

'ClS-aK 

T5iTb$4 _ .l 

*“YTV40.'i 

2‘. 3 20 

99100.7 

lb .868 

SO. 148 

297.64 

263.12 

4 

cruise 

836^33.*. 

61384.7 

2.3 2? 

I7977C.4 

76.436 

107.084 

1737.49 

2000.61 

5 

CRUISE 

6S05&5.7 

69292.0 

2.320 

197527.7 

80.82? 

187.911 

1825.38 

3829.99 

6 

descent 

4748 26. V 

lSu'.O 

0.5 u> 

5678.4 

18.055 

205.969 

167.71 

3993.69 

7 

DFS-LANO 

473994.2 

o.O 

0.30* 

072.7 

1.649 

2 .614 

7.52 

4001.22 

C 

TAXI-ALL 

071498.4 

u.O 

0.0 

2295.7 

5.000 

212.614 

0.0 

4001.22 


spcYall 

648579.9 

0.0 

0.0 

23122.6 

0 ,f> 

212.614 

0.0 

4f>',1.22 


CLTO ISC 

6461 29 .4 

1500. v 

o.scii 

2446.0 

0.492 

219. 106 

2.07 

4003.29 

li 

CL8-ACC 

615902.2 

96721.8 

0.9 50 

30227.7 

10 »9o 7 

223.612 

102.66 

4109.95 


CRUISE 

6ll5l?.6 

§o7|4.7 

i.9*>0 

4389.6 

6.362 

229.974 

5b. 83 

4164. 77 

b 

w-aCEND 

607741 

lOOOfl.n 

0.344 

3771.6 

12.115 

242.689 

88.51 

4253.29 

14 

LOITER 

572720. S 

10000.0 

0.350 

35320.5 

30 .'‘00 

272.089 

t.G 

4253.29 

~n> 

oIPIand 

570972.6 

J.J 

0.3w 1 

1747 ,9 

3.4b«. 

275. 969 

lb. 66 

4269.95 


TOTAL FUEL UStl» 


563176.4 



TABLE 34. - GE21/J11 B3 BASELINE 


LEG OPERATION WEIGHT ALTITUDE 

NO* KG. METERS 




INITIAL WEIGHT 

629203.4 



i 

WU I TO 

622324.4 

0.0 

o S 

2 

CL T019C 

619027.4 

437.2 

H'J G-: 
, 

3 

Cl 6 -ACC 

573132.9 

18712.9 






ijri ^ 

4 

CRUISE 

A 67000.2 

20118.3 

ft*' 

5 

CRUISE 

378848 . 1 

21470.4 

i:^ ( 1 

6 

OESCENO 

374795.3 

497.2 


T 

DCS -LAND 

3741 72.1 

0.0 


a 

TAXI-ALL 

372379.6 

0.0 


9 

9PCT ALL 

359560.4 

0.0 


10 

CLTO 19C 

358252.1 

437.2 


11 

Cl B -ACC 

347471.4 

14392.7 


12 

CRUISE 

340263.0 

14935.0 


13 

DESCEND 

337646.0 

3041,0 


14 

LOITER 

318150.9 

3048.0 


19 

OCS-iAND 

316980.6 

0.0 


TOTAL 


6 


/ 

i 

DESIGN MJSSIOn/sUKWARY, INTERNATIONAL UNITS 


iCH NO. PUCL ukeo 

kg; 


TINE 

MIN. 


TOTAL TIN! 
MIN. 




USED* 3 1 23 1 N . 8 



TABLE 35. - CE21/J11 B3 BASELINE DESIGN MISSION SUNMAKY, ENGLISH UNITS 


TIG 

OPERATION 

WEIGHT 

ALTITUCf MACH NO. 

RUEL USED 

TINE TOTAL TINE 

RANGE TOTAL RANOf 

NO. 


ROUNDS 

REIT 


ROUNDS 

MIN. 

MIN. 

N« M. 

N. M. 


\ 

ISiTSST.O 








. A 

JWL6J0 

1ST 1441.0 

r.o 

0.244 

19561.0 

10.000 

10.000 

0.0 

0.0 

a 

CL T01»C 

1364727.0 

1500.0 

0.500 

7244.0 

0.414 

10.414 

3.00 

1.00 

T" 

TERrar 

T3 63516.0 

61544.0 

2.520 

101156.0 

12.147 

23.911 

175.40 

177.40 

4 

.CRUISE 

1029555.3 

66009.1 

2.520 

235441.7 

10.733 

104.944 

1123.34 

2000.74 

S 

CRUISE 

835217.1 

70441.0 

2.5 20 

14ST57.2 

•1.401 

185.445 

1141.30 

3842.04 


olittND "" 

826282.2 

1500.0 

0.500 

•434.4 

16.191 

202.046 

192.70 

3444.74 

_,X. 

JCS-LMO 

824908.2 

0.0 

0.500 

1574. 1 

1.413 

205.974 

4.73 

4001.47 

• 

TAXI-ALL 

820956.6 

0.0 

0.0 

3492.4 

9.000 

201.974 

0.0 

4001.47 

1" 

wm~ 

79269<*. 9 

0.0 

0.0 

21241.7 

0.0 

208.974 

0.0 

4001.47 

10 _ 

CL TO 15 C 

789810.6 

1500.0 

0.500 

2SS4.4 

0.340 

209.414 

l.M 

4002.06 

11 

CL ■ -ACC 

766043 . 2 

47084.1 

0.490 

23766.3 

3.458 

214.397 

47.34 

4030.20 

TT" 

Hulii 

~ 750151.6 

47416.4 

0.490 

19S42.T 

14.096 

221.413 

124.44 

4100. IS 

u 

OESCENO 

76.4 3 3 2.) 

10000.0 

0.520 

5748.4 

10.391 

251.7*6 

72.67 

4292.19 

14 

LOITER 

701 602 .7 

19000.0 

0.504 

42*80.4 

50.000 

26S.764 

0.0 

4212.09 

19 HI-t«fcD 

698822.7 

0.0 

0.500 

2SE0.1 

2.421 

271.616 

13.01 

4266.6* 




TOTAL RUEL 

USED* 

688536.3 







TABLE 36. - BOQNONIC MISSION CHARACTERISTICS 
WITH REFINED VSCE S02B ENGINES 



Baseline 

TOOf. kg (lb) 

254 141 

(560 284) 

Fuel weight, kg (lb) 

95 229 

(209 943) 

Design range, km (nmi) 

4630 

( 2500) 

Far 36 T.O. dist ■ (ft) 

1753 

( 5750) 

Bal field length, m (ft) 

2081 

( 6826) 

Airframe weight, kg (lb) 

112 768 

(248 611) 

Engine weight, kg (lb) 

23 128 

( 50 988) 

Mission fuel wei^it, kg (lb) 

76 501 

(168 655) 

Block tine, nin (min) 


( 167) 

Engine thrust, n (lb) 

251 160 

( 56 463) 


TABLE 37. - EOONGMIC MISSION CHARACTERISTICS 
WITH VCE 11 2C ENGINES 



Baseline 

ff 

H- 

cr 

320 966 

(707 608) 

Fuel weight, kg (lb) 

121 739 

(268 390) * 

Range, km (nai) 

4650 

( 2500) 

Far, 36 T.O. dist n (ft) 

1724 

( 5656) 

Bal field length, ra (ft) 

2070 

( 6794) 

Airframe weight, kg (lb) 

138 186 

(304 645) 

Engine weight, kg (lb) 

37 932 

( 83 626) 

Mission fuel weight, kg (lb) 

96 639 

(213 048) 

Block time, min (min) 


( 167) 

Engine thrust, n (lb) 

367 178 

( 82 549) 




TABLE 38. - BCGHOMIC MISSION CHARACTERISTICS 
WITH GE21/J10 B1 ENGINES 


* 

Baseline 

TOGW, kg (lb) 

417 609 

(920 670) 

Fuel weight, kg (lb) 

171 070 

(377 144) 

Design range, km (mu) 

4630 

( 2500) 

Far 36 T.O. dist a (ft) 

1816 

( 5959) 

Bal field length, ■ (ft) 

2159 

( 7083) 

Airframe weight, kg (lb) 

166 931 

(368 020) 

Engine weight , kg (lb) 

56 360 

(124 253) 

Mission fuel weight, kg (lb) 

130 675 

(288 088) 

Block tiae, win (ain) 


( 169) 

Engine thrust, n (lb) 

426 300 

( 95 836) 


TABLE 39. - BCONDNIC MISSION CHARACTERISTICS 
WITH GE21/J11 B3 ENGINES 



Baseline 

TOGW, kg (lb) 

519 218 

(1 144 680) 

Fuel weight, kg (lb) 

214 709 

(473 353) . 

Design range, kn (rani) 

4630 

( 2500) 

Far 36 T.O. dist, a (ft) 

1855 

( 6087) 

Bal field length, a (ft) 

2232 

( 7323) 

Airframe weight, kg (lb) 

207 300 

(457 018) 

Engine weight, kg (lb) 

73 812 

(162 728) 

Mission fuel weight, kg (lb) 

162 334 

(357 886) 

Block tiae, min (min) 


( 166) 

Engine thrust, n (lb) 

524 557 

(117 925) 


122 





TABLE 40. - DOC, IOC, AND HOI 



123 









ORIGINAL PAGE IS 
OF POOR QUALITY 


TABLE 41 


VSCE 502 B BASELINE ECONOMIC MISSION SUMMARY, INTERNATIONAL UNITS 


h-» 

K) 


t?r. 

OPERATION 

we IGHT 

ALTITUDE pach nc. 

FUEL USEO 

TIME TOTAL T|Ki 

RANGE TOTAL RANW 

NO. 


K 

M( TM $ 


KG. 

PIN. 

MtN. 

KM. 

RN. 


IT I*L WEIGHT 

21616C.6 








1 

tall f. TO 

250337.2 

0.0 

0.289 

3803.6 

10.000 

10.000 

0.0 


2 

CL/ 1500 

268016.7 

657.2 

0.500 

1262.5 

0.569 

10.569 

6.75 

6.73 

3 

CL IMP 

262163.7 

1C060.1 

0.900 

6931.0 

3.526 

16.173 

63.12 

KIT 

6 

CBU1SP 

230565. 3 

1C6U.0 

0.900 

11598.6 

61.502 

56.176 

665.97 

760.51 

5 

CL/ACCEL 

22C068.7 

15850.7 

2.320 

10696.6 

5.667 

61.561 

162.36 

903.16 

6 

C»IIISE 

150063.8 

2C387.5 

2.320 

6J006.5 

51.576 

163.61? 

' 3411.21 Qli.fi - 

7 

CfSCfNr 

17(512. 8 

657.2 

0.500 

1531.0 

17.303 

160.720 

300.22 

6616.50 

a 

OES/IANC 

175275.1 

0.0 

0.300 

237.7 

1.596 

162.316 

13.65 

6630.03 

q 

TAX I 

177638.8 

0.0 

0.0 

635.2 

5.000 

UT.3U 

o 

. 

O! 

— 4194. SI 

10 

5T PBS. 

173F16.8 

0.0 

0.0 

3525.0 

0.0 

167.316 

0.0 

6630.03 

11 

CL / 1 500 

173258.6 

657.2 

0.500 

525.3 

0.662 

167.975 

5.07 

6635.10 

12 

CL/ACCFL 

16F760. 1 

12687.3 

0.900 

6699.5 

5.596 

176.573 

WfTfl 

— 4T44.44 

13 

C«UISF 

16623*. 6 

12586.5 

0.900 

2550.7 

12.975 

159.565 

210.53 

6975.36 

16 

CFSCFNO 

165367.8 

3C68.0 

0.650 

871.5 

9.597 

199. 165 

122.73 

5096.11 

15 

LOITfP 

15863C. t 

3C66.C 

0.660 

5929.2 

30.000 

229. 165 

~«?r~ 

- swum - 

16 

CES/l AMR 

15(816.6 

0.0 

0.300 

522.3 

3.726 

232.569 

>3.51 

5131.62 




TOTAL FUEL 

UStP. 

95223.7 







TABLE 42. - VSCE 502B BASELINE ECONOMIC MISSION SUMMARY, ENGLISH UNITS 


LEG 

OPERATION 

WEIGHT 

ALTITUDE MACH NC. 

FUEL USEO 

TtMf TOTAL TIME 

RANGE TOTAL RANGE 

NO. 


POUNDS 

FEET 


POUNDS 

MIN. 

MIN. 

N. M. 

•N. M. 


“umui Wight 

56C284. 1 








1 

WU A TO 

551899.1 

0.0 

0.299 

8385.1 

10.000 

10.000 

O.C 

0.0 

2 

CL/ 1800 

549115. 6 

15C0.0 

0.500 

2783.2 

0.849 

10.849 

3.63 

3.63 

3 

~ClTm8 

533825.4 

32439.8 

0.900 

18280.2 

3.524 

14.373 

2 5.98 

29.61 

4 

CRUISE 

508245.4 

34156.7 

0.900 

25570.2 

41.802 

56.174 

370.39 

400.00 

5 

CL /ACC El 

485124.5 

62174.3 

2.320 

23140.9 

5.667 

61.841 

87.66 

487.66 

6~ 

CRUISE 

394928. 5 

66920.9 

2.320 

88 195.6 

81.576 

143.417 

1842.45 

2330.61 

7 

DESCEND 

393553.6 

1500.0 

0.500 

3375.4 

17.303 

160.720 

162.10 

2492.71 

8 

DES/l AND 

343029.4 

0.0 

0.300 

524.1 

1.596 

162.316 

7.26 

2499.97 

9 ~ 

“TAXI ' 

391624.1 

0.0 

0.0 

1400.4 

5.000 

167.316 

0.0 

2499.47 

10 

5« RES. 

383196.3 

0.0 

O 

• 

o 

8432.7 

O 

• 

o 

167.316 

0.0 

2499.97 

11 

CL / 1 SOO 

382038.2 

1500.0 

0.500 

1158. 1 

0.662 

167.978 

2.74 

2502.71 

Ti~ 

“CL/ACCEL 

372118. 6 

41001.6 

0.900 

9919. T 

8.594 

176.573 

70.06 

2572.77 

13 

CRUISE 

366495. 1 

41321.5 

0.900 

5623.4 

12.475 

184.548 

113.67 

2686.45 

14 

OESCENG 

364573.9 

10000.0 

0.450 

1921.2 

9.597 

199. 145 

66.27 

2752.72 

15 

"loiter 

351502. 1 

10300.0 

0.440 

130T1.7 

30.000 

224.145 

0.0 

2752.72 

14 

OES/LANC 

35035C.7 

0.0 

0.300 

1151.4 

3.724 

232.869 

18.09 

2770.81 


TOTAL FUEL USEC* 209933.4 



TABLE 43. - VCE 11 2C BASELINE ECONOMIC MISSION SUMMARY, INTERNATIONAL UNITS 


o? 


LEG 

URFRATION 

WEIGHT 

ALTITUOE MACH NCI. 

FUEL USFi' 

TIHE TOTAL T1HE 

RANGE TOIAL RANGE 

NO. 


KG. 

HETERS 


KG. 

MIN. 

MIN. 

KM. 

KM. 


INITIAL WEIGHT 

32"9ft 5.7 








1 

WO L TO 

316559.6 

3.0 

0.299 

4406.2 

10.900 

10.000 

0.0 

0.0 

2 

CL/1V' 

31529V. 3 

457.2 

:.5 05 

1267.2 

5.895 

U.835 

6.40 

6.40 

j * 

CL I HE 

308794.2 

17G27.C 

0.9 CO 

6505.1 

3.411 

16.217 

65.95 

52.19 

4 

CRUISE 

2929110. 3 

11294.8 

0.9 (VI 

1581%. 9 

62.313 

56.530 

688.46 

740.11 

5 

CL/ AC CEL 

28133 >.** 

20226.5 

2.3 2> 

11649.4 

6.296 

62.786 

195.13 

915.96 

6 

CRUISE 

227 , >2 , j,5 

21726.: 

2.3 25 

514^.5 

80 .944 

143.729 

3392.36 

4128.12 

7 

OtSCFNO 

225662.4 

457.2 

'>.i oo 

2263.1 

16 .444 

169. 179 

289.79 

4618.07 

e 

DES/LANU 

/2S313.* 

o.c 

0 .3 V 

349.4 

1 .5U 

161.689 

12.72 

4610.79 

V 

TAXI 

2243ZS.7 

r.,7 

0.9 

984.3 

5 .OOP. 

1*6.689 

0.0 

4*30.nr 

u 

51 RES. 

2 1444*,. U 

o.c 

0.0 

4831.9 

0.7 

166.689 

0.0 

4610.79 

11 

CL/1*>" , ‘ ' 

211.896. 4 

457.2 

3.5 CG 

6?' *4 

3.579 

167.259 

4.40 

4615.19 

12 

CL/ACCEL 

213244.7 

13078.8 

0.903 

9601 .6 

9.226 

176,685 

127.26 

4762.45' 

13 

CRUISE 

2*4811.7 

13192.4 

0.9 00 

3453. 0 

13.164 

189.669 

213.59 

4976.04 

14 

HESCE NO 

20R4CS.1 

3048.0 

0 .428 

1323.6 

9,483 

199.132 

122.07 

9098.11 

1* 

LOITER 

1444to7.> 

3348.0 

0.4 IF 

8523.8 

3J.00G 

229.132 

0.0 

5098.11 

lo 

OES/LANO 

199226.1 

0.0 

0.33- 

741.2 

3.414 

232.536 

33.45 

5120.55 




TOTAL FUEL 

USED* 

121734.4 









TABLE 44. - VCE 112C BASELINE ECONOMIC MISSION SUMMARY , ENGLISH UNITS 


lO. n3 
c -‘ 3 > 


0 

1 


If 1 

UPERLTICN 

WEIGHT 

ALTITWE MACH NG. 

FULL USFO 

TIME TOTAL T IMF 

RANGE TOTAL RANGE 

NO. 


ROUNfiS 

FIET 


FCUNOS 

MIN. 

MIN. 

N. M. 

N. N. 


INITIAL WEIGHT 









1 

wu t in 


o. 0 

0.244 

4713.4 

IV. COO 

13. jOO 

0.0 

0.0 

? 

CL/ 1500 

6V41 16.1 

15**.*, 

0.5 77 

2778.4 

.7 .474 

IP.8O5 

3.46 

3.66 

3 

CL IMF- 

<iB>774.S 

34R66.5 

: .4 :* 

14141.2 

3.411 

14.217 

24.81 

28.Y7" 

4 

CRUISE 

6454U.1 

37756.3 

0.4 >'• 

34863.7 

42.313 

56.633 

371.73 

600.00 

•> 

CL /AC CEL 

t>?»J2?7 mb 

66360.0 

2.7 20 

26682.4 

6.256 

62.786 

105.37 

505.37 

6 

CRUISF 

577684.7 

71274.5 

2.3 20 

117738.4 

80.444 

143.724 

1831.70 

2337.06 

7 

OESCfND 

447503.4 

1500.0 

0.4 07 

4444. a 

16 .445 

160.175 

136.63 

2643.31 

8 

UES/LANU 

44673* *2 

0.0 

0.3 00 

770.2 

1 .514 

161.684 

6.87 

2500.38 

4 

TAXI 

4V454.C.1 

0.3 

0.0 

2170.1 

5.909 

166.684 

0.0 

2300.38 

10 

•>* RFS. 

4834*7.7 

n." 

0.7 

10662.4 

0.0 

166.684 

0.0 

2300.38 

11 

CL/ 1600 

482693.4 

1500.0 

0.600 

1323.7 

0.570 

167.254 

2.38 

2502.76 

12 

CL/ ACCEL 

470234.4 

42404.3 

0.400 

12344.5 

6.226 

176.485 

6*8.72 

2571.68 

13 

CRUISE 

462654.7 

43282.3 

0.430 

7678.7 

13.164 

184.664 

115.33 

2606.80. 

14 

OFSCENO 

454637.7 

10000.0 

0.4 28 

2419.0 

4.483 

144.132 

65.61 

2752.71 

15 

LOITER 

440842.4 

10000.7 

0.418 

18784.2 

30.000 

224.132 

0.0 

2752.71 

16 

OE S/LAND 

434218.4 

0.0 

9.3 00 

1634.1 

3.404 

232.536 

16.66 

2764.13 




TOTAL FUEL 

USED* 

268340.0 






fsi 
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TABLE 45. - GE21/J10 B1 BASELINE ECONOMIC MISSION SIJNMARY, INTERNATIONAL UNITS 


M 

S 


LSG 

OPERATION 

WEIGHT 

ALTITUDE PACH NC. 

RLEL USED 

TIME TOTAL TINE 

RANGE TOTAL RANGE 

NO. 


KG. 

METERS 


KG. 

PIN. 

MIN. 

KM. 

KM. - 


INITIAL NEIGPT 

4176C5.C 








1 

NU £ TO 

412807.1 

C.O 

0.299 

6801.1 

10.000 

10.000 

0.0 

0*0 

2 

CL/1500 

61C826.7 

657.2 

0.500 

1882.6 

0.805 

10.805 

6.38 

6.38 

3 

CLIPS 

399927.5 

10355.3 

0.900 

10997.3 

6.033 

16.838' 

56.12 

62.71 

6 

CRUISE 

37739C.7 

10765.7 

0.900 

22536.7 

61.519 

56.357 

678.11 

740*81 

5 

CL/ACCEL 

354615.6 

15079.6 

2.320 

20775.2 

9.602 

65.759 

282.65 

1023.46 

6 

CRUISE 

25C890.6 

20266.6 

2.320 

65725.2 

78 . 576 ~ 

144. 333 

3287.02 

4310.46 

7 

OESCENO 

288352.2 

657,2 

0.500 

2538.1 

17.781 

1*2.116 

306.12 

6616.59. 

8 

OES/LANC 

217957.5 

0.0 

0.300 

396.7 

1.645 

1*3.758 

13.86 

4*30.47 

9 

TAXI 

266936.3 

C.O 

0.0 

1023.2 

5.000 

1*8.758 

0.0 

6*30.47 

10 

5* RES. 

28C600.6 

C.0 

0.0 

6533.7 

0.0 

168.758 

0.0 

6*30,47 

it 

CL /1 500 

279357.5 

657.2 

0.500 

1063.1 

0.663 

169.221 

3.61 

4634.07 

12 

CL/ACCEL 

266916.7 

17285.0 

0.950 

12660.8 

10.914 

180.135* 

187.66 

4821.73 

13 

CRUISE 

266557.7 

17291.9 

0.950 

1919.0 

6.675 

186.810 

114.32 

49|**04l 

16 

OESCENC 

263318.3 

3068. 0 

0.356 

1679.4 

11.899 

190.709 

161*04 

5097.09 

19 

LOITER 

267298.6 

3C68.0 

0.333 

16019.7 

30.000 

228.709 

0.0 

5097.09 

1* 

OES/LANC 

264*31.1 

C.O 

0.300 

7*7.5 

3.366 

232.073 

24.*9 

5126.77 


TCTAL FUEL USEO 171077.6 



TABLE 46 


GE21/J10 B1 BASELINE ECONOMIC MISSION SUNWARY, ENGLISH UNITS 


LE6 

OPERATION 

ME IGHT 

ALTITUDE 1 

AACH KC. 

FLEL USEC 

tine total TINE 

RANGE 

TOTAL RANGE 

NO. 


POUNOS 

FEET 


PCUNOS 

FIN. 

NIN. 

N. N. 

N. N. 


~TnI Ti A L~N EIGF T 

92C67C.3 








1 

HU « TO 

9 10064. 1 

O.C 

0.299 

10506.2 

10.000 

10.000 

0.0 

0.0 

2 

CL/ 1500 

9C5934. 1 

1500*0 

0.500 

4149.9 

0.005 

10.005 

3.45 

3.45 

3 

CL INS 

081689. 2 

33974.1 

0.900 

24244.9 

4.C33 

14.030 

30.41 

33.06 

..4 

CRUISE 

•32004.2 

3S32C.6 

0.900 

49605.0 

41.519 

56.357 

366.14 

400.00 

5 

CL/ACCEL 

706202.9 

62597.0 

2.320 

45001.4 

9.402 

65.759 

152.62 

552.61 

6 

CRUISE 

<41303. 6 

66410.6 

2.320 

144099. 2 

70.574 

144.333 

1774.02 

2327.43 

7 

CESCENO 

6357C6.C 

1500.0 

0.500 

5595.6 

17.701 

162.114 

165.29 

2492.72 

• 

OES/LANO 

634037.7 

C.O 

0.300 

070.2 

1.645 

163.750 

7.49 

2500.21 

9 

"TAR I 

C325E2.C 

0.0 

0.0 

2255.7 

5.000 

166.750 

0.0 

2500.21 

to 

9« RES. 

618177.6 

C.O 

0.0 

14404.4 

0.0 

160.750 

0.0 

2500.21 

It 

CL/ 1500 

615070.0 

1500.0 

0.500 

2299.6 

0.463 

169.221 

1.95 

2502.16 

12 

CL/ACCEL 

500450.7 

56709.4 

0.950 

27427.3 

10.914 

100.135 

1C 1.33 

2603.49 

13 

CRUISE 

504219.9 

56732.0 

0.950 

4230.7 

6.675 

106.010 

61.72 

2665.21 

14 

DESCEND 

500517.6 

1C0CC.0 

0.354 

3702.4 

11.099 

190.709 

06.95 

2752.16 

19 

TO ITER 

545200.2 

1C0C0.0 

0.333 

35317.4 

30.000 

220.709 

0.0 

2752.16 

1* 

OES/LANC 

543SC0.1 

0.0 

0.300 

1692.1 

3.364 

232.073 

16.03 

2760.19 


TOTAL FUEL USED* 


377162.2 



TABLE 47. - GE21/ . 11 B3 BASELINE ECONOMIC MISSION SUWARY, INTERNATIONAL UNITS 


M 

S 


LEG 

Nr. 

DMtKATlON 

WL IGHT 
KG. 

ALT ITUOl MACH NO. 
METERS 

FUEL USED 
KG. 

TIME TOTAL TIME 

MIN. MIN. 

RANGE TOTAL RANGE 

KM. KM. 


INllIAL WEIGHT 

M4218.1 








l 

WO r. TO 

5 122*. 7. 2 

c>..» 

0.2 W 

6V 70.8 

ID. ODD 

1 C. 300 

O.C 

0.0 

2 

LL/150G 

510007.9 

457.2 

r .SC* 

223<*. <* 

O.t.22 

10*622 

4.74 

4.74 

3 

CL I Mb 

<***0311.7 

10771.1 

o.w:c 

1318V. 1 

3.312 

13.934 

45 . 08 

50.47 

4 

CRUISE 

**6 8 © 7 2 .7 

11072.3 

0.9 CD 

28146.1 

42 .<*7 3 

56.407 

640.34 

740.81 

s 

CL/ ACCEL 

<•*•7330 .*> 

? j3v j. 2 

2 .3 20 

21342.2 

o .511 

62.918 

198.93 

939.74 

o 

CRUISE 

3c*3?v7 .** 

217-3. C 

2.3 2C 

F4032.6 

FI .066 

144 . 00 © 

3399.06 

4338.80 

7 

DESCEND 

3542 v7.*. 

457.2 

c..5r> 

-000.5 

15.434 

159 . 9*0 

274 . 3 fc 

4618.15 

8 

OES/LANO 

358o76.6 

i . t 

0 .3 0 > 

620. C 

1 .*7C 

16 1.418 

12.40 

4630.55 

9 

taxi 

3*6eb3.7 

J • V 

0 .V 

1742.9 

5.CGC 

166. *16 

C.O 

4630.55 

10 

Vi RES. 

3<*E7o7.0 

:> 

VI . J 

tile. 7 

C .0 

166.418 

0.0 

4630.55 

11 

CL/lf'T 

3<*7*> r .8 .7 

457.2 


1258.4 

0.327 

166.744 

2.46 

4633.00 

12 

CL/ACCE L 

337vG6.? 

1 <* 021.4 

o.vfo 

1C 502 .4 

5.164 

171.929 

83.87 

4716.87 

13 

CRUISE 

3244*,©. 1 

l-7©8.4 

0 .v 5v 

7057.2 

1**.347 

186.276 

245.72 

4962.59 

K 

UESCF NO 

3273 V*.. 1 

3 v 46.0 

0,3*.. 

2555.1 

10.132 

196.408 

133.45 

5096.04 

If- 

LUlTER 

3< 5oft6.s 

3u<*E • j 

0.3 21 

21727.2 

J* v • OC 0 

22 o .406 

0.0 

5096.04 

lo 

LtS/LANl 

3 i'*,**©? .4 

V. . 0 

0.3 0*' 

1173.V 

2 6^0 

22 V . 343 

25.74 

5121.78 




TOTAL FUEL 

use o= 

214724.7 







TABLE 48. - GE21/J11 B3 BASELINE ECONOMIC MISSION SUNMARY, ENGLISH UNITS 


LPG 

NO. 

UPFRATION 

WP IGHt 
POUNDS 

ALTITUDE PACH Nli. 
FEET 

FUEL USED 
POUNDS 

TIME TOTAL TIME 

MIN. MIN. 

RANGE TOTAL RANGE 

N. M. N. M. 


INITIAL WEIGHT 

J 1 440 EC. .0 








1 

WU C TO 

1129312. G 

u. o 

0.2 v» 

15368 • ' 

i: .030 

1C. COO 

O.C 

0.0 

2 

CL/1 SOI 

1124375.0 

1*00. . 


•*7 37.0 

0.47? 

10.o22 

2.SV 

2.5V 

3 

CLIME 

1 

3433d. 1 

C 

24 77 7. C 

3.312 

13.434 

24.64 

27.25 

4 

CRUISE 

lCJiMn.f 

363?©. f. 

\j.9 0 

♦•70S 1 .5 

42.473 

6 6. **07 

372.75 

400.00 

5 

CL/ ACCEL 

98fcl VS 

66bv7.C 

2 .3 2- 

47 >51 .5 

4 .SI 1 

42.418 

107.41 

507.41 

6 

CRUISE 

U0C«?4 .4 

V133S.2 

2 .3 2 V 

n •'240.1 

Ll . -CL 

144.004 

1835.31 

2342.72 

7 

DEStENO 

74/1 IS. < 

ISC'. “ 

C.SC-' 

0014.7 

IS. 434 

159.440 

150.84 

24V3.56 

8 

UES/LAND 

7VC 74© •♦• 

C-.u 

0 .3 ©G 

1342.4 

1 .476 

141.414 

6.69 

2500.25 

9 

taxi 

,'S*.7V4.1 

*J 0 J 

0 .v. 

3452.0 

S. ?\> c 

166.418 

c.o 

2500.25 

10 

St RES. 

7©fct4V.7 

• • - 

'• •*.’ 

17094.3 

c • « 

1 46. 414 

w*w 

2500.25 

11 

CL/ ISO' 

7o6l25.o 

ISC... ■ 

C-.SL 

2'* >4.2 

v.327 

146.744 

1.33 

2501.5b 

12 

CL/AC CP L 

742471.7 

47V7C .5 


23153.8 

S.1E4 

171.V2V 

45.29 

2546.87 

13 

CRUISE 

727-13. 4 

•n^si.i 

. .V*', 

ISSSh.*. 

14 . 3*. 7 

146.276 

132.67 

2679.54 

14 

DESCEND 

7?I VfcC .4 

icor 

C .3 «/> 

*>«»3 ? .v 

K.IS? 

1 4 ©. 4oB 

72.06 

2751.60 

IS 

LOITER 

t.73C?C.l 

1 'j 0 © C* . C 

C.3 21 

47 c *'.0.4 

3v .LOO 

22 6 * 408 

0.0 

2751.60 

14 

DE S/LAND 

♦•71 2**2 . 1 


f ,3T' 

2sea.*' 

< .436 

224.343 

13. VC- 

2765.50 




tctal fuel 

U J F I s 

473367,7 







analyses discussed previously. The calculations said results are suanarized in 
table 49. The nacelle drag increments are those used in the detailed analysis 
and are shown in table 49. Specific fuel consumption increments were obtained 
from plots of installed performance (figures 31 through 34, 37, and 38). 

Takeoff thrust increments were obtained from figures 30 and 36. Propulsion 
weight increments were obtained by adding 198 kg (445 lb) per engine for 
miscellaneous propulsion systems to the nacelle and engine weights of tables 
19 through 22. Included in the table 49 is a column for the "revised" base- 
line. 

Table 49 indicates good agreement of the sensitivity method relative to 
the detailed analysis for engines with small changes relative to the baseline 
(V5CE 502B and VCE 112C). The error increases as the total takeoff gross 
weigh 1 ratio, increases. Reasons for error include the following: 

(1) Reading and extrapolating the sensitivity curves may produce same 
errors. 

(2) Sensitivities for changes in transonic acceleration thrust were not 
included. Tables 24 through 27 indicate a wide range in thrust -drag ratio at 
mach 1.2. Thus, acceleration times and fuel used may vary widely from the 
baseline. 

(3) The assumption that the aircraft will cruise at minimum SFC is 
optimistic; the lift-drag ratio characteristics may tend to drive the operating 
point to a high SFC. 

(4) Scaling factors of the detailed analysis are determined for small 
changes and are therefore of questionable accuracy for large changes in air- 
craft characteristics. 


While the differences in takeoff gross weight using the two methods are 
large for the large aircraft, the accuracy of the sensitivity method. is suffi- 
cient to indicate when changes may be of interest and when they are definitely 
not advantageous. Thus, the sensitivity method allows the user to easily iden- 
tify problem areas. For exanple, '.tan table 49, the largest contributor to 


vehicle weight increase of the VCE 112C is effective takeoff thrust (Rp.^ „ 
is 1.10) while the largest contributor of the GE21/J10 B1 is propul- 


propul- 


sion system weight (R^ is 1.16). 
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TABLE 49. - TAKEOFF GROSS WEIGHT SENSITIVITY CALCULATIONS 


Vehicle/ci gine 

Variable 

Baseline 
VSCE S02B 

Revised 
Baseline 
VSCE S02B 

Refined 
VSCE 502B 

VCE 112C 

■ 

■ 

. Propulsion weight, 
kg (lb) 

31 498 
(69 493) 

30 882 
(68 081) 

33 374 
(73 576) 

35 014 
(77 192) 

38 298 
(84 432) 

33 300 
(73 412) 

SR\ kg/hr, daN (lb/hr/lb) 
Supercruise 
Subcouse 

1.387 (1.36) 
1.005 (0.985) 

1.387 (1.36) 
1.005 (0.985) 

(0.387 (1.36) 
0.981 (0.962) 

1.474 (1.445) 
0.964 (0.945) 

1.469 (1.44) 
1.137 (1.115) 

1.448 (1.42) 
1.086 (1.065) 

Effective takeoff 
thrust, daN (lb) 

27 899 
(62 723) 

26 957 
(60 604) 

28 480 
(64 028) 

23 187 
(52 130) 

23 677 
(S3 231) 

18 713 
(42 070) 

Nacelle drag coefficient 
niach 1.2 
mach 2.32 

0.00057 

0.00047 

0.00040 

0.00030 

0.00061 

0.00033 

0.00085 

0.00048 

0.00202 

0.00095 

0.00276 

0.00115 

R wr 

1.000 

1.000 

1.032 

1.059 

1.16 

1.031 

^SFC 

Supercruise 

Subcruise 

1.000 

1.000 

1.000 

1.000 

1.000 

0.998 

1.071 

0.994 

F W» M ■ 

1.050 

1.008 

Sne TO 

1.000 

1.000 

0.978 

1.100 

1.077 


r cd 

1.000 


1.003 

1.021 

1.083 


R total 

1.000 


1.010 

1.266 

1.464 


TOGW, kg (lb) 
Sensitivities 
Detailed 

Error, 1 

323 048 
(712 188) .. 
323 048 
(712 188) 

316 783 
(698 375) 

316 783 
(698 375) 

320 146 
(70S 568) 
320 228 
(705 790) 
-0.03 

401 092 

(884 258) 

402 62S 
(887 622) 
-0.4 

463 708 
(1 022 213) 
514 450 
(1 134 149) 
-9.9 

S10 136 
(1 124 631) 
629 306 . 

(1 387 359) 
•18.9 

















DISCUSSION OF RESULTS 


Bigine shape, airflow- lapse rate with mach number, thrust- lapse rate with 
■ach number , SC, and noise characteristics have large effects on aircraft take- 
off gross weight. «s an example of engine shape effects, a comparison of cross- 
sectional area variation of nacelles with VSCE S02B and VCE 112C engines is 
shown in figure 3. The only significant difference in shape is that the VCE 112C 
has a smaller nozzle exit area. This results in drag and takeoff gross weight 
changes as shown in table 4. Thus, a nozzle that is 0.043 m (1.7 in.) smaller 
in diameter results in takeoff gross weight increment of 3 200 kg (7 000 lb) 
just due to the nacelle drag change. As shown in figure 25, a one-drag -count 
change results in about a 1 -percent change in takeoff gross weight. 

Engine airflow lapse rate with nach number directly affects inlet capture 
area. For exa^le, the GE21/J10 B1 has approximately 16-percent lower super- 
sonic cruise airflow relative to takeoff airflow than does the VSCE 502B. The 
smaller capture area results in approximately 4 -percent lower inlet recovery 
at static conditions and therefore reduced takeoff thrust. In addition, the 
smaller capture area results in a more rapid increase of nacelle cross-sectional 
area with nacelle length, and therefore higher drag. 

Bigine takeoff thrust and thrust lapse rate with nach mmber has a signifi- 
cant effect on engine size required to meet takeoff distance requirements. For 
exaaple, the VCE 112C has 6-percent lower takeoff thrust at static conditions 
and 20-percent lower thrust at nach 0.3 (at reduced power to meet noise re- 
quirements) than the VSCE 502B for a given static takeoff airflow. This 
resulted in a considerable increase in engine size to meet balanced field 
length requirements 

As indicated in figure 27, a change of 1 percent in SFC at supersonic 
cruise results in a 1- percent change in aircraft takeoff gross weight or about 
3 200 kg (7 000 lb). 

Engine exhaust noise characteristics have a significant impact on aircraft 
takeoff gross weight. All four engines were assumed to employ thrust cutback 
at the takeoff noise measurement point. However, all the engines did not take 
full advantage of the extra growd attenuation while the aircraft was still on 
the ground. For exaaple, the GE21/J11 B5 has 20-percent lower thrust -per-unit 
airflow than the VSCE 502B. Thus, the OiZl/Jll B3 yields approximately 6 db 
lower sideline noise at takeoff on the ground, but it must be sized larger to 
meet the takeoff distance requirement. 


REOOfMENDATIONS 

The sensitivity method has been shown to be a valid method for preliminary 
assessment of propus ion system modifications, and it is therefore recommended 
to be used for this purpose. Continued airframe/propulsion integration studies 
and coordination effort between engine and airframe manufacturers in the afore- 
mentioned high-sensitivity areas are also recommended. 
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AFPBOEX 


ADVANCED SUFBBOMGC TRANSPORT 
NACELLE DSTALLATICN DRAG INGRBOfT 
COMPUTER PROGRAM 




This appendix describes a computer program developed to determine nacelle 
incremental drags for the NASA arrow-wing supersonic transport configuration. 
Program inputs include fr e es trea n Mach readier and parameters defining nacelle 
geometry. Data points are stored internal to the p rogra m to r epr esent wave 
and friction drag values with up to six independent variables. Given a set 
of input parameters, the p r ogra m sets up arrays so as to p e r f orm an internal 
search. Then, a cubic is fit to these data points, and the drag increments 
are determined. 

The computer program requires 40K bytes of me m o ry, and each case, 
consisting of four evaluations of nacelle drag, requires about .1 seconds 
of execution time on a IEM 370/168. 

This appendix relies heavily an references 1 and 2; thus, the reader should 
have a copy of those references as he reads the following. 
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SYMBOL LIST.- The following symbols appear in the coding for Subroutine 
NDRAG: 


AC 
AC 24 
Ad 

AMAX - 

AMAXM - 

AN 

ANM 

0X312 - 
0X3232 - 
CDCM - 
CDF12 - 
CDF232 - 
CDW12 - 
CDW252 - 
DC 
DER 

F2TM2 - 
FTM 
IP 
M 

M2TF2 - 

MTF 

L 

LIPN - 
LM 
PI 
SI 

S(2) - 

S(3) - 

SREF - 
SRErM - 
XMAX - 
XMAXM - 


Capture area, sq. ft. 

Capture area, sq. a. 

Input value of AC saved, sq. a. (sq. ft.) 
taxinx area sq. ft. 

Maximum area sq. ta. 

Nozzle area sq. ft. 

Nozzle area sq. m. 

Nacelle drag increment at Mach 1.2 
Nacelle drag increment at Mach 2.32 
Nacelle drag increment at input Mach lumber 
Friction drag increment at Mach 1.2 
Friction drag increment at Mach 2.32 
Wave drag increment at Mach 1.2 
Wave drag increment at Mach 2.52 
Diameter of capture area, m (ft.) 

Straight line slope for curve defining Mach number 

effect on total drag 

Conversion factor from sq. ft. to sq. m. 

Conversion factor from feet to meters 
Set to 1 for use in equation defining CDCM 
Mach number 

Conversion factor from sq. m. to sq. ft. 

Conversion factor from meter to feet 
Total nacelle length, M. ( f t.) 

Set to 2 for use in equation definig CDCM 
Total nacelle length, m. 

'ir / 

Set to 929 sq. ra. (10000 sq. ft.) for basepoint wing area 
Slope at point where Mach is 1.2 
Slope at point where Mach is 2.32 
Reference wing area, sq. ft. 

Reference wing area, sq. m. 

Length to maximum area, ft. 

Length to maxinun area, rc. 


Refer to the discussion of subroutine NDRAG for definitions of the 
symbols associated with data statements in that subroutine. 
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COMPUTER PROGRAM 


Problem Description 

Nave and friction drag increments were determined for a range of para- 
metric nacelle shapes in reference 2. In order to determine drag increments 
for any nacelle shape of interest, values oust be interpolated and/or extra- 
polated. A major part of the program is involved with the determination of 
the data points to use when calculating the drag associated with a given set 
of values of the independent variables. Once these data points are found, a 
cubic equation (reduces to a linear fit for two points) is fit, and the desired 
dependent variable is calculated. This sorting process and the resulting cubic 
fit occurs in a definite order and is repeated many times before the final 
answer is determined. 

The numeric method for fitting a cubic is similar to Heraite interpolation 
in that the coefficients of the cubic are determined by two points and the 
derivatives of this cubic also satisfy the slopes evaluated at these two points. 
The derivatives are defined by passing parabolas through those points using 
a total of four points for the derivative evaluation. (A cubic is also de- 
fined using a total of three points with a modified definition for the two 
slopes). These points are chosen so that the desired independent variable 
is in the middle interval. 


Method of Solution 

There are four possible situations that must be handled by the program: 

1. Given four points with the value of the independent 
variable in the middle interval, 

2. Wien the desired value is in the last or first interval, 

3. Wien there are only two values, and 

4. When the value is outside the range of data. 

Given four points with the desired value of the given independent 
variable in the middle interval, as shown in figure A-l, assime the equation 

Z = AX 3 ♦ BX 2 + CX * D (1) 




( 2 ) 


will go through the points and P3. Also, the derivative 

V * 3AX 2 ♦ 2BX ♦ C » A 1 X 2 ♦ B X X ♦ C 1 

is evaluated at P 2 using the points P^, P 2 and P- for the calculation of the 
coefficients A 1 , B 1 and C 1 . It is also evaluated at P 3 using the points P2, 

P3 and P4 so that equation (1) passes through points P2 and P3 and the 
derivation equation (2) , also satisfies the slope at these two points. There- 
fore, four equations result to determine A, B, C, and D, in equation (1). 

This is the basic numeric method used by the program to determine the de- 
pendent variable Z given Xq. The detailed equations are contained in a 
following section. 

Wien the desired value is in the last interval, the cubic is defined 
by points P3 and P4, and the derivative at P 3 is determined by passing a 
quadratic through P 2 , P3, and P 4 while the derivative at P4 is defined by 
passing a straight line through points P3 and P 4 . Similarly, when is in 
the first interval, the cubic is defined at points Pj and P2. The slope at 
P2 is derived by fitting a quadratic through Pj, P2 and Pj. The slope at 
Pi is the slope of the straight line passing through points Pi and P2. In 
both cases, only three points are used. In the event that only two points 
are defined, then linear interpolation is used. When constraints are violated, 
the process determines points as if the value of Xp is in the first or last 
interval. Then, a linear or cubic curve fit results in the calculation of 
an extrapolated value for S. 


Program Description 

OPERATING ENVIRONMENT.- This program was written using standard Fortran 
statements using a IBM SYS370 model 168 conputer with the operating system 
0S/VS2. 

PROGRAM SPECIFICATIONS. - Source listings of the main program, NDRAG, 
and subroutine NDTLAE are at the end of the appendix; memory requirements are: 
NDRAG - 9874 decimal bytes 

NDTLAE - 3806 decimal bytes 

The program uses 40K bytes including system subroutines, and no common is 
used. 


PROGRAM MODULE DEFINITIONS. - The main program calls subroutine NDTLAE 
to initiate the interpolation for the total drag given values for the inde- 
pendent variables. Four calls are made in order to calculate the necessary 
quantities at each of two mach numbers. 
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Main Program Module (NDRAG) : The main program performs the followiqg 
functions: 

(1) Setup of data for subroutine NDTLAE, 

(2) Conversion of units, 

(3) Calculation of the effect of Mach number on drag increments, 

( 4 ) Input and output and checks on independent parameter values, 

(5) Setup of data in data statements. 


A large portion of this module is involved in setting up the input data 
in the correct format for subroutine NDTLAE. The input variables required 
for this subroutine are 


of units. 


*MAX , \ , ^SttX 

L Ac Ac 


L , Ap and M. Here, 

dc 




..JTjlEc 


(Refer to figure A- 2). 


and these are set up and used in the English system 


The main program also includes the conversion of data from English to 
International units (and vice-versa) for input -output convenience. 


Hie calculations for the Mach number influence on total drag increment 
start at statement lumber 83. HER, the slope, is determined from the 
equation 


DER = 


CD0232 - CD012 

2.32-1.2 


The slope at Mach i.2 then has a value of 2.0 * DER, and the slope at 
Mach 2.32 is 0.35*DER. (These values for the end derivatives are estimated 
from a study of the basic shapes or trends inherent in curves such as 
Figure 61 on paje 87 of reference 2). Using these two points and the two 
derivatives, a cubic is fit, and the total drag increment at any Mach number 
is determined. 

Additional functions include the input and output of data and the 
necessary calculations to determine if a constraint has been violated. If 
this is the case, the value of the variable is printed out and the program 
extrapolates. 

The ranges are: 

Indep. var. #1 .4 * *MAX - 1. 


L 
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Indep. var. 12 


1.0 ± \ 2.0 

Indep. var. #3 1.25 - ^VOC * 2.0 

AC 

Indep. var. #4 5.5 -JL s 7.0 

“c 

Indep. var. #5 20. - Aq - 40. 

Indep. var. #6 1.2 - M - 2.32 

Independent variables 3, 4, 5, and 6 are required for the calculation of 
friction drag. 

Included in the first part of the module are several data statements. 
These contain values of the independent and dependent variables as determined 
from Figure 24 through 59 for wave drag and from Table XIII for friction drag 
of reference 2. Each array has a special definition and the order of the 
variables is specified by the interpolation method from subroutine NDTLAE. 

Values of the independent variables are: 

todep. var. fi. W - .4, .5, .6, .7. .75, .8, .85, .90, 95, 1.0 

Indep. var. #2, \ = 1.0, 1.25, 1.50, 2.0 

Indep. var. 13, ^MAX = 1.25, 1.50, 2.0 

Indep. var. #4, L = 5.5, 7.0 
d C 

Indep. var. #5, Aq =20., 30., 40. 

Indep. var. #6, M = 1.2 and 2.32 

Each data statement array is defined below. 

NX is the array of the number of values for each independent variable 
in the order as specified above. 

X is the array of individual values of the independent variables in the 
order as specified above. 
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Z1 is the array of values of the dependent variables for friction, drag 
from reference 2, Table XIII. 

The order of these variables in the Z1 array will be described by using 
notation defining the independent variables used for the calculation of 
friction drag. Thus, variable X 3 has three values X,(l) , X 3 (2) and X 3 (3); 
variable X 4 has two values X 4 (l) and X 4 (2) ; variable X5CI), X$ (2) and Xij(3); 
and the last variable Xg has two values X$(l) and X^(2) . The first variable 
in the Z1 array corresponds to the following order for the independent vari- 
ables: 

for Zl(l) use X 3 (l), X 4 (l), X 5 (l), ^(1), 

Then, for the second element , 

21 (2) use X 3 (2), X 4 (l), X 5 (l), XgCl), 

21(3) use X 3 (3), X 4 (l), X 5 (l), X^l). 

Then, for Zl(4), Zl(5), and Zl(6) repeat the first three lines with X 4 (l) 
replaced by X4(2). Now a total of six Z1 values are defined for Zl(l) 
through Zl(6). Values for Zl(7) through Zl(12) are obtained when these first 
six lines are repeated with X 3 (l) replaced by X$(2) . In the first 6 lines 
replace X 3 (l) by X5(3) then Zl(13) through Zl(18) are determined. So far, 

18 lines or values for Z1 have been defined using Xg(l). If these 18 lines 
are repeated with X^(l) replaced by X^(2) then Zl(19) through Zl(36) are 
specified. These 36 lines represent the Z1 matrix for the friction drag. 
This is the definition of the data sequence required for the Z matrix as 
used in subroutine NDTLAE. 

ARW1 through FKW1 are arrays of values of the dependent variable for 
wave drag, for Mach 1.2. The basic method for ordering the dependent vari- 
ables is the same as for the friction <* 0 rr^y Zl, with two additional 
independent variables. 

Arrays ARW2 through FRW2 are values of the dependent variables for wave 
drag for Mach 2.32. Again, the method for ordering is the same as for 
friction drag only with two additional independent variables. 


Subroutine NDTLAE: The main function of this subroutine is to determine 
the value of the dependent variable Z for given values of the independent 
variables X(I). These points are contained in the DATA statements appearing 
in the main program. 
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Hie calling sequence for this subroutine is 

CALL SUBROUTINE (NIV, Xq, X, Z, NX, RES, A, B, C, D) 

where 

NIV Number of independent variables (4 or 6) 

Xq Array of given values of independent variables 
X Array of independent variables 

Z Array of dependent variables 

NX Array of number of values given for each independent 
variable 

RES Final required Z for given XG using X and Z arrays 

A,B,C,D Coefficients of the cubic fit to the data points thus 

defining Z values (In some cases this reduces to a linear 
fit with two coefficients C and D) . 

Initially, various arrays are defined for the purpose of determining 
the location of the first Z to be used from the Z array. The object is to 
determine which interval of input points in the X array will be used to 
calculate the Z array given Xq. After these subscripts and indicators have 
been set, the evaluation the necessary equations for the Z array can be 

made. The cciposition of this array, and the sequence of steps leading to 

the final Z value will be outlined here. (Refer to the section 'Method of 

Solution" for a description of the number of intervals used given the X array 
and to the section "Equations For Curve Fitting" for the equations used for 
the curve fit.) 

The procedure, for calculating the location of ;he first Z (LISTZ) to 
be used, starts with defining the ICF array where 

ICF(I) = Number of values of each independent variable (I) to be 
used for curve fit (i.e. 2, 3, or 4) 

Next, the IGAT(I) array is set to the following values, if 

IGAT(I) = 1, X Q (I) is in one of the middle intervals 

= 2, Xq(I) is in the first interval 

* 3, X G (I) is in the last interval. 
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At the same time, the L0C (I) array is defined as the location in the X array 
of the first value to be used for independent variable I. Thus, the subscripts 
for the X array for the polynomial curve fit are obtained. 

The location of Z( I] given values of X(I) array are then determined. The 
equations to do this are contained between statement numbers 32 to 35 in the 
listings at the end of the appendix. The formula, 

NIV-1 

Z(I) index = 10C(1) + y^NX(l)*NX(2). . .*NX(k)*(L0C (k+l)-l] (3) 

specifies the location (or index) of the element number in the Z(I) array 
associated with each sequence of values for the independent variables X(I) . 

(That is, an analytical gearing between the X and Z array is defined.) 

Statement number 50, which sets LL, through statement number 2000, where 
the polynominal for the curve fit is evaluated, is described by a step by step 
analyses of the coding in the section "Step By Step Analysis of Subroutine NETLE." 
That section also describes, by specific example , the evaluation of the LLCTR 
sequence at statement number 410. Thus , the location of the next Z to use is 
developed. Figure A-3 describes the flow within this subroutine with enphasis on 
control during the evaluation of the Z matrix. (Refer to listing of subrouting 
NDTLAF for complete symbol definition.) 

An exanple of the use of equation 3 will be demonstrated with the friction 
drag data. For convenience, number the variables 1 through 4 (instead of 3 
through 6) as described in the previous discussion of independent variables; 
therefore, 


Variable 

n 

is 

MAX with NX, 

■*r 1 

= 3, 

Variable 

#2 

is 

L with NX2 

d c 

= 2, 

Variable 

#3 

is 

Aj, with NX3 

• 3, 

Variable 

#4 

is 

M with N X 4 

= 2. 


Let subscript refer to variable number: 


Z(l) is 

dependent variable for Xj(l) 

x 2 aj 

X 3 (l) 

x 4 d) 

2(2) 

Xj_(2) 

x 2 (l) 

x 3 (l) 

x 4 d) 

2(7) 

Xi(l) 

x 2 (l) 

X 3 (2) 

x 4 (i) 

2(21) 

Xl(3) 

X2(l) 

X 3 (l) 

X4(2) 
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Fran equation (3) for Z(7) and for element 1121, 


UJC(l) ♦ 3 * [L0C(2)-1] + 3*2 [I0C(3)-1] ♦ 3*2*3[I0C(4) - 1] 

Then the eleaent is: 

1 + 0 ♦ 6 ♦ 0 » 7 

For Z(21) t values in equation (1) for eleaent 3112 give the eleaent 
3 + 0*0 + 18 - 21 

Because of the number of variables in each array, the value of equation (3) 
would always be 1, and new Z values are specified by LL arrays vdien modified 
by LLCTR values. That is, the L0C appearing in this exaqsle is aodified by 
LL array. 

PROGRAM LOGIC.- Figure A- 4 shows a prograa flow during the calculation of 
the Z array. The last Z calculated is the final result. Note that in this 
figure the required niraber of Z values is determined by the number of values 
contained in the next array of independent variables. Figure A- 5 shows in more 
detail the series of Z as calculated for the friction drag. ZQGG1 follows as 
a result of ten previous curve fits where three are linear and seven are cubic. 

INPUT-OUTPUT SPECIFICATION.- Standard input -output functions are utilized 
to input case data and to output the calculated total drag for each mach nuaber. 
The input data is printed out for identification purposes. 

Input: A read tape S is used to read 1 case of 1 card containing seven 

numbers with a E10.6 format. In order. 


Card Column Identification 


Word #1 - A c 

Capture area 

(CC2-10) 

Word #2 - A,y^ 

Maxinun area 

(CC12-20) 

Word #3 - Ajy 

Nozzle area 

(CC22-30) 

Word £4 - 

Length to maxinun area 

((£32-40) 

Word - L 

Total length 

(CC42-50) 

Word 16 - SREF 

Reference wing area 

(CC52-60) 

lford #7 - M 

Mach number 

(CC62-70) 
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Figure A- 4 


PROGRAM FLOW FOR Z ARRAY 
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Normal input is asstmed to be in the English units with the basic length 
measured in feet. A minus sign in colimn 1 indicates the units are in the 
International units with the basic length measured in meters. 

Output: A write tape 6 with an F format is used for printing data. The 
above input data and the friction drag, wave drag and total drag for Mach 1.2 
and 2.32 are printed. Also printed is the total drag for an input Mach number. 

RESTRICTIONS. - Each figure for wave drag in reference 2 contains values 
of Aj^/Aj. for a given value of A^^/A^. Physically, A^/A^ cannot exceed 

Amax/Ac* ^ >ut a ■athematical sense a numerical value for the drag must be 

determined in order to sustain the interpolation process. In this case, when 
a limit is exceeded, the extrapolation assumes that the limiting values of 
drag are used for values of A^/A^ greater than A^^/A^.. That is, on figure 24 

(reference 2) with two curves for A^/A^ of 1.0 and 1.25, drag values for the 

1.50 and 2.0 curves are assumed to be identical to ♦'hose for the 1.25 curve. 

The number of points used in the curve fit depends on where the desired value 
of A^/A^ falls relative to the index values of 1.0, 1.25, 1.50 and 2.0. If 

1.5 cA^j/Aj. S2.0 then last three points would be used. Wien 1.25<A^A^<1.S0, 

a total of four points would be used for the curve fit since this is the 
middle interval. When extrapolation for the other independent variables is 
required, no limiting values of the variables are substituted. 

The terminal point for each wave drag value is specified by the condition 
that the boat tail angle does not exceed 10°. In order to satisfy this require- 
ment, the final drag value for each curve (i.e. each Ajj/A £ value) is repeated 
in .05 increments for Xj^y/L until X?ftx/L is 1.0. Thus, the proper drag level 
as set by the maxi nun boattail angle will be meet. 

DIAGNOSTICS. - When a desired value for an independent variable does not 
fall within the given data range, the variable and its value will be printed. 
The calculations will use this value for the extrapolation as described above. 

Sub rout ine NDTLAE contains several tests to determine if calculated index 
values are set properly. The statement number of the test generating the 
diagnostic is printed out. 

TEST CASES.- Cases 1 and 2 (tables A-I and A-II are presented herein and 
the values of CDO have been verified by cenparing with values read directly 
from table XIII for friction drag (CDF) and figures 24 through 59 for wave 
drag (CDK) in reference 2. Cases 3 and 4 are presented for the NASA LTV 
nacelle, described in reference 3. These cases were checked by cross plot- 
ting the curve data. The parametric study is based on a nacelle that had 
straight line radius connections of the inlet face, maximum cross section 
and the nozzle exit. The NASA LTV nacelle deviates from these straight line 
connections and, therefore, it is necessary to simulate the nacelle with 
straight lines. Case 3 (table A-III) simulates the LTV nacelle using the inlet 
capture area, actual maximum area, and the nozzle exit area. A considerable 
difference in CDO may be noted. It has beat determined that foe best simu- 
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lAtim is obtained by using the na ri a w cross-sectional area that occurs at 
tne intersection of straight lines ori ginating fn» the inle t «i the n ozzle 
and those slopes nearly natch the slopes of the actual nacelle as illustrated 
in figure 5. The irina area deviation will be aepitudnately 101. Case 4 
(table A- IV) represents the later siwulation and the CIO’s agree with the cal- 
culated CDO’s. The difference in QM is attributed to the fail-ring consid- 
erations: 


1. The wing was resized from 10,996 sq.ft, for the reference con- 
figuration to 9819 sq.ft, for the paranetric study. 

2. the narH fa* were relocated further inboard and forward in 
accordance with the ground rules of reference 1. (See figure 
21 of reference 1). 

The above discussion concerns only the 1.2 and 2.32 Mach maber data points. The 
additional CTO at any specified Math nuaber is derived frtm a cubic curve fit 
described on page S. The deviation of CTO nay be as auch as -.0002 between 1.3 
Mach and 1.5 Mach diK to the fact that the true drag coefficient as a function 
of Mach nuaber does not necessarily follow a cubic. 
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TABLE A- 1. -Case 1 Output 

NACELLE CkChfTRY 

CAP HIRE AREA 

MAX 1*11*1 AREA 
NOZZLE AREA 

LOC. OF MAX* AREA 

TOTAL length 

HIN6~REFFXENCE AREA 929.T)3 — 5GN 1 lOOtJO.TJO — SO FT I 

INCREMENTAL NACELLE 1»PA6 COEFFICIENTS 

MACH 1.2 COF* ' 0.00 C*5 3 COM* O.OGG6G ' COP* 0.000*3 

MACH 2.32 CPF* C.OOCK7 CPW* O.GCOOT CLC* C.C0054 

M ACH 1*20 COO* O* 0 0093 

XMAX/L* 0*600 AN/AC* 1.250 AMAX/AC* 1*500 L/DC- 5*500 


~2 .79 — SC N ( 39*00 SO FT! 

4*10 SQM ( 45.00 SQ FT) 

3.40 SON I 37.50 SC FT) 

6.22 ~ M t 20*40 FT ) 

10*36 M ( 33.90 FT ) 


INCREMENTAL NACELLE I* AG COEFFICIENTS - LTV REFERENCE VEHICLE 

MACH 1.2 COF* 0*00072 COW* -0.00041 CDw* 0.0 >031 

MaCN“ 2*32 COF* - 0.00 OoO OJW^“^W>OIX CDD= — U.'iOO^Z 

MACH 1*20 COO* 0.00032 


TABLE A-II.-Case 2 Outp t 

NAC Lilt <H*kTkV 


CA*» Tl Pt AktA 

3.11 

>L M 

( 

sc 

FT) 

MAXIKI.H Akf^ 

1 

SC N 

< 

K-.l-l SC 

FT) 

NL*2/LL AM A 

\ f 

S< M 

c 

v “ • C *! SU 

FT) 

LlC. CF MAX. AM- A 


M 

< 

•♦2 .<»o 

FT ) 

TLTAl LfcPKf.TH 

1 !) Ai 3 

M 

1 

•.V .9*. 

FT ) 

WING REFFRtNCt A**t A 

*:>. s 

SC M 

1 1 

: r .>o.o; st 

FT) 

JNCFtMENlAL NAC C IL» TP AG 

COIF - 1C IT NTs 




MACH 1.? CI»F = 

:. . .■ ~2 

Cl*W= 


.c7 t:<« = 

v. 0CC99 

MACH ',.3/ COM= 

r.r ric.-l 

aw= 

CoO" 

g.>s a»c= 

to 

Ma<H « .32 




CL-P* 

f . : octo 

XMAX/L= v.*V A!I/AC- 

1 .T>v’’’ 

A MAX/ AC = 

2. ..>30 L/L>C= 7 

INCPtHENTML NACfLLt URAI 

CILF* 1C 

It NTS 

- LTV 

HtH kfcMT.t VLhIClL 

MACH 1.2 CDF- 

vi >2 

CLW= 

— C-.COO—l CPo- 

0. v>0031 

MACI-. 2.32 Cf>f = 

v. >uO t -: 

Cf'W= 

-0.3 C 

Olt CDv= 

;.:cv4 2 

MACH i .3? 

SElS 



Cl L- 

0. ~GC«*2 
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TABLE A- II I. -Case 3 Output 

NACELLE GEOMETRY 


CAPTU»f A°FA 

2.57 

SI- M 

I 27.69 SO 

FT ) 

MAXIMUM AREA 


SC M 

1 51.09 SO 

FT 1 

N027LF ARFA 

A.2« 

SC M 

1 *.6.13 SO 

FT) 

LUC. OF MAX. APcA 

g.‘o 

M 

1 31 .A£ 

FT ) 

TOTAL LENGTH 

12. bl 

M 

1 A-2 . 03 | 

FT ) 

WIHC- PEF r PENCE AREA 

929.5 3 

$f M 

( lOCfO.tf SC 

FT) 

I Nf.F EMFNT AL NACELLE 7*? AG 

CCEFF 1C It NTS 



MACH 1.2 CHF= 

'>.00'fcA 

CI»W= 

—O.vO-GS-* CPC- 

0. j GO 10 

MACH 2.X2 crF = 

o . : o <;>■ t. 

rr *= 

-*.0CC?1 CO G = 

“• I OC *5 

hACH l.e r - 



C 0 = 

: . ' Ov/26 

xre y/L= an/ac= 

1 mt C»C- 

A HA X/AC- ? . S7 7 L /DC = 7 

I».CCE.r :*ntal nacelle erae 

COEFFICIENTS 

- LTV REFERENCE VFhTCLF 

MAC- 1.7 C r » c = 


COW= 

co'= 

o. c: -i 

uae *' : ? cr»E = 

e 

cnw= 

-c.coci* cm* 

o. 

MAO' 1.°' 



rr r - 

r /f-f *.3 


TABLE A- IV. -Case 4 Output 

LL r »> * ' ‘ 1 T v V 


c A r* 1 1 U ; ; ■> r ,i 

2.-7 

MAXIH'W A r A 

S.17 

N077LE A«EA 

9 • 7 9 

LHC. r- f-iiX. tErf 

°.'E 

TOTAL LI NOTH 

l? .r l 

WTN1 ° r E r r ' E *'C C A c t A 

«*29.~? 

I»C.Rtf-FNTAL f ACFLLt. .-SAG 

CUE EE II 

MACH 1.7 CE«e= 


"ACM '' . X 7 C r >E= •; 

7 

PACr- 1.0 


X‘1aXA.= ’..T*'’ ah/ac. - 

1 .01X1 


S< P. 

( 

77. 0*- 

M FT) 

S l M 

1 


SO FT) 

St P 

( 

. 1 * 

Sf FT) 

M 

( 

X"* .<»%> 

h T ) 

J» 

( 

*?.c? 

FT 1 

Sf K 

( 

KCC .'' r 

Si- FT) 

I^-FTS 




ci k= 

** • 

"CX.. C(..' 

= 

cow= 

-c. 

o'joIa me 

= ■;. ^-3 



CEO 

- *. '<>?{< 


amax/ac= :. r i.. 1/^0= f 


1'OEMtNTAL NACFLLt I tW COEFFICIENTS - LTV EE Ft PENCE VEHICLE 


maCh 

1 .2 

cef= 

-.‘.no. 72 

C n W= -v.ci>C*»1 

C0'* = 

r, 000 X 1 

MACH 

7.3 7 

tr»F = 

" « Gv *.<■ ) 

cr w= -c.crsit 

cr : = 

c. :c:«.2 

MACH 

l.oG 




C nr ‘= 

o. vo:a6 


C7** 


07« 
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EQUATIONS FOR CURVE FITTING 

1. Slope using parabolic fit 

General equation of parabola 

Z * aX 2 ♦ bX ♦ c (1) 

Slope of parabola 

Z* - 2aX ♦ b (2) 

If given three independent variable values with corresponding 
dependent variable values then: 

*1 ,Z 1 X 2 ,Z 2 *3» Z 3 

Z 1 * + bXj ♦ c (3) 

Z£ = aX 2 2 + bX£ + c 
Z 3 * aX 3 2 ♦ bX 3 + c 

Solving equation (3) for a and b 

(z r z 2 ) (x r x 3 ) - (z r z 3 ) (x r x 2 ) 

(x^-x^) (x r x 3 ) - (x^-x 3 z ) (x 1 -x 2 ) 

b ( z r z 3 > - ^i; x 3 2 ) cyz 2 ) 

(X 1 z -x 2 z ) (x r x 3 ) - (X x 2 -x 3 2 ) (X r x 2 ) 

Then substituting in (2) above the slope at is 


S = 2&Q 


(Z r Z 2 {(X,) 2 - (X 2 ) 2 } a) 


Xj-X 2 


2. Derivation of coefficients for the cubic 
General font of cubic 
Z = aX 3 ♦ bX 2 + cX + d 
Slope of cubic 
Z'= 3aX 2 ♦ ZbX + c 


To determine the cubic passing through two given points, (X 2 ,Z 2 ) and 
( x 3 jZ 3 ), and having given slopes, S 2 and S 3 , at these points, the following 
system of equations is solved for a, b, c, and d. 

Z 2 - aX 2 3 ♦ bX 2 2 ♦ cX 2 + d (4) 

Z 3 = aX 3 3 + bX 3 2 + cX 3 ♦ d 


S 2 - 3aX 2 2 + 2bX 2 + c 
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S 3 - ZaX^ 2 * 2bX 3 ♦ c 

Solving for the four unknowns: 

- { 2(Z 3 -Z 2 ) - (S 2 -S 3 ) (Xj-Xj) } 
a = — 

(*j - v 3 

. h-h ( i 

b ' rr ■ s 2 ‘ ( v 2 X z ) (x 3 x 2>*»1 ! ^ 

V * 2 

c = S 2 -2X 2 b - a 
d = Z 2 - aX 2 3 - bX 2 2 - 

3. Final calculation of dependent variable 

The dependent variable Z is then calculated using the coefficients 
a, b, c, d and the given value of in the equation 

Z = aXg 3 + bX^ 2 + cX^ ♦ d or 

Z = d + XpCc + XjjCb + Xga)) 
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STEP BY STEP ANALYSIS OF SUBROUTINE NDTLAE 


The sequence of steps to calculate the friction drag will be defined 
by reference to the program steps contained in subroutine NETLAE. Primary 
emphasis of the following discussion is on setting values of counters and 
subscripts leading to the determination of the final interpolated value for 
Z. 


This, refering to the section on subroutine NBTLAE for definition of 
independent variables, 

NX(1) = 3, NX(2) = 2, NX (3) = 3, NX(4) = 2. 

Fra’ statements starting at statement number 3, 

ICF(l) = 3, ICF(2) = 2, ICF(3) = 3, ICF (4) = 2. 

Progr a m steps ending at statement number 23 define the L0C array as: 
L0C(1) = 1, L0C(2) = 1, L0C(3) = 1, L0C(4) = 1. 

Also, for IGAT assume IGAT(l) * 3, IGAT(2) is linear, IGAT(3) =2, 
IGAT(4) - 3. Refer to equation (3) of this appendix for the location of jd, 
first Z to be used, i.e. 

LISTZ = L0C(1) = 1 

The basic logic for performing the calculations, starts at statement 
number 50. 


NIV = 4 

L = 0, LL = 0, LLCTR (I) = 0, I = 2, 3, 4, 5 
50 LL = LIST Z = 1 
K = 1 

L = L+l = 1 
IC = ICF(l) = 3 

For J = 1, 2, 3, and M = 1, 2, 3: 
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ZPR(l) = Z(l) - Z ml , ZPR(2) - Z 2m , ZPR(3) - Z(3) - Z 31u 

The derivative of X (XPR) with vaules from first X array are then 
determined: 

XPR(l) = X^l), XPR(2) - X 1 C2). XPR(3) - X x (3) and IIX * NX X = 3. 

Then, IL * IGAT(l) * 3 for XG X in last internal for variable XI; then 
go to 12. 

At statement #72, 

IPN = 1, IIN = L = 1, IS = 2, M = 3, IRX = 1 and go to 80 

At statement number 80, 

IN = IIN+IRX =2, IP = IPN+IRX = 2 

Calculate slopes S(2), S(3) and A, B, C, D and final ZPR(l) = AX 3 + 
BX^+CX+D for X = XG X (The equations for these quantities are defined in 
Appendix A) 

Note that LLCTR(2) = 0 here 
K = K+l = 2 

LLCTR(2) = LLCTR(2) +1=1 
Test LLCTR(2) - ICF(2) = l-2< 0 and 
K-2 = 0, go to 310 for next Z to be used from Z array 
Find LL = LL+LLCTR(2) *LP where LP = LP*MX(1) = 3 
Note that LLCTR(3) = 0, then 
LL = 3 
go to 50 

Therefore, here Zq XXX in ZPR(l) has been calculated 
and since LL = LL+LISTZ = L0C(1) +3=4, the next 
cycle will calculate Zq 211* (Note LLCTR(2) = 1) 

Reset, K = 1, L = L+l = 2, IC = ICF(l) = 3 

For ZPR, J=l, 2, 3, M=2, 3, 4 and 
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ZPR(2) - 2(4), ZPR(3)«Z(S), ZPR(4)«Z(6) 

Again, for X array, 1IX=3 and since k=l , IX=L0C (1)=1 so that 

XPR(l) = X^l), XPR(2) = X 1 (2), XPR(3) = X x (3). 

Again, IGAT(1)=3, go to 72 and calculate slopes. A, B, C, and D 
for fit so ZPR(2) is defined. 

k=k+l=2 

LLCTR(2) = LLCTR(2) +1=2 
LLCTR(2) = ICF(2) , ICF(2) = 2 
go to 301, K^-NTV, go to 500 
L=L + 1 - ICF(2) = 1 
go to 55 

Note that here Z^-^ and Zq211 in ZPR(l) and ZPR(2) have been calculated 
and next is a linear interpolation with respect to X 2 values. 

Subscripts for Z are defined in Figure 6. 

At 55, IIX = NX(2) = 2 and since K > 1, go to 49 where 
I=K-1=1 and IX=NX(1) +L0C(2)=3+1=4 
IXP =IX + J-1 = 4,5 for J = 1, 2 so that 
XPR(l) = X 4 (l), XPR(2) = X 4 (2) 

IIX=2 so go to statement number 64 for linear curve fit 
for C § D and go to 2000 for ZPR(l) or Zqqj calculation. 

K=K+1=3, LLCTR(3) = LLCTR(3)+1+1, K>2, to to 315 and 
Set KK=K-l+2, LLCTR(2)=0 

Find next Z to be used from array, redefine LL where NX 1 =3, NX 2 =2 
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z 2 1 1 2 ~J (Replace variable 4 by 2, l.e. second value of last Independent variable) 

23112 “ 1 


Note: Subscripts refer to values of independent variable and G refers to 

the desired value of each Independent variable. Thus Hq2 1 1 m®* 0 * 
the H for Xg for variable U\ , the second value for variable £2, the 
first value for variable #3 and the first value for variable #4. 


Figure A-6. 1 CALCULATION FOR FRICTION DRAG 



LL=NX 1 *NX 2 «6 
go to 50 

The next 6 elements of the Z matrix will be processed for Z™ 2 i snd 
Z G221 the same logic as for the first six elements, but now LL=LL+ 
LISTZ * 7 for Z the prime values. The X values are the same as for the 
first six lines. 

The next series of calculations ends with the calculation of Zgg21 
from Zqi 21 and Zq 2 21* Similarly, Zqjji and Zq 2 31 are used to determine 
Z GG31. Since three values of Z for variable number 3 have been determined 
a cubic can be fit and the resulting answer is Zqqgi. Thus one value of 
Z has been determined corresponding to the first value for variable number 
4. 


Next, replace variable 4 by its second value and proceed through t) - 
calculations the same as above. This determines ZqGG2 and one more line 
interpolation produces Zqqgg the desired value of the dependent variable 
given values of the independent variables. 
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SOURCE DECK LISTINGS FOR NDRAG AND fOILAE 
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( JAN 76 » 


IS /.o'.. KiMMtJ *1 EXItNDiK 


DATE ?o.l 62/15.1.4.43 


R£ OUES 160 OPTIONS: 
3RTIONS IN EFFECT* 


ISN 0002 


ISN 0003 

ISN 0004 
ISN 0005 


ISN 0006 
ISN 3*. 07 
ISN 0008 
ISN 0009 
ISN 0010 

ISN 3011 


ISN 0012 



L 1 ST, XPfcF * t.bCK 

NAME I MAIN I UWTlHlZfMI L INf Cl'OuT ( I SIZE (MAX) AUTwth l ( WON fc » 

SOURCE EbCDlf- LIST DECK OhovLl MAP NDPJKmAT GOSlKl XKIP ALL MIANSF Nl)Tt.RMlNAL FLAG 1 1) 


XGI6) ,X(2 4) ,ZIN4( ) 

♦ OP W1 1 14 j ) f tRWUl**0> 


,NXttJ t AH W I 4 1 4 C ) ftPWUl4*>) , 

, F<tWl<2.» ,ARW2I 12 >1 ,bl*W2U2XI, 
, FT W2 ( I i * ) 

, N X X 1 4 ) , XR R ( »% 1 i iVHikl ,S(3» 


DIMENSION 
1CRW11 l4j) 

2CRW2( 1 2 V | ,l>kW/< 1 2C I , (,kHk( LiiC) 

DIMENSION XblKI ,XHU-> ,Zl(3o) 

bOUT VALENCE < XG1 1 1 f ,XM 3M , (XlllSMIM) » < NX 11 I > , NX < 3 ) » 

EC U1 VALENCE I Z < 1 1 , AKWl Cl ) ) , ( 2 (1 *»1> , E F MM \) ) , CZI2CII ,CRMMllt 

i(zi 42 i),r>Rwim) , m&oi ) ,erwii n > , u n *. 1 >, prwi ( 1 )) , 

2(Z(72i) ,arw2(1 n , mfcAii »mrk2U)) , i/CNon.LKHzm > » 

3(211081 1 ,CRW?(1H, ( H 1701 l,FRW? (1 1) , ( Z ( 1 V 1 ) , PKW2<1 ) 1 
REAL L,LM,MTF,M2TF2*M,LW»LR»MLT, LTVW? LTV<‘ 

DATA LW /.33C32#,:C046 t ,3C.0b2* .2 CS46 , ,C "43 , .ICC 4j , .0004 1/ 

DATA ML T /l. 2 . 1 ,3, 1 ,4, 1 .t, , l.e, 2, 32 , 2 ,70/ 

DATA NX /in, 4, 3, 2, 3, 2 / 

DATA X /. 4,. 5, .6, .7, .75, .80, .85, .96, .95,1 . 0 , 1.0,1.25,1.5,2.0* 
11,25,1.5 1,2.0, 5 .5,7.0 ,20 ,,30., 40*, 1.2,2.32 / 

DATA Z1 ✓ . 30034, .30036, .30039, .OOC42,.COC44,.;CC48, 

1.00150, .COC 53 , .00057, .00 :61, .000 of , .0007", .00060, 

2 •00*' 69, .00 1 75 , .00080, .00086, .00092 , .00031, .00032 , 

3.13035, .03038, .C034C, .3*04 3, *0045, .C0C47, .CC051, 

4.00053, .00056, .C/OOvl, .C3C6b, «0'Mc*l, .00066, .0 0071, 

5.00175, .00081 / 

DATA ARWl / . 00 03 8 , .0 > )43, .o0046 ,.00066 ,.00049 ,.00049 , 


l .00048 , 

2- . 01 f . 17 

3- .03C0* 

4- .00026 

5- .06019 

6.C010 , 

7 .00143 , 
8.3 >D3C9, 
9-. 00032 
1-. 0001 07 
l-.00v5>3 
2.00385 , 


00047 , .00046 ,.00)46 , ,".00010 ,".O0vl4 , 


-..10019 ,-.30022 
-.00.10 ,-.00014 
-.003266, -.00 330 
— .0D>?2 , — .00624, 
00106 ,.00115 
C0140 , .00 »4V 


,— .00024 , — .01C26 ,0C0 2C>6 00030 , 

,-.00017 ,-.01019, -.0 3022 , -.10 2n ♦ 

.0 PC 08 ,-.Cn01C ,-.0'1*l4 ,-,70017 , 
,—.000 26 , , 0102 86 *— * Of/1 3C , 

,. 00123, .03128 >.00133 , .00137, .00140 , 

, , jO'.i 47 ,.00045 , *Ou04i> , *t>OU377, .1)00345 , 
0121 ,-.010C34,-.C;0lC7,-.*.'001V,-.a0u27 


3 

4 

5 


00 :,?.n , .00121 ,. 

-•0003o2» — .0004i ,-.00*'45 5,-.C.lC5(*3,— .000545,— • 100034 , 

-.00019,-. O0027 ,-.16032 ,-.000362 ,-. 00040,, 030455, 
-.ICC 546, .00318 , .C03IH ,.00^24 , .00J47 ,.00368 ,. 001*85 
03385 ,.00385 , .00365 ,.30216 ,.Oc*2l4 ,.002lo ,.00223 • 


ZFR1-2 
ZFRl-2 
Zfkl-2 
ZFR1-2 

ZFRi-4.* 
ZFRi-2 o 
ARWl 0001 
ARWl 7002 
AKWl 0003 
ARWl 3 005 
AKWl 0005 
ARWl 0006 
ARWl 0007 
AKWl 0006 
,AKW1000V 
ARW10010 
AKWl 0011 
, ARWL0G1? 
ARWl 1013 



Main 


uS/.U-C* HJKlkAN H ‘-XUN'UO 


DAl £ ?o. ic,2/l‘. 04.43 


O' 

K> 


i Jam 7s I 


ISM 0013 



ISM 0014 


3. r if>224 

4,00128 


5.0A34& * 

• on -ia « 

*.o>0w6 » — 

. •*/.>> * 

— • 1/s, wa2 » — » 

* — #0Cta>6 3 , 


6-.0'*7 a 


,-.” 1 C41 

* . C * C 1 3* 

. 0 w •< J * . . 

OO C A * — » w C i t> 6 * 


7-.0.K>ll 

,-.70015 

*— . a -jo i. f 

n':w< 

*— .Ojv71 * 

-•oco2i *-.:*■; >i 2 * 


b— • 0501 4 

y-.OC-i'llb 

,-.:oc24 

r.r/ 24 

,-.C*r.2 7 

*— .0' '/y * — • o c o 2 »• 

* 

9-.0O 27 

*— . 3CC2<l- 

/ 





OAU tPWl / -. 

**GC12 

I'^lA * -. 

cjcie 

CO 2a *»,!,.vil> * 


l-.O'J >27 

,-.“0026 

• in-c rt- 

, -.00 >7.7 

* — . 0 1 0 ? o 

* — * i.< 0017 * — • j v-< 14 

* 

?— • 0701 ft 

*-.0C V4 


,-.0CA27 

,-.01028 

,-.C0'2fl *-.00027 

* 

3-.0-1O2c. 

* .00075 

r.or>75 , 

O. tk* * 

•OCOaB , 

• vvi *.W'-<J2 , 



4.07'25 T 
4-.W-H4 

?-,o:as3 

8- .0-.353 
9.03223 
1.00160 
1 .00A5B 
2.00011 
3-. 030 IS 
A-.f,?2VV 
S. 03019 , 

6- . 0:3 , 

7- . 0-001 3 
P-.0004 » 

9- . 0^03 6 


*2 3t, ,.4*244 , ."fV44 
, • ‘.*01 1ft 


. C 2 a* * . * 3 ✓ 44 , « ’.> C 1 v f • 0 0 1 2 f 

• • l b l U 94 f . ' A> 1 \ * * CO 1 V f.Ofll/ f 


-'0017 * . V C 1 7 . • - jY( * .t*0' 24 ,*:'.»<>21 ♦••Hl-l *-.'OV)l * 

“•00'’ 22 »•*.(' *— .0 y '• b’l *— .3'.32 -32 *— .wOl * 

-.00032 »-.0 '0a5 »-.<;(u 50 f-,0'63 »-,i.t 60 *-,»0u65 * 
-•OOOS-j *-.0001 » 0, C.2? ,-.00032 t-.OOC/AS ,-.Oou3w » 
-•J»'*055 *-.00066 *-.000 53 *-.vJ05' * 

03215 * .3020 5 *.::i9 «.<?U83 *.C-C>l7r. *.l'l69 OloV « 

00l6« *. r >01o ,.yG14* ,. -12* *.031:3 ,.i 04 , .* j >74 * 
nop 5 *.0' 1 Ct * ,<•.'■ "5 *.J3 10 5 t.wOObf * . jOOoO *.»0:Ci0 * 

• OOOuB *-.3002 *-.0Cv 27 ,-. C003 *-,.‘003 *.'--1015 * 

— •<10(43 * — .-10070 * — *.rOC-8? * — »C 04 a * — .v-M *—.■•'•>102 * 

-•^0085 *.0004' *./0056 *.>005 *.«»OOa *.;».*3? *.vi0u2 fc - 
00' 12 *.0' 05 f.00'06 *-.0>Cl3 *-.<0:14 * 0-:?5 * 

Vl-UCSS , — •* Of’A *■'■“* 

— .*.* 0 ' 1 4 * — • *. j 02** 

— *'iG ja3 *— • >0 04 * ( 

— ,OC*;4 »— .0*i0«.*> « 


.O' 06 *-. G1C13 *-.\C-:iV ,-.0025 * 

- •>'■ v0a 5 *— .t-vw45 *— *v>W'i**3 * — * 

*— I'. 3 *-«C0Ca *— .^ 'OAb * 

r — • 013 »— »— .wCt'25 , — .. >03 * 

, -,r Ohs * - . f r^A 3 *-.r,or>H ✓ 



aRWIOOIa 
AMWt 3016 
ANWl cOlft 
AN Ml 001 7 
A 6 Ml (/ -J 1 1 
AkWl 0014 
AN Ml 0021 

ftKMl OC 2 1 
ftkMluO** 
JrRWl 0C23 
bkWl 0024 
6kWl 0025 
bftWlOG2*. 
OK ML 0 02 V 
tJKWlOu2ft 
tt«Wl0024 
ftKMl s.0 JO 
ftKMl Oil 3 l 
HKMl 0032 

th Win A3 A 
bftMl OOJ a 
PHM i 0036 
6NWIO036 
pmWI O03T 
tiKMi 003b 
t'KWl 303.9 
bPMlCiC4'^ 

ON Ml uv4 1 
Ch Ml 0042 
ilKwlvOA.» 
C.kWl 0044 

CM Ml 0045 
CKM10w46 
LW Ml u0a 7 
CKMl 064ft 
CNMl 0049 
, UkMlOCftvi 
CKMl 0041 



< JAN 75 I 


MAIN 


L S /.» u f j 


K K 1 b AN M tXTt NUt- 


UATfc 7to. lo2/15.v4,*3 


1.0074 f 

2- . 07012 

3- .000*? 

4- .C7732 

5- .0*:*i 

6- . 0*041 
7.0003 ,• 
8-. 0001 
4-. 00 040 


.70037 ,.O7 0U ,. 
,-.0001 ,-.001; iv 
.00041 .'l»0w4 

,-.00735 ,-.0 7:*, 
,-.0)017 ,-.00 02* 
,—•0004 ,—.00034 


7 C ’ * 

,-.C' 


,-.0701 , — . “001 o ,-.j001B ,-. 00314, CRWl 0052 


?i ,-.000 32 ,-. "iO ij.it ,-.w)w , CKW10053 

,-• 0034 ,— .COCl ,-.owOl? .(.ji»25 , CRW1005* 

, “.0"' i 'a? t“ .0 'OAl ,“.vCt»4 »— *00034 , CRW10055 

;*w32 ,-*0**35 ,-. .1004 ,-.'CG*2 , CRWlOoS*. 

,.00114 , .H011 ,.00045 ,.tvO tot' , .00051 ,CRWl005? 

• 0C025 ,.*7*36 ,.n0O3 to ,.*7C3e , .00051 ,.0"039 ,.00016 , CRW10C5B 
,-•00024 ,-.0004 ,-.00055 ,-*G005o ,-.000*6 , OKWlOf-59 

/ CKW10C6G 


1SN 0015 


s") o 

TC- 

-0 O 

O 

Si > 

*e. 

<£ 

C-, * > ga 

- r * 

»- 

*- 


1SN 0016 


DATA l)RWl / -.00005 ,-.00124 ,-.0C0* ,-.000o , 

1- . 07774 , -.00*^74 ,-.00074 ,-.00j?4 ,-.0)06 ,-.000*7 

2- . 0704 ,-.0006 ,-,00065 ,-.00074 ,-.00079 ,-.00074 

3- .0006 ,.00356 ,.0033 , .00 302 ,.03 279 ,.0026 5 t mV j 
4.00240 , .03248 ,.00246 ,.30265 ,.00220 ,.00185 
5.3W85 ,.07061 ,.00091 ,.00041 , .3C091 ,.0vl63 
6.00337 f.OCOOto , -.CO 024 ,-.000 36 
7.00055 ,.00003 ,-.00045 ,-.0009 

8— .0)127 ,470012 ,—."00103 , .00015 

9.0031 , 0. ,- .COCOS ,-*00'V>3 , 

1-.0301 ,-.0002 , —.03 C3 ,-.03036 
1-.00028, -.00017 ,-.0001 ,-.0001 

?— .07037 ,-.00736 ,-.0303? ,-.000 28 ,-.07017 ,-.00* 
3— .0702 ,—.0003 ,—.00036 ,—.00037 .00036 , — .00033 

4.0019 ',.0016 , 70013 5 ,.03107 , .00095 ,.00085 ,.0.7 

5.00*95 , .00C95 ,.0:08 ,.000 5 ,*0001 9 ,-*00312 

6- * 00046 ,-.00045 ,-.00036 ,-.000 3e , 0. ,-.0003 

7— .07077 ,-.00085 ,-.0T092 ,-.00046 ,-.00045 ,-.Gs7C 

8 0. ,-.0003 ,-.00058 ,-.0j177 ,-.0 0 0 8 5 , -.OOG*»2 

9- . 03095 ,-.00095 , -*00)V , -.00^76 / 


, *0 

, m2 

,-.00735 , — .00027 
*-.0011 ,-.C0123 , 
,.00032 , .00036 , 
0 • ,0. ,-.uun 

*— ,00037 ,-.00036 
,-.0002 ,-.0003 


—.00365 , 

5 ,-.00024 , 
*0“G74 , 
255 ,.00248 , 
Cl 35 ,*CCU , 
014 ,.00091 , 
,-.C'.*027 , 
-.70128 , 
.00023 , 

17 , O00 1 , 

, — • 000 35 , 

-.00036 , 

1 , — .0001 , 
,-•00026 , 
083 ,.00095 , 
on 26 ,-.ooo4 
,—•00056 , 

*0 ,—.03076 , 
« 


DATA 

fcRWl 

/ .00536 

,.005 33 

f 

•0046* 

, *C04t>9 

, .00470 

, • OC 47 5 

1 .00475 

,*00475 

,*00475 

,*00475 

, *00436 

, *003o5 

,.2 030 7 

,.37252 

2.00233 

,.00221 

, *002 32 

,*002 32 

V 

•0 023? 

, *00232 

, .00313 

,.0023 

3.001 54 

,.CO09C 

, *00062 

,*02.35 

♦ 

•000 22 

, *00u25 

, uO 2 5 

, *00025 

4.G0121 

, *03733 

, -• JC045 

, — •0 01 01 

,— *0>1?1 ,— mJ 

01* ,—#001*5 , 


5— .0*145 ,-.00138 ,-.00114 ,.00061 ,.00075 ,.0)060 
6.00019 ,-.0*0'S ,-.0071*. ,-.07011 ,.000(3 ,.o*cc:i 

7- . 0001 ,-.00023 ,-.00035 ,-*0;G4 ,-.000*5 ,-.U00? 

8- .00045 ,—.90034 ».00" , '4 , -.0001 , -*0uwi23 , -.00035 


, . C 0 03 6 , 
fiCCCC** , 
,— .00049 , 
, — • C .» ) * , 


OK Wit 061 
DRWi 0062 
DRW1 006.3 
OR Wl 0064 
DKW10065 
DR Wl 0 066 
DRWi 0067 
UKW10066 
DRWI 0069 
DRWI 0070 
OR Wl 0071 
DRWI 0072 
OR Wl 0073 
OR Wl 007* 
DRWI 0075 
DRWI 0076 
DRWI 0077 
DR Wl 0078 
DRWI 0079 
DRWlOOBw 


ERW10CS 1 
fcRWl 0082 
OR Wl 008 3 
fcRWl 0064 
6RW10086 
ERW10086** 
fcRWl 0067 
fckWl 0068 
fcRWl 0C89 


o* 
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< JAN 75 ) 


MAIN 


l S/36 j FORTRAN H IXTfNOlO 


( Alt 76. Ic2/15.CA.43 


ISN 0017 
ISN 0018 


ISN 0019 


4-. 00046 ,-.0005 ,-.OOC'44 ,-.0O0Af. ,-.C1'*39 * .GCG04 ,-. CO ) l , 

1-.C7..23 ,-.00035 ,-.0004 ,-.G0u46 ,-.0005 ,-.u0C»*4 ,-.0004.6 , 

1— .00934 ,.0'2u7 ,.C 1°3 ,.0 lo(> , .1.0126 ,.vCC95 *.<X0o4 ,.C-wC>6 , 

2.00392 ,.00092,.: -092,.: Ill , •- C.9 , .0* 25V ,-*00 mI » 

3— .03 33 7 ,-.0005 25 , .' O' 34 ,.30 34 ,-.;'>0l5 ,-.30019 , 

4 — .J '* j 727 ,-.30049 t-.r-vrt.^ ,-.r.( 7«2 , -. 0(095 ,-.ouo3 ,-.oot 9 , 

5— .0006V ,-.0'i015 ,-.GuGl‘* ,.000/7 ,—.00044 ,-.0v>64 ,— .f'COE? , 

6.0(1045 ,— .7 >103 , — •C'CC.4 ,— . • 064 ,.Co6S9 , .(w3t> , .0C31 ».v0453 , 
7.00417 ,.(3391 ,.0.4i’ ,.0L4^ , .i C**!d ,.0342 ».vC46' r * * * Oc** v «* , 
6.00346 ,.0?o ,.0-7216 , .OVl 77 , ,0017t , 

9.0021 # .0 rt 21 ,.0021 / 

DATA FRW1 / .0034 , .0027a , .0 0*C? y.CCm ,.'0*64 ,.*?02? , 
1.00307 , .0002 ,.0*t(3 | (OOt. 3 , »0v.l6 , »OvcF , — *vWwj 5 | — *vf Ov7 , 

2— .C7798 ,-.00125 ,-.CCl40 ,-.00140 ,-.9)134 »-.9ull / 

DATA ARW2 / .0-*114 , .0 M , . .00* it ,.0>::7 , .00016 , .03*751 , 

1.000 3 4* , .00004 , .00 O'jJi, ,. 003036 , -.OX >49 , -. , 

?-. 00008 , -.010091 , -.0001 , -.000103 , -,MM1 , -.Of. MU ♦ 

3— .0.0118 ,—.00012 , -.000.44 ,— .G';C , — • ' M )Ct ,— *i.C’u09l , 

4— . oo:i , -.oom:; 3 , -.: n , -.s-'ftiu ♦ -. uuin , -.uooi 2 , 

4— .000349 , — .0 ) 06fc , — , —.000041 , —.00*1 , — .wC .1:3 , 

6 “. 00011 , -.000112 , -.0>-’ il6 , — • 0 • 0 1 2 , • jG 031 , .u'.OPbt- , 
7.00726 , .01023 ..(.10-216 , . r>M2 , .OOCll! , .00(17 ,.70017 , 

fi.0<»17 , .000105 , .CC0-.7 , .500 732 , C.C , -.0000? , -.00u03 , 

9— ,0)C .4 , -. ^ 1075 , — .00006 , — * 3000 o , — « J 0 OO 6 , — «tH)01 , 

1 — • *C >1 3 , -. 01116 , -.0V.17 , jCCle , -.CCOIL*,-. 70019 , 

1- .0XH97 , -.0002 , -,07.0c , -. O.'Ol , -.00.13 , -.*..016 ,-.00017 

2- . 0031ft , — • 000 1 66 , -.00014 , ~.00:»1V7 , -.( .02 , 

3.00081 4 , .000826 , .CP**&i* f „f, f 0 64 , . 7* 0037 , , .000837, 

4.905837 , .000837 , .000837 , .070615 , .0uC57 , • • v06l , 

6.00044 , .0304 , .0:c-3>. ♦ .::0032 , • * C C 3 2 ,.00032 ,.00 j32 , 

6.00038 • .030314 , .0**025 , .OC517 , .0 7l? fc , .C*C>Jts2 , . uC0t4 , 

70.0 , 0.0 , 3.0 , .0*0045 , -.00005 , -.0 *13 , -.0002 

8— .0702 3 , — »9102to , —.0002 6 , — . (7(3 , —.3" "3/ , 

9 -.00033? ✓ 

DATA KPW2 / .CCC02G , .CCCM8 , .CC001 , -.>2000?, -.OoOvl , 

1- . 0)0012 , -.00.016 , -.:c,02 ♦ -.00026 , -.(i30:* 25 , -.003 C6fc 

2- .0"C0h , -.CI^V .000102 , -.00011 , -.000114, -.C.'-QIU, 

3- .0301 ? , -.'0C1? , -.C.V.12 , -.7fMo£l , -.CfJOOt , , 


£RWl009j 
LK Wl 004 l 
tkWl 0092 
fcKWloOVi 
t KWl C094 
EKW10095 
CRW1P046 
t KWl 0047 
tk.wioova 
fckwi CG99 
tRWloluO 

FkWlGlOl 
FkMl 0102 
FRW1 Ol 03 
ARW20001 
ARW2 0G02 
ARW20003 
"AFW20004 
AKW2G0G6 
AKW2 00 06 
ARW20007 
ARW200C8 
AH W20009 

ARW2 00HT 

AKW20C11 
,AKW20012 
ARW20013 ~ 

AKW20014 

ARW20015 

ARW2GC16 

ARW20C 17 
, ARM20016 
ARW20019 
AR M2 0020 

&KW2C0Z1 

*bkW2CG22 

7RW20023 

bkW2G024 



t JAM 75 ) 


MAIN 


oS/3<>v RKlKAN H IXTtNOlL* 


DAT t 7o« 


ISN 0020 





ISM 0021 

M 


a-.oogk,? * -.fin#, -..c lit, . l.* , , hKw2.025 

, -,MC )oS , -.'v!Cb , -.JCtOv , -.''1X10 2 , Oil , fKW2'jw2o 

f-.ooon*, ,-.oooiib ♦ -.ooc-i? , -.«:c<>i 2 #-.*x~i 2 , .* 2 t*> , t>w2-*Cv7 

7.J0 72b , .0002 3 , .00^2.6 , .*X14 , .COCK' , .'.001? , U<W2o028 

0.07016 , .03016 * .02- -It* , .CliCfc , .0' r '4 , .7C'C1 , -.000)25 , 6RW29C29 

9— .{.'>30«* , — .uC0."5 , — . .rj'»Ot> , —«»/“■ 07? , — . GuG07fe , -.jO'GTb , bKW2003C* 

1-. 07008 , -.000125 , — .OOi lt» , -.tO'-lUi , -.C 702 » -.Cu02vib , hK#2O03l 

1- . 0*021 , -.000215 f -.OGOri* , — . 0- r, C2 1 ? , -.<%0&7fc , -..'0012 5 , &KW20G'A* 

2- .O/Clo , —.000180 y — »wiw2 f — ■ WOZi}* 1 , — »vv721 , — .000215 , blW20di3 

3- .0732I* , -.OC2212 , .CCutu* , .OC-'SGU » »j75C , .*‘O0<>9 , 8«W2l>034 

4. 02265 # .CO 0(> 1 » .3005b , .:33 5b , .07056 ,.Xt5b f .t: .-5bb , bKW2G&35 

5.000507 t .000425 * * .007295 * .ui'JJo , ,u0024 , .00624# bKK20vJ3fc 

b. C 0024 , .C0C24 , .00035 , .C" 02*5 , .C' Olt. » .Of C76 , .0000 3 , bR*20037 

7-. 03031 , -.00003 t -.00004 , -.00004 ,-.00004 , 0.0 *-.00013 ♦ 6kW2003b 

e-.0032? « -.U00287 , -.0JO312, -.00-3338 # -.00035 , -.000362 , bkW2C-039 

9-. 37 337 ,-.000369 / bRW2v040 

DATA CRW2 / .30015-, , . ' 3* K , .003123 , .0301 , .COCO 8 , CKVCOOhI 

1 .00007 , .030 62 , .000055 , .00t<»51 , .000051 , -.0001 , CRW20042 

2- .C33128 , -.003149 , -.CtClo , -.00017 , -.000172 , -.000176 , CRW20043 

3- .OVH6 ,-.00016 #-.00016 , -.0.01 * -.00012b * -.000149 , CRW20044 

4- . 00016 , -.00017 « -.000172 ,-.007 17b , -.OOClt , -.07C1P, , OKW200-»5 

5- .0;0l8 , -.0001 , -.00128 , -.(a.0149 , -.00016 , -.00017 , CKW20046 

6- .C33172 , -.C73178 , -.07018 , -.00018 , -.0)016 » .70948 , CRW20047 

7.00046 #". 00042b , .0^37^ , .00346 , .CO 7 32 , .20324 , .2C02B , Ck*20B48 
8.0C328 , .C-7028 , ."0017 , .00012 , .3X07 , .CC0012 , -.0^002 » CR»C0044 
9— .0707 49 , — .©0'*07 , -,0000b , -.000078 , -.000u7e ,“*00008 # CHW2G05" 

1-. 00014 , -.0701«» , -.000232 , -.00)256 , -.00027 , •*. 00 3286 , CHW20051 

1- .031295 , -.000295 , -.CGC29 , -#C‘.0C8 , -. 23614 , -.02019 , CRW2O052 

2- .00023? , -.000254 , -.00027 , -,C702bb .-.*00295 t-.00>295 , CHW20G53 

3- . 00024 f .001214 , ,.001163 , .0"U5 , .00111 , .001092 tftW20054 

4.00107 , .00137 , .00107 , .C/71o7 ,.001-77 , .00082 , .0007*. , CKW2 0055 

5.00066 , .00058 , .00054 , #0. ,i j9 , .OC’C'49 , .00049 ,.00049 , CRW2G056 

6.00349 , .073495 , .-33C186 , .020262 , .C0C175 , .0uul?6 » CRW20057 ‘ 

7.007075 , .000032 , P. . ,0. , 0. , .00002 . -.03011 , CRK20O58 

h-. 00022 , -.C30325 , -.00037 , -.CO041 , -,C'G43R , -.30045 , CRW20059 

9-.0O346 » -.030469 / CRW2OO60 

OATA 09W2 / .070155 , .00C13o , .00012 , .000092 , .C3Jv35 , 0RW2O061 

1.30-3 370 # .CC0C6G , #CCv042 , .CO. >023 ,,300'y23 , — >0v76 , DRW2uOt2 



( JAN 75 > 


MAIN 


ilS/iov FORTRAN H iXU-NDtl 


HATE 7c. 162/15. 04.43 


S 


c 

ISN 0022 


C 


2- .O301 -.000173 « -.rCTl A f -.r*C)A7 . -.•"f'15 , -.(*07154 , 

3- .G.70155 , — • Of : 1 •? 5 , , -.COO 75 , -.00>1 , -.100123 , 

4- .C731A , — .0 >01 a 7 , -.<i3ul5 , 0>154 , -,>00155 * -»*i0ul55» 

5- .030154 t -.00007*i , -.0C01 ,- .CO*' 123 , -. O' CIA -,?0 j147 « 

6- . 07*15 f -.00015 a , -.070 155 , -.0^0155 , -.00t.lt A , .tCCA’5 t 
7.000 38 * .0003a , .0003 * .70028 , .0CC2o ♦ ..‘0*2* * t .0'«722 * 

8 .003 2 2 » »O r 022 , .GOvlJO , .0C“SC7 » .000030 » -.C*j002o , 
V-.073042 t -*00006 3 t -.*>**0.80 , -.»OCl * -. '00115 * .".'JK3 11 5 , 

1-.00008 , -.030150 , 00237 , — .0 CC2A5 , -.CC':2 O t ,'j2 7a 9 

1- . 000280 * -.000284 , -.00078* .-.OOCJfcO , -.OOOO* , -.000150 , 

2- .00T207 f -.C03 ?a 5 , -.‘007o , -,C/ 01*7 a , -,i" **-02hC , -.OouJBa « 

3- .00028 5 * -.000280 , .001008 » .00103 , .10097 , .0<70v>5 « 

4.00087 ,*.07084 , .030815 , .v 00 El 5 , .700C15 , .000815 ,.00374 , 
5,007638 , .C 7*53 , .'..C^2 , .C 3036 , #jC0315 , #yuu31 , «oo032 i 
6.00732 , .00032 , ,0K4a , *0 ! ‘ 31 , .0:>lt* , .."Ojc>5 , ,j9C305 < 

7- .C3004 , — • OGQC 7 , -.CO-" 1 07 , -.OOOOc » -.00006 , -.03002 , 

fi— .000175 , — . 000275 , —.000385 , —.000425 .-•OOua. , — . 0 0 ‘i A 8 , 
*>—•00044 , -.00044 , -.CO -4b 5 / 

DATA ERW2 / .000102 , .OCOll , .070114 , . C 1 1 , .COOlw 5 , 

1.^071 y .0Cf)09*> , .PO0084 , ,.7070 <;0 , — .1. vOOao , 

2- . 07078; , -.000055 , -.0001 , -.00011 , -.00012 » -.00013 , 

3- .000137 , -.00014A , -,C0015t! , —.00043 , — .*i&CGt.*' , — .uu>‘.»85 , 

4- . 0071 , -.00011 , - • J C 7 1 2 , -.c ::i3 , —•"'30137 , -.70C1-A , 

5- . 0*0158 , -.030040 , -.070360 , -.370085 , -.0701 , -.00011"', 

6- •0 , "'0l2 , — . 0001 3 , -*(1*"'137 , —,000 1 4a ,~i',i071M> , *070472 t 

7.000433 , .00 >40, .000375 , .00 34.7 , .000 3c. , .*" 037 , .700365 , 
8.003365 , .000385 , .00 'll , .00005 , —.007*1 , —,*.0005 , 

9— .0. .‘OCA 8 , — . 007 8 , — .OGC09 , — • U >01 , —.00011 ,— .10011 , 

l-.0001t>fe , -.00 •.'225 , -.0302 7 , 0u3i , -.00C 329 , -.C0034 , 

1— .00035 , -.00036 , -.'00036 , -. C70357 , -, :001oE , -.0CC225 , 

2— . 00027 , -.00031 , -.000329 , -.T** 3a , -.C7035 ,-..*0036 , 

3— .03036 , -. 0 30352 » ,0015a , . 0>l46E , .* CIa , 01338 , 

4.001305 , .07128 , .00128 , .07128 ,.C;i2B ,.CC123 , .001025 , 

5.0004 , .000760 ,' .30066 , »9C*Jo35 , .'iOC't»i , * 0 * 0 o * ».C076*> , 
6.00065 ,.‘*0065 , .'"''062 y .tOOAo , .."0.33?; , ,00' 1 a , .707 13 , 
7.000375 , .00-7755 , .00.775 , .007075 ,.000075 , ..jO 706 , 

8- .000115, —.000255 , — .vO’736 , —.*>00435 » — . jOuab , —.*>0051 , 

9- . 00-753 , -.0C053 , -.00052 / 


DRW2C063 

DRW2006A 

DR W 2 0065 
OR W 2 0066 
DRM 200 O 7 
DPW 2 J 066 
DR W 2 0069 
DRW 2007 ^ 
(KM 2 C 071 
UR M 2 0072 
DKW 20073 
DR M 2 C 07 A 

DRW20075 

OH W2 00 76 
LRW20077 
DRW2007R 
URW2 0074 
OR M2 0080 

6 RM 2 0 C 61 
LRM 20 C 62 
ER M 2 0083 
E* M 2 0084 
ERM 20065 

FRW70086 

ERM 20087 
EKM 2 C 088 
ER M 2 0089 
tKW 2 o 09 o 
EKW2GCVI 

FRW200V2 

fcRW 20 C 43 
ERW 20044 
ERW 20095 " 

E RW 2 0096 
EKM 20097 
ERW 200 W*- ~ 
6 RM 20 C 99 
6 RMZ 010 U 



I JAN 75 » 


MAIN 


is/7t»‘ 


FOklKAN li iXn.NUtL) 


OAlfc 


ISN 

0023 

DAT* FRV»2 / • j i 3v?26 * • '*<•* 2> 2 t #000175 , • i’COI 4** , «uO«l?£ 

1 *000 1 1 t *U0 OC'4 , *0)0*4 » • C >1 .>? ,. *v» '1 2 » — 0 Gt* 7 , 

2-. 33010 3 » -.07:132 , , -.‘CC17; , -*.uK: , -. 'Clk. 

0 

o c 

to 

3-.o>nft7 t •••jQ'j it? , -..c m , -.:■:?■*#? « — • :-o • 1 . l - * .*;;>13? 

A.-.ooni6t , - .000170 , -.c'civb , -.cooifc.' » -.oo^no , -.jj-.it; 

s-. o)oi8i * -.oo'-c©? • -.or.nios v -.*.7Ci3* , -.coci*,* , -.conn 

6— .000175 , — .077180 , — J0lfc2 » — .vC01£2 , — .iGClGl * .*>w067> * 

7.20354 , .07)49 » .00 343 , .CC03 92 * .000455 , .0> r * 36 • .lev*.! » 



8.00041 ,.00041 , .00023 , .TO -lc ♦ .00*86 , .'„-*;?716 , -025 



9-.000.l6 t -.00009 t -.0*' ' OVO,-. 00 O.oO , -.vCCiCfcv , -.030126 , 
1-. 0)0167 , -.00021 , -.C00?6 , -.000287 , -.000316 , -.3003*. , 



1-.W036 , —.07034* , —.00033 ,—.0.0125 , — .uGOl**/ , — .00021 , 



2- . 37026 , -.0)3287 , -.CO . 316 , -. 00024 , -.."C03o ♦ -.D.'.C36 , 

3— .00033 , .00163 , .00166 , .00l4c5 , .7013u , .0713.' , .001242 
4.00127 ,.00127 , .07127 ,.C 127 , .00117 , .C'103 , .OOOhH , 
5.00073 , .00065 , .000575 , .000565 , .00u&2 ,.07062 ,.300«*2 , 
6.00778 , .000545 , .6042 , .00024 , .000165 , .00008-7 , 



7.007057 , -0??08C , .00712 ,.00012 , .C.C2 , -.UC'oOlS • -.6002 , 
B— .07 735 , -.005414*, -.07 *'47 , -.00)616 , -.00064 . 



4-.0005*. , -.00052 / 

ISN 

0024 

PI«3. 141543 

ISN 

0023 

FTM « .3048 

ISN 

0026 

F2TM2 « .092903 

ISN 

0027 

MTF * 3.28064 

tSN 

0028 

M2TF2 *10.76391 

ISN 

0024 

S1-1C090. 

ISN 

0030 

1 READ <5,101 AC, AMAX , AN, XM AX,L,SREF,M 

ISN 

0031 

10 FORMAT! 7F 1C. 6) 

ISN 

0032 

IF (AC I 5,7,7 

ISN 

0033 

5 SREF«SR6F*M2TF2 

ISN 

0034 

7 COF12*ft. 

ISN 

0033 

CtJWl 2*0. 

ISN 

0036 

00012*0. 

ISN 

0037 

CI>F232=C. 

ISN 

0038 

CDW232=C. 

ISN 

0039 

CD0232**‘». 

ISN 

0040 

CD0M*0. 

ISN 

0041 

SR AT 10 *Sl/SRtF 

ISN 

0042 

XM1)*XMAX/L 

ISN 

0043 

Xf (2 )*AN/AbS (AC 1 


O' 


FRW20101 
^ W2 01 v2 
FKWj civ 3 
»'Fkw2 Olv** 
»FKW2010t 
FRW2 01 0c 
FRW7D107 
F k W2 0 1 u fc 
FKW20109 
FRw2 0ll-5 
FRW2 0111 
FKW2 7112 
FKW2 01 a 3 
FR W20 1 1 4 
FRW20118 
FkW2011t> 
FKW20117 
FRW201 1 fa 
FKW20119 
Fk *2 120 



§ 
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MAIN 


US /So.' FORTRAN H fcXUNOtD 


DATE 7o. l62/15.U*.43 


? SW 0044 
ISN 0045 
ISN 0040 
ISN 0047 
ISN 0048 
ISN 0049 
ISN 0050 
ISN 0051 
ISN 0052 
ISN 0053 
ISN 0054 
ISN 0055 
ISN 0056 
ISN 0057 
ISN 0058 
ISN 0059 
ISN 00 60 
ISN 0061 
ISN 0062 
ISN 0063 
ISN 006* 
ISN 0065 
ISN 0066 
ISN 0067 
ISN 0068 
ISN 0070 
ISN 0071 
ISN 0072 
ISN 0074 
ISN 0075 
ISN 0076 
ISN 0078 
ISN 0079 
ISN 0080 
ISN 0082 
ISN 0083 
ISN 0084 
ISN 0086 
ISN 0087 
ISN 0088 


XG(3 ) * AMAX/AB 5 ( AC ) 

0C*S0RT(4*A8SIAC )/ Rl) 

xGf4)*t./nc 
IF 4 AC I 15,17,17 
15 AC.M*AbS(A() 

AC*ACM * M2TF2 
AC1-AC 
AMAXM* AMAX 
AMAX*AMAX * M2TF2 
ANM sAN 

AN *AN * M2TF2 
XMAXM *XMAX 
XMAX *XMAX * MTF 
LM *L 
L*L * M1F 
SRCPM*SREF*F2TM2 
GO TO 20 
17 ACM*AC6F2TM2 

AMAXM sAMAX ♦ F2TM? 

ANM * AN * F2TM2 
XMAXM* XMAX * FT M 
LM* L* FTM 

SRFFM * SP.fcF * F2TM2 
20 XG (5)*AbS(AC ) 

IF iXG(l>.LT..4 .OR. XG(1) .CT.l. )GU Tt 2*- 
GO TO 30 

25 WRITE 16,101) XG<1) 

30 IF (XM2I.LT. 1. .OR. XG(2).CT.2. ) CO TO 35 
GO TO * 3 

35 WRITE (6,102) XG(2) 

40 IF ( XGI 3 ) .LT . 1 .2 5 .OR. XG( 3) .GT.2. I GO TO *5 
GO TO 5C 

45 WRITE (6,103) XG(3) 

5v IF (XG(4).LT .5.5 .OR. XG(4).GT.7.) GO TU 55 
GO TO 60 

5* WRITE (6,104) XG (4 ) 

60 IF (XG(5) .LT.20. .OR. XG( 5 I.GT .4 0 . ) GO TO 70 

65 GO TO 75 

70 WR ITt (6,105) XG( 5) 

75 XG(6)*t.i 





( JAN 75 ) 


main 


US/3t>s, FUKIRaN H IXHNOH) 


DAT £ *?b. It>2/i5.09. A3 


ISN 0089 

ISN 0090 
ISN 0091 

' ISN 3092 
ISN 0093 

ISN 009A 
ISN 0095 
ISN 0096 

ISN 0097 
ISN 0098 
ISN 0099 
ISN 0100 
ISN 0101 

fSRnno? 

ISN 0103 
ISN 0109 
ISN 0105 
ISN 0106 
ISN 0107 
~ 15RT0138 
ISN 0109 
ISN 0110 
ISN 0111 
ISN 0112 

~ ISN 0T13 


ISN 0119 
ISN 0115 
ISN Olio 
ISN 0117 


ISN 0118 


C 

C 


C 


C 


C 


... CDW AT 1.2M 

CALL NDTLALI6 ,Xf>, X , Z , NX ,U>W 1? * A ,b » t I 
... crw AT 2.32M 
XG(t>)* 2.32 

CALL NDTLaFU ,XG,X,7,NX,CDW232, A,C.,C,t»l 
... COF AT 1.2M 
XG(61*1.2 

CALL NDTLAH I 9 ,XC1 ,X I , Zl , NXA ,C PF12 ,A , b,C ,1 I 
... COF AT 2.32 M 
XG (6 1*2.32 

CALL NDTLAEI9 ,XG1 ,X1 , 7 1 ,NX1 ,CDF 23 2, A ,B, C, 0 I 
CDF12*Cl)F12 * SR AT 10 

CPW12«CCWr2 * SRATIO 

COF232*CDF232 * SRATIO 
CDW232* CDW232 * SRATIO 
C0012*CDFI2 ♦ C0W12 
CD0232*COF232 ♦ CDW232 
• . • COO AT INPUT MACH NO. 

IF fM I 63,95,83 * " 

83 OER*ICD0232— CD012I/1.12 
S( ?1*2.*0ER 
S<3)>.35* DER 
LIPN*2 
IP*l 

zpRin*coonr — — 

2PR(2 )*CD0232 




XPRI 11*1.2 
XPR( 2 I =2.32 

A = -I2.*IZPR(LIPN! - ZPRILIPN - 1 >1 - 4SI3I ♦ SI 2 1 1 *IX Pft I IP NDT4»v28l 

2* II- XPR (IP 1 1 I/I XPR 4 IP ♦ II - X PR 4 IP 11**3 NDT00282 

B * CIIZPRILIPMr "=• IPRTLIPN - in / fXPRCTP * II - XPRIir IINtJTtn>2&T' - 

21 - SI2I- (XPR 41 P ♦ II ♦ 2.*XPR 41 P ll*IXPR41P ♦ II - XPR4 IP ll*AN0T00289 
31/IXPR I IP ♦ II - XPR4IP II N0T002B5 

C * SI 2 1 - 2,*XPR(IP l*b -3.WRUP |**2* A NDT00286 

0 * ZPRILIPN - II - - XPRI IP l*(C ♦ XPR I IP 1*18 ♦ XPRIIP i*Al I NUT30287 
CDOM«D*M*(C*M*<b+M*A) ) 

95 WRITE (6,100! ACH ,AC“7~XMAXr ANA X , ANM ,AN , XMAXM , XMAX , 

1LM , L , SRfcFM ,SREF , C0F1 2 »CDW12 ,CD01? , C0F232 ,CUW23? , 

2C00232 

100 FORMAT 1 1 HI 20X, 16MNACFLLE GEOMETRY //1H 2oX, 19HC.APTURE AREA 



o 
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MAIN 


US/;>6.' FORTRAN H EXTENDED 


DATE To, 162/15.04.43 


1SN 0119 
ISN 0120 
ISN 3121 
ISN 012' 
ISN 0123 

ISN 0124 
ISN 0123 


ISN 0126 

ISN 0127 
ISN 0128 
ISN 0129 
ISN 0130 
ISN 3131 
ISN 0132 
ISN 0133 
ISN 0134 


ISN 0135 
ISN 0136 
IS N 0137 
ISN 3138 
ISN 0139 
ISN 0160 
ISN 0141 


1 Fo,2 , 2X » ) 2H Su M t FG , 2 X, fcH SO FT),/ lh 2oX ♦ lkhMAX IMUM AREA 

2 F6.2, 2X ,1 2HSC M ( Fc> .2 ,2X , fcHSP FT ) / 1H 2oX» 19hNU2ZLt 

3AK IA F6.2,2X,1?HS(. M ( Fo.2,2X,6hSC HI/ lh 2oX,lRHLO 

4L. OF MAX. ARE A Ffe.2 ,*»X, 1 PhM ( Fu.2 ,-X, *»hFT >/ lh 2oX,19HT 

SOT AL LENGTH Ft.. 2 *4X, 1 tiHM ( F6.2 ,*X ,4HFT )//lH 2<‘X,2? HW 

toING REFERENCE ARfcA FV.2,2X,VMSt M ( F9.2,?X,oHSL FTI//1H 20X,3 
77H INCREMENTAL NACELLL DRAG CUE FF 1C It NTS/ In i6X , 17F.MACH 1.2 CDF 
6 = Fb.5*3X,5HCDW= Fb.* ,3X, *'HCU = Ffc.5/ lh 2oX»17FMACh 2.32 CDF* 
4F6.5.3X ,5HC0W= F b. 5 , 3X « *.HCDw»* FF .0 ) 

IF <M) 1P9, 112,10 

109 WRITE *6,110) M.CDCM 

110 FORMAT * 1H 26X,4FlMACh F5 . 2, 35 X, 5HC0 0 = Fb.5) 

112 WRITE *6,115) XG(l),Xr.(2),XC(3) ,XC*a) 

115 FORMAT */lH 2CX, FHXMAX/L* F7.3,3X, ThAN/A C* F7.3.3X, VHAMAX/AC* F 
17.3.3X, 6HL/DC* F7.3) 

WR 1TE (6 , 140 ) 

140 FORMAT *//lH 2CX .ElHlNCREMENTA L NACELLE DRAG CcLFF'ICUNTS - LTV Rfc 
1FERFNCF VEHICLF/lh 2c>X , 25FiMACh 1.2 CtF* ^ . 0CC72 ,3X, ISfiLDW* -0. 
200041, 3X,!3hCOO* C.00031/IH 2oX,25«»MACH 2.32 tUF* 0 . .• ) ‘tbZ , 3X , l 
33HCDW* — 0 .CT018, 3X, 13 mCD0* C . 0 J 2 ) 

IF <M) 201,310,201 

C LTV REFERENCE VEHICLE NACELLE WAVE AND FRICTION DRAGS 
201 IF (M-1.4) 204,234,20t> 

204 1*1 

GO TO ?20 

206 IF *M-1.8I 206,238,210 

208 T *3 

GO TO 223 
210 1=5 

220 LTVW* I H—MLT * I + l ) )* *M-MLT * 1 »2 ) ) /* *MLT ( 1)-M’ T < 1+ 1 ) )* (MLT ( 1 )— MLT *1 + 2) 
1) )*LW* I ) ♦*M-MLT* 1 ) )**M-MLT ( 1 + 2) )/ **MLT* 1*1 )-MLT*l ))*tMLl * I >1 )-*MlT 
2(1*21) )*LW(1+1) ♦ ( M—MLT ( 1 ) )**H- MLT* I ♦!))/» I MLT 1 1 *2 )-MLl » 1 ))*IMLT 
3 * I +2 ) -MLT 1 1+ 1 ) ) ) ♦*. W ( 1 ♦? ) 

LTVO*LTVW 

230 WRITE <o,145) M.LTVO 

145 FORMAT * 1H 2oX,4hMACh F5.2,33X,7H CUC* Ft. 5) 

131 FORMAT <1H1 2CX, 22HXMAX/L OUTSIDE RANGE *FO .2 ) 

132 FORMAT *1H1 2CX, 21HAN/AC OUTSIDE RaNCF =F6.2) 

103 FORMAT *lhl 20X, 23HAMAX/AC OUTSIDE RANufc *Ft>,2) 

104 FORMAT *1H1 2UX, 20HL/DC OUTSIDE RANGE *Fo.< ) 
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DATt 7o. X t»2/ 15 .o*.*3 


ISN 01*2 
1SN 01*3 
ISN 01 -»* 


105 FORMAT 1 1H1 2 jX t ItHAL UJTSiOl P /.NO 
310 GO 10 1 
fcNO 


O 

'll 


b 

7C 

© 

Cl 

r 

►—» 

K 


§ 

O 

2 


o 

cd 

a 





( JAN 7*1 > 


0S/360 FORTRAN H F XT ENDED 


DATE 76.009/08.56.13 


RE0UFS1ED OPTIONS: L 1ST ,XRF F ,OFCK 

OPTIONS IN EFFECT: NAME(MAIN) OPTIMIZE!!) L I NfCDUNT ( 4? ) SIZE(MAX) AUTODBLINPNE I 

SOliRCE EBCDIC LIST DFC* OBJECT MAP NOFORMAT GOSTMT XREF ALC NOANSF NOTERMINAL FLAG! I) 


ISN 0002 SUBROUTINE NOTLAF ( NI V* XG , X, Zt NX* RES* A* 6t C. D) N0T00010 

C NDT00020 

C NDT00030 

C NOMENCLATURE NDT00040 

C NIV * NUMBER OF INDEPENDENT VARIABLES N0T00050 

C XG = ARRAY OF GIVEN VALUES OF INDEPENDENT VARIABLES N0T00060 

C X = ARRAY OF INDEPENDENT VARIABLES NDT00070 ......... 

C 7 = ARRAY OF DEPENDENT VARIABLES NDT00080 

C MX = ARRAY OF NUMBER OF VALUFS GIVEN FOR EACH NDT00090 

C INPf PENOENT VARIABLE NDT00100 

C RES * FINAL INTERPOLATED VALUF OF Z AT XI GIVEN* X2 NDTOOilO 

C GIVEN, XNI V GIVEN NDT00L20 

C A*B * C * D = COEFFICIENTS OF THE CUBIC DEFINING THE Z VALUESN0T00130 

C TCF * ARRAY OF NUMBER OF POINTS TO BE USED FOR CURVE NDT00160 

C. FIT FOR FACH INDEPENDENT VARIABLE N0T00I50 

C XXG * ARRAY OF GIVEN VALUFS OF INDEPENDENT VARIABLES NDT00160 

C UNLESS VALUF WAS OUTSIDE RANGE. IN THIS CASE NDT00170 

C CLOSEST VALUE IN THE TABLE IS ASSIGNED XXG NDT00180 

C IGAT * ARRAY OF INDICATORS OF LOCATIONS OF X GIVEN IN NDT0Q190 

C X ARRAY NOT 00200 

C l* IF XG IS IN ONE OF THF MIDDLF INTERVALS N0T00210 

C 2* IF XG IS IN THE FIRST INTERVAL NDT00220 

C 3* IF XG IS IN THE LAST INTERVAL N0T00230 

C LUC * ARRAY OF LOCATION IN X ARRAY OF FIRST VALUE TO NDT00240 

C BE USED FOR EACH INDEPENDENT VARIABLE . NDT002S0 

C XPR * TEMPORARY ARRAY FOR X VALUES BEING USED NDT00260 

C ZPR * TEMPORARY ARRAY FOR Z VALUFS BFING USFD NDT00270 

C LISTZ * LOCATION OF FIRST VALUE USED IN THE Z ARRAY NDT00280 

C LLCTR * COUNTERS FOR THE ZS WITHIN A SUBSET NDT00290 

C L SUBSCRIPT FOR Z -PR IMF NDT00300 

C LL * SUBSCRIPT FOR Z NDT00310 

C * SUBSCRIPT INDICATING THF SUBSET NDT00320 

C NDT 00330 

C NDT 00 340 

C DIMFNS1GN NX(l) f L0CIN1VI, XG(1»,X(1), Z(l), Z PR 14 ♦ 3*(NIV - 2>>,NDT00350 



OS/360 FORTRAN H EXTENDED 


DATE 76 .009/06. 56. 13 


ISN 0003 


„ISN 0004 


ISN 000*> 


ISN 0006 
ISN 0007 

ISN 0008 


XPR<4), S ( 3 ) » XXG(NIV), TGAT(NIV), ICFINIV), LLCTR(NIV+ 1) NDT00360 

NDT00370 

NDT00380 

OIMENSION NX(1), L0C<10>, XG <1) , X<1), Z<1), ZPR128), XPR(4), SI3INDT00390 

i, xxg ( in, igat(io), icf(ic>, llctrud ndtooaoo 

00 1 I * 1, N1V ..NOT 00410 

NDT00420 

TFST TO SFE IF FACH VARIABLE HAS MORE THAN ONE POINT NDT00430 

NOT 00 440 

TF(NXm-l) ?, 2, 3 NDT00450 

NDT00460 

... .VARIABLE, l DNLY HAS ONE POINT, THEREFORE PROGRAM CANNOT CONTINUE, , NDT0Q470- 

N0T00480 

2 WRITE ( 6»l r '0) I,NX(II NOT 00490 

IOC FORMAT ( 13H1 VARIABLE 13,10H ONLYHAS I2,41H POINT, THEREFORE RRNDI00500_ 

20GRAM CANNOT CONTINUE) N0TC0510 

STOP NDT00520 

... MOT 00530— 

SFT UP ARRAY ICF WHICH IS THE NUMBER OF POINTS REQUIRED FOR CURVE NDT00540 

FIT N0T00550 

_ NOT 00560 - 

3 IFtNXin - 3) 4, 4, 5 NDT00570 

4 ICFf 1 1 * NX(I) NDT00580 


3 IF(NXm - 3) 4, 4, 5 

4 ICFI I 1 * NX(I) 

*>~ICF( I ) * 4 
1 CONTINUE 

ICFfNIV ♦ 1) IS ARTIFICIAL VALUE TO CONCLUDE CURVE FITTING FOR 
ENTIRE ROUTINE 


NDT00600 
NDT00610 
NO TOO 620 
N0T00630 
NOT 00640 


r c 


ISN 0014 
ISN 001*5 


ISN 0016 
ISN 0017 
ISN 0018 
ISN 0019 


ICF (NT V ♦ 1 ) = 1 NOT 00 660 

IT * 0 NDTOO n 

, NDTOLadC 

FINO LOCATIONS OF INTERVALS CONTAINING GIVEN INDEPENDENT VARI ABLE SNDT 00690 
THUS FORMING THE LtJC ARRAY. NDT00700 

... ... . _ NDTOO 710 

DO 6 I = 1, NIV NDT00720 

XXG(I) = XG ( I ) NDTOO 730 

IL ■ IT ♦ 1 _ NDTOO 740 

I IX = IL NDTOO 750 


v| 

04 
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ISN 

0020 



IT » IT ♦ NXJ II 

NDT00760 

ISN 

0021 



IC.ATU ) « 1 

NDT00770 

ISN 

0022 


1 4 

IEUJX - IT) 7 , ?, e 

NDT00780 



c 



NDT00790 



r 


ERROR OCCURRED 

NOT 00800 



c 



NDT00810 

ISN 

002* 


o 

IERR ■ 14 

N0T00620 

ISN 

0024 



go to Aroc 

NDT00830 



c 



N0T00840 



c 


C0MP\RF GIVEN VALUE WITH X ARRAY FOR A PARTICULAR INDEPENDENT 

NDT00850 



c 


VARIABLE 

N0T00860 



c 



NDT00870 

ISN 

0025 


7 

IF<XXGtn - X f I IX » ) 12 t 11 , 10 

NDT00880 

ISN 

0026 


1 0 

IF ( XXG ( I ) - X(TIX ♦1)1 14, 14, 13 

NDT00840 

ISN 

0027 


1* 

IFIIIX ♦ 1 - 1T1 13, 16, 17 

NDT00900 

ISN 

0028 


17 

IFRR ■ 13 

NDT009I0 



c 



N0T00920 



c 


ERROR OCCURRED 

NOT 00930 



c 



NDT00440 

ISN 

0024 



GO TO 5000 

NOT 00950 

ISN 

00*0 


15 

IIX « IIX ♦ 1 

N0T00960 

ISN 

OJjI 



GO TO 10 

NOT 00970 



c 



NDT00980 



c 



NOT 00990 ... ... 

ISN 

i 032 


1 6 

GO TO 20 

NOT 01 000 

ISN 

0033 


1 4 

IFIIIX - IL) 18, 11, 20 

NDT01010 

ISN 

0034 


1 8 

tERR * 14 

NDT0I020 



c 



NOT 01 030 



c 


ERROR OCCURRED 

NDT01040 

* 


c 



N0T01050 _ 

ISN 

0033 



GO TO 5000 

NDT0I060 

ISN 

0036 


1 2 

GO TO 11 

N0T01070 

ISN 

0037 


11 

TGATU) * 2 

NDT0IO8O 



f 



NDT01090 



c 


X GIVEN IN FIRST 1NTERWI 

NDTOUOO 



c 



NDT01110 . 

ISN 

0038 


8 

I T X * IIX ♦ I 

N0T0I120 

ISN 

0034 



GO TO 23 

NDT0I130 

JSN 

0040 


20 

IFINXII1 - *) 8, 24» 24 

N0TOJ 140 

ISN 

0041 


24 

TFUIX ♦ 1 - IT) ?*, 25, 26 

NDT0U5D 



15 > 0042 


>N 004? 
($N 0044 


( JAN 


1SN 0Q4 ? 
ISN 0046 
INN 0047 
ISN 0048 
ISN 004® 


ISN 0050 


ISN 0051 
ISN 0052 
ISN 005? 
ISN 0044 
ISN 0055 
ISN 0056 
ISN 0057 
ISN 0058 


ISN 005® 
ISN 0060 


1T.N 0061 


74 ) NOTlAl 

2t I ERR * 24 

ERROR UCfURRFP 

GO TO 5000 
2 * I GAT ( I ) * 3 


0S/360 FORTRAN H EXTENDED 


C 

C 

c 

c 

c 


0 

c 

c 

c 

c 

c 


c 

c 

c 


c 

c 

c 


X GIVEN IN LAST INTERVAL 

rjeTFRMTNF LOCATION OF C IPST VARIABLE TO BE USEO FOR 

23 LOCI!) * lix - 1L 
6 CONTINUE 

L’ST2 * LOCH ) 

33 TFINTV - 1) 30, 31 # 3? 

*0 IFMR * 33 

ERROR OCCURRED 

GO TO 3000 

F|NO LOCATION OF FIRST Z IN ARRAY 

3? 00 35 I * 2 » NIV 
IRROO * 1 
10 * I -l 
00 36 IXP • 1. IP 
3 6 JPROD * T PROO*NX I I XP > 

IPROO ■ I PROD*( LOC < I > - 1) 

3 5 L1STZ « LIST 2 ♦ IPROO 
31 NC * NIV ♦ 1 

INITIALIZE COUNTER TO BE USEO ?0R Z VALUES 

00 37 I • ? 9 NC. 

37 LLCTRU) • 0 

INITIALIZE SUBSCRIPTS FOR Z ANr ZPRIME 
L ■ 0 


CURVE FIT 


§5 

f s 

^ o 

C re- 

5 £ 


N0T01400 
NDT01410 
NOT 01 420. 
N0T01430 
NOT 01 440 
liDXjQl 4Sfl 
NOT 01 460 
NOT 01470 
NOT 01400 
NOT 01 440 
NDT01300 

MflTOmtA 


N0T01520 
NOVO 15 SO 
NOT 01 540 
NDT0I550 


DATE 76.004/08.56.13 


NDT0I160 
N0T01170 
NDTOlliO 
N0T0I150 
NDT0&200 
NDT01210 
NOT 01 220 
N0T01230 
N0T01240 
N0T01250 
NDT01260 
NOT 01 270... 
NDT01 200 
N0T01240 
N0T01300 
NDT01310 
NDT01320 
NO 101310.. 
NDT01340 
NOT 01 350 
NDT0136O 
N0T01370 
NDT01380 
. N0TQ13ML 


i n 
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ISN 006 2 



1. L * 0 



N0T01560 

I5N 0063 


Sf 

LL » LL * LIST* 



N0T01470 


C 





NDT015S0 


C 


SFT UP SUBSCRIPTS 

- K 

FOP VARIABLE, L FOR Z PRIME, AND LL FOR l 

N0T01590 


C 





NDT01600 

IS* 0064 



K * 1 



NDT0I610 

ISN 0065 



L • L ♦ 1 



N0T0I620 

ISN 0066 



IC « If F(1 1 



NOT 0)630 


C 





NDT01640 


C 


MOVE Z 70 ZPP1MP 



NOT 01650 


c 





N0T01660 

ISN 0067 



on 42 j * i, ic 



NOTO 1 670 

ISN 0068 



M * L ♦ J - 1 



N0T01680 

ISN 0066 



7 PR ( M ) « ? ( LL ) 



NOT 01 640 

ISN 0070 


4? 

LL - IL ♦ 1 



NDT01700 

ISN 0071 


5 5 

IX * 0 



NDT01 710 

ISN 007? 



TIX » 4 



N0T0I720 

ISN 0073 



IF(NXlK) - 3) 44, 

44, 

45 

NOT 01 73d 

ISN 0074 


44 

TIX ■ NX ( K ) 



NOT 01 740 

ISN 0075 


4 4 

IF(K - 1) 47, 48, 

49 


NCT01750 

ISN 0076 


47 

IERR a 4 r 



N0T01760 


c. 





NDT017T0 


<; 


EPROR OCCURPFO 



NDTO170O 


c 





NOT 01 790 

ISN 0077 



00 TO 5000 



N0T01800 

ISN 0076 


44 

J « K - 1 



NDT01810 


c 





NOT 01020 - 


c 


SFT UP SIJPSCRIPTS 

FOR 

X VALUES 

NOTO103O 


c 





NDTO104O 

ISN 0076 



00 60 J1 • 1, J 



NOT 01 050. 

ISN 0080 


6 0 

IX « IX ♦ NX(J1) 



NOT 01 060 

ISN 0081 


4f 

TX * IX ♦ LOC < K ) 



NOTO107O 


c 





NOTO100O 


c 


MOVE X TO XP«7MF 



NDTO109O 


r. 





NDT01400 

ISN 0082 



no 61 J « 1, I IX 



NDT01910 

ISN 0083 



TXP • IX ♦ J - 1 



NOT 01 920 

ISN 0084 


61 

xpp ( j > ■ xuxp) 



NDT01430 

ISN 0085 



A * 0*0 



NDTO1940 

ISM 0086 



8 *0,0 



NOT 01450 




( 

JAN 

7A 1 

1 NPTLAF nS/3«0 

FORTRAN H FXT ENDED 

DATE 76 .009/08.96.13 



C 




N0T01960 



C 


te.>t numbfr pf points 


N0T0I670 



C 




NDTOI90O 

ISN 

0087 



IFUIX - 2) 63, 64, 65 


NDTOl 990 

ISN 

0088 


6? 

C ■ 0* 0 


NDT02000 

1 f.N 

0080 



on TO 66 


N0T02010 

ISN 

oooo 


b4 

C * (ZPR(L+11 - ZPR(L) )/( XPP <21 

-XPR ( 1 ) ) 

NDT07020 

ISN 

0091 


66 

P « ZPR (L) - OXPR 111 


NOT 02030 

ISN 

0092 



0,0 TO 2000 


NPT02040 

ISN 

0093 


6 8 

TL * I OAT ( K ) 


NOT 02050 



C 




NDT02060 



C 


DFTERMJNF what interval X GIVFN 

IS IN 

NOT 02070 - 



c 




NDT02080 

ISN 

0094 



IFUL.LT, 1 .OR. TL.GT.?) G C Tf 

68 

NOT 02 090 

ISN 

0096 



TO TO (70, 71, 72), IL 


N0T02100 . . . 

ISN 

0097 


68 

WRITE (6,4000) 


NDT02110 

ISN 

0098 

<.rnn 

FORMAT (• IN VAL TO INDEX FOR 

COMPUTED GO TO •) 

N0T02120 

ISN 

0099 



STOP 


NOT 02130 



C 




NOT 02 140 



C 


X GIVEN IN FIRST INTERVAL 


NDT02150 



c 


ONLY FIRST THRFf POINTS WILL bF 

USED FOR CURVE FIT 

NOT 02 160 



c 




NDT02170 

ISN 

0100 


71 

IS ■ 3 


NDT02180 

ISN 

OlQl 



TPN * 1 


.... NOT 02190 

ISN 

0102 



TIN * L 


N0T02200 

ISN 

0103 



M « 2 


NOT 02210 

ISN 

0104 



IRX • 0 


NOT 02 220 

ISN 

010* 


80 

*N • UN ♦ TPX 


N0T02230 

ISN 

0106 



TP ■ JPN ♦ IRX 


N0T02240 


. ... 

c 




NOT (12210. 



c 


CALCULATE THF SLOPE OF STRAIGHT 

LINE 

N0T02 260 



c 



o e> 

N0T02270 

ISN 

007 



S(M) • ( ZPR ( IN ♦ 1) - ZPR ( IN ))/( XPR ( IP ♦ 1) - XPRdPM ^2 

N0T022S0 

ISN 

0108 



IC * 1 

’v O 

NOT 02 290 

ISN 

0109 



0,0 TO 1050 

& g? 

N0T02300 

. ...... 


c 




NOT 023 10 



c 


a r,IVEN IN LAST INTERVAL 

£ ^ 

NDT02320 



c 


ONLY THREF POINTS WILL 8E USED 

y 2: 

N0T02330 



c 



f ' '•* M 3 

t wT3 

NDT02340 

ISN 

0110 


72 

I PN » 1 

^ tf) 

NOT 02 330 
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-4 

CO 


ISN 

ISN 

ISN 

TSN 

ISN 


0111 
01 12 
0U3 

0114 

0115 


ISN 0116 
ISN 0117 

isn one 

ISN 0114 


ISN 0120 


ISN 0121 

ISN 0122 
ISN 0124 
ISN 0125 
ISN 0126 
ISN 0127 
ISN 0128 
ISN 0129 
ISN 0130 
ISN 0131 
ISN 0132 
ISN 0133 


ISN 0134 
ISN 0135 

ISN 0136 


C 

C 

t 

C 


I!N « L 
IS » 2 
M * 3 
1»X * 1 
GO TO RO 

X GJVFN IN MlDDi.F TNTFRVAL 
FOUR POINTS WILL RE USED FOR 


CURVE FIT 


70 


IS « 
TPN 
1 IN 

i G 


2 


* 2 


C 

C 

c 


c 

f 

c 


CALCULATF SLOPE USING PARABOLIC FIT 

1050 AP • I 1 2PR (JIN) - ZPRUIN ♦ 1))*<XPR!IPN) - XPRUIN ♦ 
21 1 IN) - ZPRUIN ♦ 2) ) * 1 XPR 1 1 PN ) -XPR ( 1 PN* 1 )) ) /! ( (; PR! 
3PP<IPN*1))**2)*1XPR!IPN) - XPR ( I PN ♦ 2)1 - UXPR1IPN) 
4PN ♦ 2 ) )**2) * (XPR1 IPN) - XPRUPN ♦ 1))) 

S(IS) * 2.*AP*XPPUPN ♦ 1) ♦ 1ZPR1IIN) - ZPRUIN ♦ 1) 
2*2 - XPRUPN ♦ 1)**2)*AP)/1XPR(IPN) - XPRUPN ♦ II) 

IF UG.LT.l .OR. IG.GT.3) GO TO 1052 
GO TO (90, 92.92), IG 

1052 WRITE(6.*r01 ) 

500) FORMAT! • INVALID INDEX FOR IG ') 

STOP 

92 IS a 3 
IPN * 2 
X1N ■ L + 1 
IG * 1 
TP * 2 
GO TO 1050 

CALCULATF COfFFIClENTS 


NDT02360 
N0T02370 
NOT 02380 
NDT02390 
NOT 02400 
NDT02410 
NOT 02420 
NDT02430 
NOT 02440 
1 NDT02450 

NDT02460 
NOT 02470 
NDT02480 
NDT02490 
NOT 02 300 
NOT 02910 
2)) - IZPR N0T02520 
IPN) )**2 - (XNOT 02530 
1*42 - ( XPR UNDT02940 
NDT02590 
- (XPR1 IPN)*NDT02960 
N0TC2970 
N0T02S80 
-N0I0255C 
NOT 02 600 
N0T02610 
NOT 02 620 
NOT 02 630 
NDT02640 
NOT 02650- 


9 0 L I PN » L ♦ IP 

A ■ -I 2.* 1 ZPR ( LTPN ) - ZPR1LIPN 
2 * !)- XPR 1 IP ) ) )/ f XPR I IP ♦ 1) 
R * 11 IZPR(LIPN) - ZPRILIPN 


NOT 02 660 
NOT 02670 
NOT 02480- 
NOT 02 690 
NO f 02 700 
NOT 02710. 


1)) - IS!: 

XPR! IP ) )**3 
1 )) / ! XPR ( IP 


N0T02720 

) ♦ S( 2))4(XPR UP NOT 02 730 

NOT 02740 
♦ 1) - XPR! IP ) )NOT 02750 



NOUAF 
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?.) - $1?)- ( X p R ( I P 2.*XPR(IP >)*(XPR(IP 4 

11 - XPR (IP 

) |*AN0T02 760 




3 ) / ( XPR i IP ♦ 1) - XPPI1P n 


NOT 02770 

JSN 

0137 


e * S<?> - ? . *XPR ( IP )*P - 3 . *XPR ( IP ) **2* A 


NDT02780 

1SN 

0138 


D * ZPR1LIPN - n - X PR ( IP )*<c 4 XPR ( I P »*(8 4 

XPR (IP )*A> ) 

N0T02790 



C 



NDT02800 



C 

CALCULATE Z 


N0T02810 _ ... 



C 



NDT02820 

ISN 

0139 

?OCf. 7 PR t L ) * r 4 XX C ( K ) * ( C ♦XXr.(K)*(P 4XXfi<K)*An 


NOT 02 830 

ISN 

0140 


K « K 4 1 


N0T02840 

ISN 

0141 


LLCTR(K) r LLfTRfKI 4 1 


NOT 02 850 



C 



NDT02660 



C 

TEST Z COUNTEP WITH NUMEEP CP POINTS REQUIRED FOR 

CURVE FIT 

NDT02870 ... 



c 



NOT 02 880 

ISN 

0142 


I F ( L LCTR ( K ) - ICPIKI) 300, 301 , 301 


NOT 02890 

ISN 

0143 


300 1 F (K - 2) 3 1 p t 3io, 318 


NDT 02«00 

ISN 

0144 


315 KK * K - 1 


NDT02910 



c 



N0T02920 


- w - 

c 

RE-INI TIAL IZE PREVIOUS COUNTERS 


NOT 02930 



c 



N0T02940 

ISN 

0145 


00 320 I « 2, KK 


N0T02950 

ISN 

0146 


320 LLCTR(I) * C 


NOT 02 960 ... 

ISN 

0147 


31P LL « o 


N0T02970 

ISN 

0148 


LP * 1 


NOT 02980 


- 

c 



NDT 02990 



f. 

FIND NFXT l TO PE USED FROM ARRAY 


NOT 03000 



c 



NDT030I0 

ISN 

0149 


DO 400 1*2, NIV 


NDT03020 

ISN 

0150 


IJ « I - 1 


N0T03030 

ISN 

0151 


00 410 J « 1 , TJ 


NDT03040 

ISN 

0152 


410 LP * LP*NX(J1 


NDT 03030. 

ISN 

0153 


LL * LL 4 LLCTPU)*LP 


NDT 03060 

ISN 

0184 


400 LP * 1 


N0T03070 

ISN 

0155 


00 TO 50 


NDT03080 

ISN 

0156 


301 IF (K - NIV) 500, 500, 6 00 


NOT03090 



c 



NDT03100 



c 

FIND SUBSCRIPT OF NEXT Z p P ELEMENT 
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